solving systems of ODEs...(ikinci dereceden türevi olan sistemleri çözmek)

%ikinci dereceden türevi verilen denklemleri çözmek için o ikinci dereceyi

%iki tane birinci dereceye parçalayarak çözeriz y dediğim birinci
%dereceden türev,z ikinci dereceden türev değri
%f dosyasını açmayı unutmayın
disp('baslangıc degerlerni griniz:');

y(1)=input('y(1):');

z(1)=input('z(1)');

h=input('h:');

a=input('üst sınır degri:');

b=input('alt sınır degeri');

n=(b-a)/h;

x=a:h:b;

for i=1:n

y(i+1)=y(i)+z(i)*h;

z(i+1)=z(i)+f(x(i),y(i))*h;

end

fprintf('ikinci türev degeri=%f',z(n+1));

f FONKSİYONU

function y=f(x,t);

y=x-t;

OUTPUT
baslangıc degerlerni griniz:


y(1):2

z(1)0

h:0.1

üst sınır degri:0

alt sınır degeri5

ikinci türev degeri=3.127871>>