% algorytm ekstrapolacyjny eulera clc;clear; n=100; h=0.1; x0=zeros(1,n); x1=zeros(1,n); x2=zeros(1,n); t=zeros(1,n); x1(1)=10; % warunek poczatkowy alfa=1 for i=1:n x1(i+1)=x1(i)*(1-h*alfa); t(i+1)=t(i)+h; end x0=10.*exp(-alfa*t); plot(t,x0,t,x1,'+') ____________________________________________ % algorytm ekstrapolacyjny eulera clc;clear; n=100; h=0.1; x0=zeros(1,n); x1=zeros(1,n); x2=zeros(1,n); t=zeros(1,n); x1(1)=10; % warunek poczatkowy x2(1)=10; alfa=1 for i=1:n x1(i+1)=x1(i)*(1-h*alfa); t(i+1)=t(i)+h; end x0=10.*exp(-alfa*t); plot(t,x0,t,x1,'+') % algorytm interpolacyjny eulera for i=1:n x2(i+1)=x2(i)/(1+h*alfa); end plot(t,x0,t,x1,'+',t,x2,'*') eps1=abs(x0-x1); eps2=abs(x0-x2); plot(t,eps1,'+',t,eps2,'*')