به خاطر اینه که از متغیرهای نمادین استفاده کردید. این اصلاح شده اش:
clear, clc
disp('RK4 ')
h=0.3;
xz=0;
yz=1;
f=inline('y+x^2');
x(1)=xz;
y(1)=yz;
for i=2:99
k1(i)=h*f(x(i-1),y(i-1));
k2(i)=h*f(x(i-1)+(h/2),y(i-1)+(k1(i)/2));
k3(i)=h*f(x(i-1)+(h/2),y(i-1)+(k2(i)/2))...