MATLAB Denemeleri

 
clc
clear
x(1)=1;
y(1)=2;
h=0.2;
for i=1:4;
x(i+1)=x(i)+h;

 a(i)=h*(1/(x(i)+y(i)))
 b(i)=h*(1/(x(i)+(h/2)+y(i)+(a(i)/2)))
 c(i)=h*(1/(x(i)+(h/2)+y(i)+(b(i)/2)))
 d(i)=h*(1/(x(i)+h+y(i)+c(i)))
 
 y(i+1)=y(i)+(1/6)*(a(i)+2*b(i)+2*c(i)+d(i))
end
    
% clc
% clear
% x(1)=0;
% y(1)=2;
% h=0.05;
% for i=1:20
    % x(i+1)=x(i)+h
    % a(i)=h*(1/(x(i)+y(i)))
    % b(i)=h*(1/(x(i)+(h/2)+y(i)+a(i)/2))
    % c(i)=h*(1/(x(i)+(h/2)+y(i)+b(i)/2))
    % d(i)=h*(1/(x(i)+h+y(i)+c(i)))
    % y(i+1)=y(i)+(a(i)+2*b(i)+2*c(i)+d(i))/6
% end
plot(x,y) % grafik çiz

Kodları online denemek için
http://octave-online.net/