simpson.pas

来自「关于课件 好多种结构的算法 等等 大家自己看吧」· PAS 代码 · 共 30 行

program kk(input,output);
var
s1,s2:real;
function f1(x:real):real;
begin
 f1:=1/(x+1);
end;
function f2(x:real):real;
begin
 f2:=x*exp(-x*x)
end;
function simp(FUNCTION f(x:real):real;a,b:real):real;
const
 n=40;
 var
 h,s:real;
 i:integer;
begin
 h:=(b-a)/n;
 s:=f(a)+f(b);
 for i:=1 to n-1 do
  if odd(i) then s:=s+4*f(a+i*h)
  else s:=s+2*f(a+i*h);
 simp:=h/3*s
end;
begin
s1:=simp(f1,0,1);
s2:=simp(f2,1,2.5);
writeln(s1,s2);
end;