3. Алгоритм:
3.1. Подпрограммы:
{Функция, площадь которой нужно вычислить}
m,x f
{Вычисления по формуле многоугольников}
a,b,dx,m,toch y = 0 x = (a + dx)/2 x <=
b yes y = y + f(m,x) x = x + dx no y
= y * dx prm
=(abs((y-toch)/toch))*100 prm
x <
b a,b,dx,m,toch y = 0 x = (a + dx)/2 yes y = y + f(m,x) x = x + dx no y=
(y + (f(m,a)+f(m,b))/2) * dx
trp=(abs((y-toch)/toch))*100 trp
{Вычисление
по формуле трапеций}
{Вычисление по методу Симпсона}
a,b,dx,m,toch,n y1
= 0 x = (a + dx)/2 yes no
x <
b-dx
y1
= y1 + f(m,x) x = x + 2dx y2 = 0 x = a + 2*dx
x <
b-2dx yes no
y2
= y2 + f(m,x) x = x + 2dx A
A y
= ((b-a)/(3*n)) * (f(m,a)+f(m,b) + 4*y1 + 2*y2)
sim=(abs((y-toch)/toch))*100 sim
m a=0 b=pi/2 toch=pi/4 n=200 dx=
(b-a)/n tab[1,0]:=200; tab[1,1]:=prm(a,b,dx,m,toch); tab[1,2]:=trp(a,b,dx,m,toch); tab[1,3]:=sim(a,b,dx,m,toch,n); n:=2*n;
dx:=dx/2; tab[1,4]:=prm(a,b,dx,m,toch); tab[1,5]:=trp(a,b,dx,m,toch); tab[1,6]:=sim(a,b,dx,m,toch,n); n=400 dx=
(b-a)/n tab[1,0]:=400; tab[1,1]:=prm(a,b,dx,m,toch); tab[1,2]:=trp(a,b,dx,m,toch); tab[1,3]:=sim(a,b,dx,m,toch,n); n:=2*n;
dx:=dx/2; tab[1,4]:=prm(a,b,dx,m,toch); tab[1,5]:=trp(a,b,dx,m,toch); tab[1,6]:=sim(a,b,dx,m,toch,n);
3.2.Главная
программа:
B
B n=3000 dx=
(b-a)/n tab[1,0]:=3000; tab[1,1]:=prm(a,b,dx,m,toch); tab[1,2]:=trp(a,b,dx,m,toch); tab[1,3]:=sim(a,b,dx,m,toch,n); n:=2*n;
dx:=dx/2; tab[1,4]:=prm(a,b,dx,m,toch); tab[1,5]:=trp(a,b,dx,m,toch); tab[1,6]:=sim(a,b,dx,m,toch,n); i=1 j=0 i<4 yes j<6 yes tab[i,j] j=j+1 i=i+1
