[Hiter]

program Runge_Kutta;
uses crt;
var dy,h,yn,sn,k1,k2,y,k3,k4,x0,y0:real;
t:integer;
out:text;
procedure f(x,y:real);
begin
dy:=3*x+y/2;
end;
BEGIN
clrscr;
assign(out,'sonuc.txt');rewrite(out);
x0:=0;
y0:=1;
h:=0.1;
t:=0;
sn:=20;
writeln(out,'k1':5,'k2':10,'k3':10,'k4':10,'Yn':10 );
writeln(out,'************************************* *********':48);
repeat
t:=t+1;

f(x0,y0);
k1:=h*dy;

f((x0+h/2),(y0+k1/2));
k2:=h*dy;

f((x0+h/2),(y0+k2/2));
k3:=h*dy;

f((x0+h),(y0+k3));
k4:=h*dy;

yn:= (k1+2*k2+2*k3+k4)/6;
y:=y0+yn;

writeln(out,'y',t,'=',k1:2:5,k2:10:5,k3:10:5,k4:10 :5,y:10:5);
y0:=y;
until t>=sn;
close(out);
END.

Buda sonuç txt'in içeriği:

k1 k2 k3 k4 Yn
**********************************************
y1=0.05000 0.06625 0.06666 0.08333 1.06652
y2=0.05333 0.06966 0.07007 0.08683 1.13646
y3=0.05682 0.07324 0.07365 0.09051 1.20998
y4=0.06050 0.07701 0.07742 0.09437 1.28727
y5=0.06436 0.08097 0.08139 0.09843 1.36852
y6=0.06843 0.08514 0.08555 0.10270 1.45394
y7=0.07270 0.08951 0.08993 0.10719 1.54374
y8=0.07719 0.09412 0.09454 0.11191 1.63814
y9=0.08191 0.09895 0.09938 0.11688 1.73738
y10=0.08687 0.10404 0.10447 0.12209 1.84172
y11=0.09209 0.10939 0.10982 0.12758 1.95140
y12=0.09757 0.11501 0.11544 0.13334 2.06670
y13=0.10333 0.12092 0.12136 0.13940 2.18791
y14=0.10940 0.12713 0.12757 0.14577 2.31534
y15=0.11577 0.13366 0.13411 0.15247 2.44931
y16=0.12247 0.14053 0.14098 0.15951 2.59014
y17=0.12951 0.14774 0.14820 0.16692 2.73819
y18=0.13691 0.15533 0.15579 0.17470 2.89383
y19=0.14469 0.16331 0.16377 0.18288 3.05746
y20=0.15287 0.17169 0.17217 0.19148 3.22947



ALINTI DEĞİLDİR KENDİ ÇALIŞMAM...