Готовые лабы в ПАСКАЛЕ (описание внутри) (лаба4)

Посмотреть архив целиком







program integral;

var

i,j,n:integer;

intr,int1,int2,ints,int3,int4,h,a,b,eps,xi,r:real;

BEGIN

writeln('':4000);

writeln('Tochnost eps=0.0001');

eps:=0.0001;

a:=0;

b:=4;

{metod 1}

n:=trunc((b-a)/sqrt(eps)+1);

h:=(b-a)/n;

for i:=1 to n do

begin

xi:=a+(i-1/2)*h;

int1:=int1+h*4*xi*xi*xi/(xi*xi*xi*xi*xi*xi*xi*xi+1);

end;

writeln(' ');

writeln('Metod pryamougolnikov');

writeln(' ');

writeln('Integral bez utochnenia ',int1:20:20);

h:=h/2;

xi:=0;

n:=2*n;

for i:=1 to n do

begin

xi:=a+(i-1/2)*h;

int2:=int2+h*4*xi*xi*xi/(xi*xi*xi*xi*xi*xi*xi*xi+1);

end;

r:=abs(int2-int1)/3;

intr:=int1+(int2-int1)/3;

writeln('Integral s utochneniem ',intr:20:20);

writeln('Ocenka pogreshnosti R(h/2)',r:20:20);

writeln('Chislo razbieniy ',n);

{metod 2}

writeln(' ');

writeln('Metod Simpsona');

writeln(' ');

n:=trunc((b-a)/(2*sqrt(sqrt(eps)))+1);

h:=(b-a)/(2*n);

xi:=0;

for i:=1 to 2*n-1 do

begin

xi:=a+i*h;

If odd(i) then int3:=int3+4*h/3*4*xi*xi*xi/(xi*xi*xi*xi*xi*xi*xi*xi+1);

If not odd(i) then int3:=int3+2*h/3*4*xi*xi*xi/(xi*xi*xi*xi*xi*xi*xi*xi+1);

end;

writeln('Integral bez utochnenia ',int3:20:19);

xi:=0;

n:=n*2;

h:=h/2;

for i:=1 to 2*n-1 do

begin

xi:=a+i*h;

If odd(i) then int4:=int4+4*h/3*4*xi*xi*xi/(xi*xi*xi*xi*xi*xi*xi*xi+1);

If not odd(i) then int4:=int4+2*h/3*4*xi*xi*xi/(xi*xi*xi*xi*xi*xi*xi*xi+1);

end;

r:=abs(int4-int3)/15;

ints:=int4+(int4-int3)/15;

writeln('integral s utochneniem ',ints:20:19);

writeln('Ocenka pogreshnosti R(h/2)',r:20:19);

writeln('Chislo razbieniy ',n);

readln;

END.








Результат работы программы.







Tochnost eps=0.0001


Metod pryamougolnikov


Integral bez utochnenia 1.56689011660

Integral s utochneniem 1.56689011160

Ocenka pogreshnosti R(h/2) 0.00000000498

Chislo razbieniy 802


Metod Simpsona


Integral bez utochnenia 1.56675104780

integral s utochneniem 1.56683249840

Ocenka pogreshnosti R(h/2) 0.00000509066

Chislo razbieniy 42





Случайные файлы

Файл
CBRR4279.DOC
9972-1.rtf
47690.rtf
44702.doc
10757-1.rtf




Чтобы не видеть здесь видео-рекламу достаточно стать зарегистрированным пользователем.
Чтобы не видеть никакую рекламу на сайте, нужно стать VIP-пользователем.
Это можно сделать совершенно бесплатно. Читайте подробности тут.