سلام
نماز و روزه شما هم قبول باشه
به یک همچین فرمی باید بنویسید:
clc
h = 0.5;
x0 = 0;
y = [4; 6];
for i=1:4;
x=x0+i*h;
k1=RK(t,y);
% k2=...;
% k3=...;
% k4=...;
y=y+(h/6)*(k1+2*k2+2*k3+k4)
end
و ورودی تابع دینامیک هم باید درست بشه:
function eq=RK(x, y);
eq =...