program ex92;
{$apptype console}
const
n=80;                        
h=1/16;                      
function f(x:real):real;
begin
   result:=exp(-x*x/2)/1.414/pi;    
end;
var
  j:integer;
  k1,k2,k3,k4:real;
  x,y:array[0..n] of real;            
begin
y[0]:=0.5;                        
for j:=0 to n do                   
  x[j]:=0+j*h;
for j:=0 to n-1 do          
   begin
     k1:=f(x[j]);
     k2:=f(x[j]+0.5*h);
     k3:=f(x[j]+0.5*h);
     k4:=f(x[j]+h);
     y[j+1]:=y[j]+h*(k1+2*k2+2*k3+k4)/6;
   end;
for j:= 0 to n do
writeln('j=',j,' x=',x[j]:4:4,' y=',y[j]);
end.                   

老师:能否注解一下,谢谢

您看一下数值分析的书，上面有专门一章是讲这个的。书上有思思路和流程，我是按照它来写的。