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

dufi

jak ktos posiada dzialajacy program runge- kutty-formula3/8 to niech wrzuci z łaski swej Smile

z gory dzieki

vh

prawie to czego potrzebujesz

function formula38()
y0=2
x=linspace(-5,5,11);
n=length(x)
h=abs(abs(x(2))-abs(x(1)))
rozw=formula(y0,x,n,h);
plot(x,rozw)
pause

function df=fun(x,y)
df=0.5*x*y;

function df=dot(y,x)
df=0.5*x*y;

function f=dokl(x)
f=exp(x^2/4);

function rungekutta=formula(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)*k1);

k3=h*fun(x(i-1)+(2/3)*h,rungekutta(i-1)-(1/3)*k1+k2);

k4=h*fun(x(i-1)+h,rungekutta(i-1)+k1+k2);

rungekutta(i)=rungekutta(i-1)+(1/ Cool

*(k1+3*k2+3*k3+k4);
end