Пример: Транспортная логистика
Я ищу:
На главную  |  Добавить в избранное  

Кибернетика /

Вычисление интеграла методом Ньютона-Котеса (теория и программа на Паскале)

←предыдущая следующая→
1 2 3 4 5 6 7 8 9 



Скачать реферат


================================================

{$N+}

unit k_graph;

interface

uses

crt,graph,k_unit,graphic;

procedure winwin1;

procedure proline(ea:word);

procedure winwwodab(ea:word);

procedure error1(ea:word);

procedure helpwin(ea:word);

procedure error(ea:word);

procedure newsctext(ea:word);

procedure newsc(ea:word);

procedure win1(ea:word);

procedure win2(ea:word;var k:word);

procedure wwodn(ea:word;var n:integer);

procedure wwodab(ea:word;var a,b:real);

procedure wwod1(ea:word;var y:array of double;var n:integer;var a,b:real);

procedure wwod2(ea:word;var ea1:word;var n:integer;var a,b:real;var st:string);

procedure win3(ea:word;n:integer;a,b:real;int:double;f:string;h:array of double;var k:word);

implementation

procedure proline(ea:word);

{Проседура полосы процесса}

var

i:integer;

f:string;

c:char;

begin

newsc(ea);

setcolor(15);

setfillstyle(1,7);

bar(160,150,460,260);

rectangle(165,155,455,255);

rectangle(167,157,453,253);

case (ea mod 2) of

0: outtextxy(180,170,' Идет работа .Ждите..');

1: outtextxy(180,170,' Working.Please wait..');

end;

setfillstyle(1,12);

setcolor(0);

rectangle(200,199,401,221);

for i:=1 to 9 do

line(200+i*20,200,200+i*20,220);

delay(20000);

for i:=1 to 100 do

begin

if ((i-1) mod 10)=0 then

line(200+((i-1) div 10)*20,200,200+((i-1) div 10)*20,220);

bar(round(200+2*(i-0.5)),200,200+2*i,220);

delay(1100);

setcolor(15);

setfillstyle(1,7);

bar(280,230,323,250);

str(i,f);

f:=f+'%';

outtextxy(290,235,f);

if (i mod 25) =0 then

bar(170,180,452,198);

if (ea mod 2)=0 then

case (i div 25) of

0:

outtextxy(170,190,'Подготовка ');

1:

outtextxy(170,190,'Расчет коеффициентов в многочлене');

2:

outtextxy(170,190,'Расчет коеффициентов Ньютона-Котеса');

3:

outtextxy(170,190,'Расчет интеграла');

end

else

case (i div 25) of

0:

outtextxy(170,190,'Prepearing');

1:

outtextxy(170,190,'Calculation of mnogochlen coeff.');

2:

outtextxy(170,190,'Calculation of Newton-Cotes coeff. ');

3:

outtextxy(170,190,'Calculation of integral');

end;

setfillstyle(1,12);

setcolor(0);

end;

end;

procedure winwwodn(ea:word);

{Окно ввода числа узлов интерполяции}

var

c:char;

f:string;

begin

helpwin(ea);

if (ea mod 2) =0 then

begin

outtextxy(360,140,' В этом окне необходимо ');

outtextxy(360,155,' ввести количество узлов ');

outtextxy(360,170,' интерполяции, от которого ');

outtextxy(360,185,' будет зависить точность ');

outtextxy(360,200,' вычисления интеграл и ');

outtextxy(360,215,' количество зн чений функции.');

outtextxy(360,240,' ВНИМАНИЕ : НАСТОЯТЕЛЬНО ');

outtextxy(360,250,' РЕКОМЕНДУЕТСЯ НЕ ВВОДИТЬ ');

outtextxy(360,260,' ЗНАЧЕНИЕ N БОЛЬШЕ 12 !! ');

end

else

begin

outtextxy(360,140,' In this window you have to ');

outtextxy(360,155,' put into the number. ');

outtextxy(360,170,' The accuracy of calculation ');

outtextxy(360,185,' and the number of function ');

outtextxy(360,200,' parameters will depend on ');

outtextxy(360,215,' this number. ');

outtextxy(360,240,' WARNING: IT IS HARDLY ');

outtextxy(360,250,' RECOMENDED NOT TO PUT IN ');

outtextxy(360,260,' NUMBER MORE THEN 12 !! ');

end;

setcolor(2);

setfillstyle(1,14);

bar(70,200,340,300);

rectangle(75,205,335,295);

rectangle(77,207,333,293);

if (ea mod 2) =0 then

begin

outtextxy(90,227,'Введите количество узлов(n):');

outtextxy(80,270,'ВНИМАНИЕ: При больших n возможна');

outtextxy(80,280,'некорректная работа компьютера!!');

end

else

begin

outtextxy(80,217,'Put in number of');

outtextxy(80,227,' interpolation units:');

outtextxy(80,270,'WARNING:if you use big number ');

outtextxy(80,280,'of units,PC wont work properly!');

end;

setfillstyle(1,0);

bar(190,240,230,255);

end;

procedure wwodn(ea:word;var n:integer);

{Процедура ввода узлов n}

var

ec,p:integer;

k,f:string;

x:integer;

c:char;

begin

newsc(ea);

winwwodn(ea);

repeat

repeat

winwwodn(ea);

gotoxy(25,16);

read(k);

val(k,p,ec);

if ec0 then

begin

error1(ea);

readln;

end;

until ec=0;

n:=p;

if n>12 then

begin

if keypressed then

c:=readkey;

c:='r';

setcolor(15);

setfillstyle(1,12);

bar(140,210,490,300);

rectangle(145,215,485,295);

rectangle(147,217,483,293);

if (ea mod 2) =0 then

begin

outtextxy(150,227,' Предупреждение!');

outtextxy(150,237,' Вы дейcтвительно хотите использовать');

outtextxy(150,250,' большое значение N ???');

end

else

begin

outtextxy(150,227,' Warning!! ');

outtextxy(150,237,' Do you realy want to use a big ');

outtextxy(150,250,' number interpolation units(N)??? ');

end;

sound(600);

delay(4000);

nosound;

setfillstyle(1,2);

bar(320,260,350,280);

setfillstyle(1,12);

bar(250,260,280,280);

repeat

if keypressed then

begin

c:=readkey;

if (c=#80) or (c=#72) or (c=#77) or (c=#75) then

x:=x+1;

setfillstyle(1,2);

if (x mod 2)=0 then

begin

bar(250,260,280,280);

setfillstyle(1,12);

bar(320,260,350,280);

end

else

begin

bar(320,260,350,280);

setfillstyle(1,12);

bar(250,260,280,280);

END;

end;

if (ea mod 2) =0 then

begin

outtextxy(255,267,'ДА');

outtextxy(325,267,'НЕТ');

end

else

begin

outtextxy(255,267,'YES');

outtextxy(325,267,'NO');

end;

until c=#13;

if abs(x mod 2)=1 then

begin

n:=0;

setcolor(15);

setfillstyle(1,2);

bar(160,200,460,280);

rectangle(165,205,455,275);

rectangle(167,207

←предыдущая следующая→
1 2 3 4 5 6 7 8 9 



Copyright © 2005—2007 «Mark5»