Паскаль решение интегралов

[Bordi]

Посмотрите пожалуйста, возможно ошибся в написание кода из формулы интеграла. Подскажите пожалуйста решение. Заранее спасибо. сам интеграл (arcsin(x)^2+1)/(sqrt(1-x^2)) метод трапеции.

 program Integr6;

uses crt;
var n,i:integer; a,b,h,s,x:real;

  {Подинтегральная функция}
  function f(x:real):real;
  begin f:=(arctan(x/sqrt(1-x*x)*(x/(1-x*x)))+1)/(sqrt(1-x*x)) end;

begin
  clrScr;
  a:=0; b:=sin(1);
  n:=10;
  {Ввод нижнего а и верхнего приеделов интегрирования b,
  а также кол-ва отрезков разбиения b}

  h:=(b-a)/n;
  x:=a;
  s:=2*f(x);
  for i:=1 to (n-1) do begin x:=x+h; s:=s+2*f(x); end;
  s:=h/2*(s+f(b));

  writeln (n,' ',s:0:6); {Вывод на экран результатов счета};
  readln
  end.          
или 
program Integr6;

uses crt;
var n,i:integer; a,b,h,s,x:real;

  {Подинтегральная функция}
  function f(x:real):real;
  begin f:=(x*x*arctan(x/sqrt(1-x*x)))+1)/(sqrt(1-x*x)) end;

begin
  clrScr;
  a:=0; b:=sin(1);
  n:=10;
  {Ввод нижнего а и верхнего приеделов интегрирования b,
  а также кол-ва отрезков разбиения b}

  h:=(b-a)/n;
  x:=a;
  s:=2*f(x);
  for i:=1 to (n-1) do begin x:=x+h; s:=s+2*f(x); end;
  s:=h/2*(s+f(b));

  writeln (n,' ',s:0:6); {Вывод на экран результатов счета};
  readln
  end.

[Bordi]

Интеграл к коду

[ForenLi]

Т.е. вы не видите разницы между arctan и arcsin? Между х и x/sqrt(1-x*x)*(x/(1-x*x)) и x/sqrt(1-x*x))? Откуда во втором коде взялся x^2?

[Bordi]

ForenLi, arcsin(X) = arctan(X/sqrt(1-sqr(X))), возведение в квадрат возможно ? я и спрашиваю поэтому. Спасибо

[ViktorR]

Возможно это будет через чур, но всё же.
Откройте для себя - Т.А. Павловская, Программирование на ЯВУ Паскаль (есть в разделе литературы).
Там есть прекрасный пример по оформлению интегрирующей функции и использованию параметра процедурного (функционального) типа.
В моём варианте книги стр.103.

[Bordi]

ViktorR, Спасибо за совет, я посмотрю! но все же, как интеграл записать ?

[ForenLi]

Сообщение от Bordi

arcsin(X) = arctan(X/sqrt(1-sqr(X)))

Ну если хотите так, то и идите: у вас (arcsin(X))^2, а значит arctan(X/sqrt(1-sqr(X))))*arctan(X/sqrt(1-sqr(X))

[Bordi]

ForenLi, огромное спасибо Вам! а как можно было еще написать ? если Вам не сложно значение arcsin?

[ForenLi]

В PascalABC.NET, Delphi arcsin уже есть. Если хотите сами, то лучше вынести в отдельную функцию с контролем http://pascal.net.ru/ArcSin

[Alex314]

Что за проблема? Подключите NumLibABC, там есть превосходный численный метод для вычисления интегралов. Конечно, это не метод трапеции... намного серьёзнее.