Vi vill lösa följande randvärdesproblem y″(x) = 2y − 1 och och skriver följande program med resultatutskrift till höger.
a=0; b=1; N=9; h=(b-a)/(N+1);
A=zeros(N,N);
for k=1:N
A(k,k)=-2-2*h^2;
if (k<N)
A(k+1,k)=1;
A(k,k+1)=1;
end
end
b=-ones(N,1)*h^2;
u=A\b
u =
0.0381
0.0670
0.0873
0.0993
0.1032
0.0993
0.0873
0.0670
0.0381Vad är y″(0.5) ungefär enligt denna approximation?
Du ser bara uppgiftslydelsen
Svarsalternativ, facit, figurer och formler — och resten av kursens uppgifter — finns bakom ett gratis konto.