Симпсон әдісінің программа листингі

uses crt;

label 1,2,3;

var

a,b,eps,y,x,h,s,k,f1,t:real;

function f(x:real):real;

begin

end;

window(1,1,80,25);

textbackground(1);

clrscr;

writeln('* *');

writeln('* Данной программой вычисляется *');

writeln('* определенный интеграл *');

writeln('* методом Стмпсона *');

writeln('* *');

writeln('****************************************************');

writeln;

textbackground(0);

clrscr;

textcolor(13);

readln(a,b,n,eps);

h:=(b-a)/n;

s:=0;

j:=trunc(k);

1: f1:=f(a);

y:=f1;

y:=(y-f1)/2;

x1:=a+h;

for i:= 1 to j do

begin

f1:=f(x1);

f1:=f(x1+h);

y:=y+f1;

x1:=x1+2*h

end;

y:=y*2*h/3;

s:=y;

j:=j*2;

window(15,20,65,23);

textbackground(6);

clrscr;

writeln('Значение интеграла I=',y:7:4);

end.