ru.algorithms

 
 - RU.ALGORITHMS ----------------------------------------------------------------
 From : Yuri Burger                          2:468/85.3     29 May 2001  19:41:08
 To : All
 Subject : Integral: вычисление определенных интегралов.....
 -------------------------------------------------------------------------------- 
 

 
 27 May 01 13:55, Michael Komm wrote to All:
 
  MK> Теперь другой вопрос - нет ли у кого сырца/алгоритма, как определенный
 
     Есть сырец, писан как минимум тремя поколениями студентов :)
     Кстати, есть еще куча лаб по "численным методам", оп всяким симплексам и
 прочее, прочее.... Как освобожусь от диплома, попробую разгрести и с разрешения
 начальства пораздовать (а то вроди народ часто спрашует)...
 
     Сорри, что сообшэния в проге на украинском... не мной писано :) есть и дока
 с описаловом, теорией и т.д. но она тоже на украинском.
 [ю]ДДДДДДДД Begin MV10_MI.PAS ДДДДДДД
 Uses Crt;
 var a,b,n,y: real;
 
 function int_1(x: real): real;
   begin
     int_1:=sqrt(0.8*x*x+1)/(x+sqrt(1.5*x*x+2))
   end;
 
 function int_2(x: real): real;
   begin
     int_2:=sin(x+1.4)/(0.8+cos(2*x*x+0.5))
   end;
 
 function proiz1(x: real): real;
   begin
     proiz1:=(1.6*x*(x+sqrt(1.5*x*x+2))/(2*sqrt(0.8*x*x+1))-
     sqrt(0.8*x*x+1)*(1+3*x/(2*sqrt(1.5*x*x+2))))/sqr(x+sqrt(1.5*x*x+2));
   end;
 
 function proiz2(x: real): real;
   begin
     proiz2:=(cos(x+1.4)*(0.8+cos(2*x*x+0.5))+
     sin(x+1.4)*sin(2*x*x+0.5)*4*x)/sqr(0.8+cos(2*x*x+0.5));
   end;
 
 procedure left_rect(a,b,n: real; var y: real);
  var h,sr: real;
   begin
     h:=(b-a)/n; sr:=a+(b-a)/2; y:=0;
     while a<=b do
       begin
 y:=y+int_1(a);
 a:=a+h
       end;
     y:=h*y+sqr(b-a)/(2*n)*proiz1(sr)
   end;
 
 procedure right_rect(a,b,n: real; var y: real);
  var h,sr: real;
   begin
     h:=(b-a)/n; sr:=a+(b-a)/2; y:=0;
     while a<=b do
       begin
         y:=y+int_1(a);
         a:=a+h
       end;
     y:=h*y-sqr(b-a)/(2*n)*proiz1(sr)
   end;
 
 procedure sr_rect(a,b,n: real; var y: real);
  var h,sr: real;
   begin
     h:=(b-a)/n; sr:=a+(b-a)/2; y:=0;
     while a<=b do
       begin
         y:=y+int_2(a+h/2);
         a:=a+h
       end;
     y:=h*y+sqr(b-a)*(b-a)/(24*n*n)*(1+proiz2(sr));
   end;
 
 begin
   textattr:=blue*16+yellow;
   clrscr;
   gotoxy(10,10);
   write('Розрахунок iнтегpалу по фоpм. лiвих пpямокутникiв');
   gotoxy(10,12); write('Уведiть a='); read(a);
   gotoxy(10,13); write('Уведiть b='); read(b);
   gotoxy(10,14); write('Уведiть n='); read(n);
   left_rect(a,b,n,y);
   gotoxy(10,16); write('Значення iнтегpалу : ',y);
   readln; readln;
   clrscr;
   gotoxy(10,10);
   write('Розрахунок iнтегpалу по фоpм. пpавих пpямокутникiв');
   gotoxy(10,12); write('Уведiть a='); read(a);
   gotoxy(10,13); write('Уведiть b='); read(b);
   gotoxy(10,14); write('Уведiть n='); read(n);
   right_rect(a,b,n,y);
   gotoxy(10,16); write('Значення iнтегpалу : ',y);
   readln; readln;
   clrscr;
   gotoxy(10,10);
   write('Розpахунок iнтегpалу по фоpм. сеpеднiх пpямокутникiв');
   gotoxy(10,12); write('Уведiть a='); read(a);
   gotoxy(10,13); write('Уведiть b='); read(b);
   gotoxy(10,14); write('Уведiть n='); read(n);
   sr_rect(a,b,n,y);
   gotoxy(10,16); write('Значення iнтегpалу : ',y);
   readln; readln;
 end.
 [ю]ДДДДДДДД End MV10_MI.PAS   ДДДДДДД
                                                  Kрюгер.
 ---
  * Origin: А хто тут есть, у кого есть за что поесть? (2:468/85.3)
 
 

Вернуться к списку тем, сортированных по: возрастание даты  уменьшение даты  тема  автор 

Тема:    Автор:    Дата:    Integral   Michael Komm   27 May 2001 13:55:40   Integral   Michail Svarichevsky   28 May 2001 16:06:01   Integral   Eugene Eremin   29 May 2001 08:13:41   Integral   Sergey Naidenov   29 May 2001 17:09:19   Integral: вычисление определенных интегралов.....   Yuri Burger   29 May 2001 19:41:08   Integral   Dmitriy Nesmachny   01 Jun 2001 09:10:00