Помощь - Поиск - Пользователи - Календарь
Полная версия: Метод Трапеции и Симпсона
Форум «Всё о Паскале» > Pascal, Object Pascal > Задачи
Selena
Помогите пожалуйста дописать программу.
Задача в приложенном документе.

Я не могу сообразить как мне сделать это:"Вычислить разность в количестве итераций"

Вот что у меня получилось:
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.
Гость
правилно
Это текстовая версия — только основной контент. Для просмотра полной версии этой страницы, пожалуйста, нажмите сюда.