Скачиваний:
27
Добавлен:
09.12.2013
Размер:
896 б
Скачать
Program Simpson;
Uses CRT;
Const a=2;
b=10;
e=0.000001;

Var a_,b_,r,rest,
int1,int2,int:real;
i:integer;

Function f(x:real):real;
begin
f:=x*ln(1/x)/(1-x*x);
end;

Function Sim(a,b:real):real;
begin
Sim:=((b-a)/6)*(f(a)+4*f((a+b)/2)+f(b));
end;

Begin
clrscr;
a_:=a;b_:=b;
int:=0;rest:=0;i:=1;
repeat
int1:=Sim(a_,b_);
int2:=Sim(a_,(a_+b_)/2)+Sim((a_+b_)/2,b_);
r:=abs(int2-int1)/15;
if r<=e*(b_-a_)/(b-a) then
begin
int:=int+int2;
rest:=rest+r;
writeln('a[',i,']=',a_:0:6);
writeln('b[',i,']=',b_:0:6);
writeln('I[',i,']=',int2:0:15);
writeln('R[',i,']=',r:0:15);
writeln('E[',i,']=',e*(b_-a_)/(b-a):0:15);
writeln;
readkey;
a_:=b_;
b_:=b;
i:=i+1;
end
else b_:=(b_+a_)/2;
until b=a_;
writeln('SumI=',int:0:15);
writeln('SumR=',rest:0:15);
readln;
end.
Соседние файлы в папке МЕТОДС~1