Задача в приложенном документе.
Я не могу сообразить как мне сделать это:"Вычислить разность в количестве итераций"
Вот что у меня получилось:
Program Lab4_2;
Uses Wincrt;
Type
func=Function(xa:Real):Real;
Var
x1,x2,T1,S1,R1,R ,T2,S2, x:Real;
A,B,Exp,Eps:Real;
n1,n2,n:Word;
w:text;
Function Int(xa:Real):Real;Far;
Begin
Int:=((sin(xa*xa+0.5))/((cos(xa*xa+0.5))))/(1+2*xa*xa);
End;
Function MTrap(Aa,Ba:Real;na:Word;Int:func):Real;
Var i:Word; S,x,h:Real;
Begin
h:=(Ba-Aa)/na;
S:=(Int(Aa)+Int(Ba))*0.5;
x:=Aa+h;
For i:=1 to na-1 do
Begin S:=S+Int(x);x:=x+h;End;
MTrap:=S*h;
End;
Function MSims(Aa,Ba:Real;na:Word;Int:Func):Real;
Var i:Word;s,x,h:Real;c:shortint;
Begin
h:=(Ba-Aa)/na;
S:=Int(Aa)+Int(Ba);
c:=1; x:=Aa+h;
For i:=1 to na-1 do
Begin
S:=S+Int(x)*(c+3);
c:=-c;
x:=x+h;
End;
MSims:=S*h/3;
End;
Begin
WriteLn('1. Вычисление значений определенного интеграла');WriteLn;
WriteLn('Введите левую границу интервала');Read(A);
WriteLn('Введите правую границу интервала’);Read(B);
WriteLn('Введите число разбиений n');Read(n);
WriteLn('Введите точность Exp'); Read(Exp);
T1:=MTrap(A,B,n, Int);
repeat
n1:=2*n;
T2:=MTrap(A,B,n1,Int);
R:=Abs(T2-T1);
T1:=T2;
until R<Exp;
WriteLn('Интеграл решенный методом Трапеции равен ',T2:9:7);
WriteLn('Необходимое число интервалов разбиения: ',n1:6);
Writeln;
S1:=MSims(A,B,n, Int);
repeat
n2:=2*n;
S2:=MSims(A,B,n2,Int);
R:=Abs(S2-S1);
S1:=S2;
until R<Exp ;
WriteLn( 'Интеграл решенный методом Симпсона равен ',S2:9:7);
WriteLn('Необходимое число интервалов разбиения: ',n2:6);
End.