اون کار دیگه با توجه به تعریف جمع ماتریسها و ضرب اسکالر در یک ماتریس کار ساده ای هستش. مثلاً برای k1 و ytmp می تونید بنویسید:
k1=RK(x(i),y(i,:));
k11 = k1(1);
k12 = k1(2);
ytmp1 = y(i,1) + k11*h/2;
ytmp2 = y(i,2) + k12*h/2;
بقیه هم مشابه به همینه.
می تونید تابع RK رو هم تغییر...