Budownictwo 2006 Wydział Inżynierii Lądowej Politechniki Krakowskiej

saddam

Ma ktos może w/w programik w Octave?

Mam rozwiązać podane równanie różniczkowe (podany przedział t oraz skok h) i wynik porównać z biblioteką oraz rozwiązaniem dokładnym.

R2Rek

daj maila bo mam tu coś, a jak nie to znajdziesz w moim katalogu w szkole. Raczej pięknie zrobione.

Dobry

Mógłby ktoś powiedzieć, co jest z tym nie tak - funkcja biblioteczna zwraca zawsze 1:/

y0=1;
x=linspace(0,4,21);
n=length(x);
h=0.2;

function df=fun(x,y)
df=x^2-4*x+3;
endfunction

function f=dokl(x,y)
f=1/3*x^3-2*x^2+3*x+1;
endfunction

function rozw_biblioteczne=dod(y,x)
rozw_biblioteczne=x^2-4*x+3;
endfunction

function rungekutta=rungekutta(y0,x,n,h)
rungekutta=zeros(n,1);
rungekutta(1)=y0;
for i=2:1:n
k1=h*fun(x(i-1),rungekutta(i-1));
k2=h*fun(x(i-1)+1/3*h, rungekutta(i-1)+1/3*h*k1);
k3=h*fun(x(i-1)+2/3*h, rungekutta(i-1)+2/3*h*k2);
rungekutta(i)=rungekutta(i-1)+h*(1/4*k1+3/4*k3);
end
endfunction

for z=1:1:length(x)
rozw_dokladne(z)=dokl(x(z));
endfor

rozw_rungekutta=rungekutta(y0,x,n,h);
rozw_biblioteczne=lsode("x^2-4*x+3",1,x);

Rozwiazania=[rozw_dokladne' rozw_biblioteczne rozw_rungekutta]

plot(x,rozw_rungekutta, '-@;RK3;')
hold on
plot(x,dod(y,x), '-;Biblioteczne;')
hold on
plot(x,rozw_dokladne, '-;Dokladne;')
hold on
pause