Курсовая для первого курса, спасайте или выгонять меня... Опции

[Floyd]

Приветсвую всех. Прошу помоч сделать мне курсовую.
Тут отксерино все мое задание

http://img185.echo.cx/img185/2035/11by.th.gif

Я тут новенький и очень надеюсь на вашу помощ, еще неделя и нужно сдавать, если не сдам выгонят. Буду жутко признателен тому кто поможет. А отплатить... единственное что могу, так это сделать аватар любой сложности...
З.Ы.Хотя сомневаюсь что это кому то нужно сдесь...

Я присоединил изображение...

Сообщение отредактировано: volvo - 7.05.2005 21:25

[Floyd]

В общем, как я говорл, благодаря Вам я сдал курсач =).
А вот мой друг не сдал, почему то не пашит процедура решения интеграла, выдает ошибку Invalid floating operation 207
А в окне пишит
Runtime error 207 at 000F:0003
____
помогите разобратся в чем проблемка....
Интеграл там рещается двумя способами, способом парралелограмма и Симпсона.

program kursova_robota_PM_41_7_variant;
{Vykonav Schutskiy Pavlo}
uses graph,crt;
type vekt=array[1..1000] of real;
var n:integer; h,itg,g,t,int1,int2,s1,sh,intg,intgs:real; NN,NK,i:integer; w:text; ci,mi:vekt;
const a=0.2; b=0.3;
{_______________________________________}
      function ArcSin(x:real):real;
begin
 if Abs(x) = 1 then ArcSin:=0 else
 ArcSin:=ArcTan(x/(1-sqr(x))); end;
{---------------------------------------------------------------------------}

function f(x:real):real;
begin
   f:=ArcSin(sqrt(x))/sqrt(x*(1-x));
end;
{--------------------------------------------------------------------------}

procedure Grafika;
var gd, gm: integer;
    x,y: real;
begin
   gd:=detect;
   initgraph(gd, gm, '');
   SetBkColor(white);
   SetColor(4);
   Line(0, 100, 640, 100);
   Line(320, 0, 320, 480);
   SetColor(green);
   outtextxy(50,10,'Kursova Robota');  outtextxy(420,10,'Grafik Funkcii');
   outtextxy(50,20,'Vikonav Schutskiy Pavlo');
   outtextxy(50,30,'PM-41');
   PutPixel(320,240,0);
   outtextxy(325,245,'0');
      x:=0.2;
      While x<=0.9    do
begin
        y:=ArcSin(sqrt(x))/sqrt(x*(1-x));
        PutPixel(round(x*155)+350,round(y*155), 4);
        x:=x+0.00125;
end; repeat until keypressed; cleardevice;  closegraph;end;
{--------------------------------------------------------------------------}
procedure Zn(x,h:real);
var  c:array[1..8] of real;
i:integer; y:real;
begin
assign(w,'dan.pas'); rewrite(w);
writeln(w,'Znachennya funcii dlya kojnogo z X');
x:=0.2;  h:=0.0125;
for i:=1 to 8 do begin x:=x+h;
 y:=ArcSin(sqrt(x))/sqrt(x*(1-x));
c[i]:=y;
write(w,'f[',i,']=',c[i]:2 :3,'  ');
end;end;
{--------------------------------------------------------------------------}
procedure MaxMin(n:integer;h:real);
var x,xStart,xStop,xmin,xmax:extended;
begin
writeln(w,'');
writeln(w,'Znahodjennya MAX i MIN funcii:');
      xStart:=0.2;
      xStop:=0.3;
      x:=xStart;
      xmax:=x;   xmin:=x;
      while x<=xStop+h/2 do
begin
      if f(x)>f(xmax) then   xmax:=x;
      if f(x)<f(xmin) then  xmin:=x;
      x:=x+h;
end;
writeln(w,'Maksumalne znachennya funcii=',f(xmax):2 :3);
writeln(w,'Minimalne  znachennya funcii=',f(xmin):2 :3); writeln('');
end;
   {----------------------------------------------------------------------}
procedure Sum(x,h:real);
var Sum,y:real; i:integer; c:array [1..8] of real;
begin
for i:= 1 to 8  do begin
       if i>=6 then begin
    sum:=sum+x; end;
    x:=x+h;
end;
writeln(w,'Summa=',sum:3:3);writeln('');
end;
   {------------------------------------------------------------------------}
{nahogdenie integrala}
  procedure integral(x,y:vekt;var int,ints:real);
    var summa1,summa2,summa3,h:real;w,k,n,b:integer;
  begin
    summa1:=0; summa2:=0; summa3:=0; n:=8; h:=0.0125;
    w:=0;k:=0;n:=0;b:=0;
   for k:=2 to n do
    summa1:=summa1+2*y[k];
    int:=(h/2)*(y[1]+summa1+y[n+1]);
  for w:=2 to n-1 do         begin
    b:=w mod 2;
  if b=0 then
    summa2:=summa2+4*y[w] else
    summa3:=summa3+2*y[w]; end;
    ints:=(h/3)*(y[1]+summa2+summa3+4*y[n]+y[n+1]);
end;

{--------------------------------------------------------------------------}
Begin
    clrscr;
         Grafika;
assign(w,'dan.pas');rewrite(w);
      n:=8; h:=0.0125;
      s1:=0.2;
      Zn(s1,h);
      MaxMin(n,h);
      Sum(s1,h);
ci[1]:=0; n:=8; h:=0.0125;
  for i:=2 to (n+1) do begin
    ci[i]:=ci[i-1]+h;
  for i:=1 to (n+1) do
    mi[i]:=f(ci[i]);
      integral(ci,mi,intg,intgs);  end;
writeln(w,'Pershiy integral= ',intg:3 :3);
writeln(w,'Drugiy integral=  ',intgs:3 :3); reset(w);  end.