Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Скачиваний:
14
Добавлен:
11.04.2015
Размер:
1.01 Кб
Скачать
program kursovaya;
var r,r2,v,int,a,b:real;
n:integer;
function f(x:real):real;
begin
f:=cos(x)*exp(-x*x);
end;
procedure simpson;
var
l,h:real;
m:integer;
mn:integer;
s:real;
begin
h=:(b-a)/n;
s:=f(a)+f(b);
mn:=4;
for m:=1 to n-1 do
begin
s:=s+mn*f(a+h*m);
if(mn=4) then mn=:2
else mn:=4;
end;
int:=s*h/3;
l:=25;
r:=((b-a)*h*h*h*h/120)*l;
end;
procedure trap;
var
l,h,s:real;
m:=integer;
begin
h:=(b-a)/n;
v:=0;
s:=(f(a)+f(b))/2;
m:=1;
for m:=1 to n-1 do
begin
s:=s+f(a+m*h);
end;
v:=s*h;
l:=25;
r2:=((b-a)*h*h*h*h/120)*1;
end;
begin
write('vvedite ni)i(niy predel integrirovani9=');
readln(a);
write('vvedite bepxniy predel integrirovani9=');
readln(b);
write('vvedite koli4ectvo IIIagov integrirovani9=');
readln(n);
simpson;
writeln('zna4enie integrala metodom simpsona na dannom otre3ke=',int:12:10);
writeln('pogreIInoctb=',r2:12:10);
trap;
writeln('3na4enie integrala metodom trapecii=',v:12:10);
writeln('pogreIIInoctb=',r2:12:10);
readln;
end.
Соседние файлы в папке sam