مهندس این جوری باید بشه؟
% Runge8clear('all'), close('all'), clc
h=0.5;
x0=0;
xf = 2;
x = x0 : h : xf; x = x(:);
y0 = [4,6];
y = zeros(length(x), length(y0));
y(1,:) = y0(:).';
for i=1:length(x)-1
k1=RK(x(i),y(i,:));
k11 = k1(1);
k12 = k1(2);
ytmp1 = y(i,1) + k11*h/2...