Forum grupy 6 WIL

Gość

function df=funkcja(x,y)
df=2/x*y+x^2*exp(x);
endfunction

function df=dot(y,x)
df=2/x*y+x^2*exp(x);
endfunction

function y=dokladne(t)
x=t;
for i=1:length(t)
y(i)=x(i)^2*(exp(x(i))-exp(1));
end
endfunction

function rungekutta=metoda_rungekuttaIII(y1,t,n,h)
rungekutta(1)=y1;
for i=2:n
k1=h*funkcja(t(i-1),rungekutta(i-1));
k2=h*funkcja(t(i-1)+1/3*h,rungekutta(i-1)+1/3*k1);
k3=h*funkcja(t(i-1)+2/3*h,rungekutta(i-1)-1/3*k1+k2);
k4=h*funkcja(t(i-1)+h,rungekutta(i-1)+k1-k2+k3);
rungekutta(i)=rungekutta(i-1)+1/8*(k1+3*k2+3*k3+k4);
end
endfunction

a=1;
b=2;
h=0.1;
t=[a:h:b];
n=length(t);
y1=0;
rungekutta=metoda_rungekuttaIII(y1,t,n,h);
y=dokladne(t);
bibl=lsode("dot",y1,t);
gset xlabel"x";
gset ylabel"y";
gset grid;
gset title"1-rozwiazanie scisle";
plot(t,y,"r");
hold on;
pause
gset title"1-rozwiazanie scisle 2-rozwiazanie biblioteczne";
plot(t,bibl,"g");pause;hold on;
gset title"1-rozwiazanie scisle 2-rozwiazanie biblioteczne 3-rozwuiazanie Runge-Kutty 3/8";
plot(t,rungekutta,"b");hold on;
pause;

Gość

Huge cocks, wet pussies, hardcore action!
[link widoczny dla zalogowanych]