Определенный интеграл Опции

[PUMA]

Для заданных границ интегрирования a и b вычислите значение определенного интеграла следующего вида (используя рекурсию)
как применить рекурсию? помогите пожалуйста

  program r;
     uses crt;
     var a,b,s,ss,dx,x,p1,p2,p3,p4,p5,p:real;
     i,n,m:integer;
     begin
     clrscr;
     write ('a=');
     readln (a);
     write ('b=');
     readln (b);
     write ('m>0, m=');
     readln (m);
     if a<b then dx:=(b-a)/m else dx:=(a-b)/m;
     x:=a;
     ss:=0;
     write ('n=0,1,2,3,...., n=');
     readln (n);
     if n=0 then for i:=1 to m do ss:=ss+exp(a*x)/a
            else if n=1 then for i:=1 to m do begin
                                              p1:=exp(a*x);
                                              p2:=a*sin(b*x);
                                              p3:=b*cos(b*x);
                                              ss:=ss+p1*(p2-p3)/(a*a+sqr(n*b))
                                              end
                        else while n>=1 do begin
                               p4:=n*(n-1)*b*b/(a*a+sqr(n*b));
                               s:=0;
                               for i:=1 to m do begin
                                                p1:=exp(a*x);
                                                p2:=a*sin(b*x);
                                                p3:=b*cos(b*x);
                                                p5:=exp((n-1)*ln(sin(b*x)));
                                                s:=s+p1*p5*(p2-p3)/(a*a+sqr(n*b));
                                                end;
                               ss:=ss+n*(n-1)*b*b*s/(a*a+sqr(n*b));
                               n:=n-2;
                               end;
     writeln ('int=',ss*dx:0:5);
     readln;
     end.

Правильно ли я вообще решила

Сообщение отредактировано: PUMA - 19.06.2010 18:23

[PUMA]

function IntPow(x: double; n: integer): double;
var
  r: double;i:byte;
begin
  r:=1;
   for i:=1 to n do
   r:=r*x;
  IntPow:=r
end;
function IntEAXSINNBX(a,b,x: double; n: integer): double;
begin
  case n of
    0: IntEAXSINNBX:= Exp(a*x)/a;
    1: IntEAXSINNBX:= Exp(a*x)*(a*Sin(b*x)-b*Cos(b*x))/(a*a+b*b);
    else IntEAXSINNBX:=
      (Exp(a*x)*IntPow(Sin(b*x),n-1)*(a*Sin(b*x)-n*b*Cos(b*x)) +
      n*(n-1)*b*b*IntEAXSINNBX(a,b,x,n-2))/
      (a*a+n*n*b*b)
  end
end;
var a,b,n,x:integer;
begin
writeln('vvedite a i b');
readln(a,b);
writeln('vvedite stepen n i chislo x');
readln(n,x);
writeln( IntEAXSINNBX(a,b,x,n));
readln;
end.

вроде как решает

[Lapp]

вроде как решает

Вау! респект тебе и твоим преподпвателям: ничего лишнего, никаких идиотских CRT, и даже печать результата непосредственно из функции без всяких перекладываний - уважаю! ))
Только надо бы сдвиги делать нормальные.. иначе увязнешь на второй странице кода..

Очень просто:
- после begin (case, repeat, while) все сдвигаем право на 2 поз
- перед end (until) - влево на 2 поз
- между функциями оставляй строку (лучше две)
- разделы декларвций тоже двигаем

Вот, посмотри и старайся так делать

function IntPow(x: double; n: integer): double;
var
  r: double;
  i:byte;
begin
  r:=1;
  for i:=1 to n do r:=r*x;
  IntPow:=r
end;


function IntEAXSINNBX(a,b,x: double; n: integer): double;
begin
  case n of
    0: IntEAXSINNBX:= Exp(a*x)/a;
    1: IntEAXSINNBX:= Exp(a*x)*(a*Sin(b*x)-b*Cos(b*x))/(a*a+b*b);
    else IntEAXSINNBX:=
      (Exp(a*x)*IntPow(Sin(b*x),n-1)*(a*Sin(b*x)-n*b*Cos(b*x)) +
      n*(n-1)*b*b*IntEAXSINNBX(a,b,x,n-2))/
      (a*a+n*n*b*b)
  end
end;


var
  a,b,n,x: integer;

begin
  writeln('vvedite a i b');
  readln(a,b);
  writeln('vvedite stepen n i chislo x');
  readln(n,x);
  writeln( IntEAXSINNBX(a,b,x,n));
  readln;
end.

Только сейчас заметил, что ты переделала IntPow .. Ага, я ошибся.. забыл вычитать n. Нужно было вот так:

function IntPow(x: double; n: integer): double;
var
  r: double;
begin
  r:=1;
  if n<0 then begin
    n:=-n;
    x:=1/x
  end;
  while n>0 do begin
    r:=r*x;
    Dec(n)
  end;
  IntPow:=r
end;

Молодец! )) Твоя функция, правда, не будет работать с отрицательными n, но оно тут и не нужно.