Дослідження методів чисельного інтегрування
Курсовой проект - Компьютеры, программирование
Другие курсовые по предмету Компьютеры, программирование
p>
if (i mod 25 =0) then
begin
readln;
bar(360,40,620,450);
end;
outtextxy(430,(i mod 25)*15+40,s0+ +s1);
end;
end;
{-----------------------------------------------------------------------}
begin
InitGraphMode;
setcolor(blue);
h:=0.1;
repeat
cleardevice;
Item[1]:=Метод Чебишева 3-го порядку;
Item[2]:=Метод Чебишева 4-го порядку;
Item[3]:=Метод Чебишева 5-го порядку;
Item[4]:=Крок обчислень=0.1;
Item[5]:=Крок обчислень=0.2;
Item[6]:=Крок обчислень=0.5;
Item[7]:=Про автора;
Item[8]:=Вихид;
i:=1;
setcolor(white);
Menu(50,30,Item,7,i);
case i of
1:begin
setbkcolor(lightblue);x:=a; sum:=0; i:=0;
while x<=b do begin
inc(i);
x:=x+h;
si:=0;
for k:=1 to 3 do si:=si+f(x-h/2+h/2*t3[k]);
si:=h/3*si;
sum:=sum+si;
sum1:=(x-((a+b)/2))*(x-((a+b)/2))*(x-((a+b)/2))*(x-((a+b)/2))*
48*exp(-2*x*x+1.5)-384*x*x*exp(-2*x*x+1.5)+256*x*x*x*x*exp(-2*x*x+1.5);
n:=i+10;
end;
delta:=((b-a)/(n*fact(3)))*sum1;
delta:=Abs(delta);
AllDelta:=abs(delta);
str(Delta:12:10,s1);
Resm[i]:=Delta;
setcolor(red);
str(Sum:12:18,s);
str(AllDelta:12:20,s1);
str(h:2:3,hs);
setcolor(blue);
Outtextxy(100,240,Интеграл: +s);
Outtextxy(45,260,Точне значення (Mathcad): 2.68069545714219);
Outtextxy(100,280,Похибка: +s1);
Outtextxy(100,300,Крок: +hs);
readln;
end;
2: begin
setbkcolor(lightblue);x:=a; sum:=0; i:=0;
while x<=b do begin
inc(i);
x:=x+h;
si:=0;
for k:=1 to 4 do si:=si+f(x-h/2+h/2*t4[k]);
si:=h/4*si;
sum:=sum+si;
sum1:=(x-((a+b)/2))*(x-((a+b)/2))*(x-((a+b)/2))*(x-((a+b)/2))*(x-((a+b)/2))*
(-960*x*exp(-2*x*x+1.5)+2560*x*x*x*exp(-2*x*x+1.5)-1024*x*x*x*x*x*exp(-2*x*x+1.5));
n:=i+10;
end;
delta:=((b-a)/(n*fact(3)))*sum1;
delta:=Abs(delta);
AllDelta:=abs(delta);
str(Delta:12:10,s1);
Resm[i]:=Delta;
setcolor(red);
str(Sum:12:18,s);
str(AllDelta:12:20,s1);
str(h:2:3,hs);
setcolor(blue);
Outtextxy(100,240,Интеграл: +s);
Outtextxy(45,260,Точне значення (Mathcad): 2.68069545714219);
Outtextxy(100,280,Похибка: +s1);
Outtextxy(100,300,Крок: +hs);
readln;
end;
3:begin
setbkcolor(lightblue);x:=a; sum:=0; i:=0;
while x<=b do begin
inc(i);
x:=x+h;
si:=0;
for k:=1 to 5 do si:=si+f(x-h/2+h/2*t5[k]);
si:=h/5*si;
sum:=sum+si;
sum1:=(x-((a+b)/2))*(x-((a+b)/2))*(x-((a+b)/2))*(x-((a+b)/2))*(x-((a+b)/2))*
(x-((a+b)/2))*(-960*exp(-2*x*x+1.5)+11520*x*x*exp(-2*x*x+1.5)
-15360*x*x*x*x*exp(-2*x*x+1.5)+4096*x*x*x*x*x*x*exp(-2*x*x+1.5));
n:=i+10;
end;
delta:=((b-a)/(n*fact(3)))*sum1;
delta:=Abs(delta);
AllDelta:=abs(delta);
str(Delta:12:10,s1);
Resm[i]:=Delta;
setcolor(red);
str(Sum:12:18,s);
str(AllDelta:12:20,s1);
str(h:2:3,hs);
setcolor(blue);
Outtextxy(100,240,Интеграл: +s);
Outtextxy(45,260,Точне значення (Mathcad): 2.68069545714219);
Outtextxy(100,280,Похибка: +s1);
Outtextxy(100,300,Крок: +hs);
readln;
end;
4: begin
sound(500);
setbkcolor(lightred);
h:=0.01;
Outtextxy(150,200,Крок h=0.01 );
Outtextxy(170,250,Натиснить будь-яку клавишу ...);
delay(10000);
nosound;
repeat until keypressed;
end;
5: begin
setbkcolor(lightred);
h:=0.02;
Outtextxy(150,200,Крок h=0.02 );
Outtextxy(170,250,Натиснить будь-яку клавишу ...);
repeat until keypressed;
end;
6: begin
setbkcolor(lightred);
h:=0.05;
Outtextxy(150,200,Крок h=0.05 );
Outtextxy(170,250,Натиснить будь-яку клавишу ...);
repeat until keypressed;
end;
7: begin
setbkcolor(lightred);
Outtextxy(160,200,Програму виконав ст.гр. 3АВ-02);
Outtextxy(210,220,Нетяга О.В.);
repeat until keypressed;
end;
8: exit;
end;
until (i=7) or (i=-1);
CloseGraph;
end.
Додаток Б