% 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,'*')
Wyszukiwarka
Podobne podstrony:
cennik system k2 kan 06 2013TI 01 09 06 T pl(2)cennik studnie 06 2013wyklad 10 09 06 2 komorka chemwyklad 10 09 06 2 komorka budowaxenon E83 od 09 06876 tabela handlowa s2 048 06 2013TI 99 09 06 B pl(1)FINAL EXAM ADMINISTRACJA I ROK 22 06 2013Referendum ws paktu bezpieczeństwa opóźnione (09 06 2009)cennik system kanal zewn pvc 06 2013pdm? 2016 09 06więcej podobnych podstron