Kod Programu:
clear all
clc
h=0.01
x=0:h:5
y=1
k=length(x)
for i=2:k
y(i)=y(i-1)+h*(-x(i-1).^2+3*x(i-1)-8)
end
f=-1/3*x.^3+3/2*x.^2-8*x+1
plot(x,y,x,f)
grid on
xlabel('x')
ylabel('y')
for m=1:k
blad(m)=y(m)-f(m);
end
plot (x,blad)
-
1 2 3 4 5
Wartośd błędu dla kroku h = 0,01 -0,01002 -0,01003 -5e-005 0,01993 0,04992
Wartośd błędu dla kroku h = 0,1 -0,1017 -0,1033 -0,005 0,1933 0,4917
Wartośd błędu dla kroku h = 0,5 -0,5417 -0,5833 -0,125 0,8333 2,292
Wartośd funkcji F dla kroku h = 0,01 -5,843 -11,68 -18,5 -28,31 -43,12
Wartośd funkcji F dla kroku h = 0,1 -5,935 -11,77 -18,51 -28,14 -42,68
Wartośd funkcji F dla kroku h = 0,5 -6,375 -12,25 -18,63 -27,5 -40,88
Wartośd analityczna funkcji F -5,833 -11,67 -18,5 -28,33 -43,17
h = 0,01
Wykres błędy dla kroku h=0,01
Wykres błędy dla kroku h=0, 1
Wykres błędy dla kroku h=0, 1
h=0,5
Wykres błędy dla kroku h=0,5