RK4 help (someone help fix this)
이전 댓글 표시
so i made the code but i only get the first two values then the values start to equal zeros wich shouldnt hapend.
this is the code :
x_start=0; %valeeur initial de x
x_end=1; %valeur final de x
n_nodes=20;
h=(x_end-x_start)/(n_nodes-1);
x=linspace(x_start, x_end, n_nodes);
y_RK4=zeros(1,n_nodes);
y_RK4(1)=3;
f=@(x,y) (-0.3*y+2*exp(-x));
for i=1,n_nodes-1;
k1=h*f(x(i),y_RK4(i));
k2=h*f((x(i)+(h/2)),(y_RK4(i)+(k1/2)));
k3=h*f((x(i)+(h/2)),(y_RK4(i)+(k2/2)));
k4=h*f((x(i)+h),(y_RK4(i)+k3));
y_RK4(i+1)=y_RK4(i)+(k1+(2*k2)+(2*k3)+k4)/6;
end
채택된 답변
추가 답변 (0개)
카테고리
도움말 센터 및 File Exchange에서 Database Toolbox에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!