Pascal [по методу симпсона вычислить значение функции

[qsccsq]

в данной программе нужно вычислить значение функции(если вы сюда зашли я думаю найдёте эту самую функцию в программе )
границы заданы [0,1] цикл For m:=1 to n-2 do (не уверен что правильно)

исправьте пожалуйста.

вот то что я смог написать

program name;
uses
  crt;
var
a,b,r,r2,eps :real;
n :integer;
function f(x:real):real;
  begin
    f:=x*exp(x)*sin(x);
  end;
function Simpson(a,b:real; n:integer):real;
var
  s,h:real;
  m,mn:integer;
 begin
  h:=(b-a)/(n);
  s:=f(a)+f(b);
  mn:=4;
 for m:=1 to n-2  do
  begin
  s:=s+mn*f(a+h*m);
  if (mn=4)
  then
  mn:=2
  else
  mn:=4;
  end;
  Simpson:=s*h/3;
 end;
Begin
{clrscr;}
writeln ('vvedite granicii a b');
Readln (a,b);
writeln('vvedite pogreshnost');
readln(eps);
n:=2;
r:=simpson(a,b,n);
repeat
 r2:=r;
 n:=n*2;
 r:=simpson(a,b,n);
until (abs(r-r2)<eps);
 Writeln;
 Writeln ('resultat po simpsony',r:6:6);
 writeln ('kol-vo razbienii* ',n,' ®в१Є®ў');
 readkey;
 end.

[k0tone]

Вот держи, подправил твой код:

program test;
uses crt;
var a,b,r,r2,eps:real;
    n:integer;

  function f(x:real):real;
    begin
      f:=x*exp(x)*sin(x);;
    end;

  function Simpson(a,b:real; n:integer):real;
  var s,h:real;
      m,mn:integer;
    begin
     h:=(b-a)/(2*n);
     s:=f(a)+f(b);
     mn:=4;
     for m:=1 to 2*n-1 do
       begin
         s:=s+mn*f(a+h*m);
         if mn=4 then
           mn:=2
         else
           mn:=4;
       end;
     Simpson:=s*h/3;
    end;

begin
writeln ('vvedite granicii a b');
Readln (a,b);
writeln('vvedite pogreshnost');
readln(eps);
n:=2;
r:=simpson(a,b,n);
repeat
 r2:=r;
 n:=n*2;
 r:=simpson(a,b,n);
until (abs(r-r2)<eps);
 Writeln;
 Writeln ('resultat po simpsony',r:6:6);
 writeln ('kol-vo razbienii* ',n,' ®в१Є®ў');
 readkey;
end.

Удачи тебе, уважаемый qsccsq!!!!!



Если использовать итерационную форму, то в вашей программе нужно поправить только цикл:

 for m:=1 to n-1 do