лаба4 (Готовые лабы в ПАСКАЛЕ (описание внутри))
Описание файла
Файл "лаба4" внутри архива находится в следующих папках: c, Lab4 (4ая). Документ из архива "Готовые лабы в ПАСКАЛЕ (описание внутри)", который расположен в категории "". Всё это находится в предмете "численные методы" из 3 семестр, которые можно найти в файловом архиве МГТУ им. Н.Э.Баумана. Не смотря на прямую связь этого архива с МГТУ им. Н.Э.Баумана, его также можно найти и в других разделах. Архив можно найти в разделе "лабораторные работы", в предмете "численные методы" в общих файлах.
Онлайн просмотр документа "лаба4"
Текст из документа "лаба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