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;
y=A\b
y =
0.0553
0.0996
0.1318
0.1514
0.1580
0.1514
0.1318
0.0996
0.0553Vad ä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.