Прошу помочь с составлением программ в Паскаль (см. внутри)

[BratStigmat]

Все просто, нужно составить программы для нахождения определенного интеграла способами прямоугольников, трапеций и симпсона. Функция: ((lg(x))^2)/2 в пределам от 1 до 4, кол-во частей разбиения 52. Убедительная просьба помочь кто чем может, от этого многое зависит)))

[lu]

а вы хотя бы половину написали?

[Asix]

На нашем сайте есть классная программа - решение интегралов, просмотрите ее код, в ней реализованы все ваши методы. Думаю вам не составит труда переделать программу под Паскаль =))

[BratStigmat]

Увы..исходный код программы  найти не удалось..

[BratStigmat]

вот мой код для метода прямоугольников:

program Integral;
var
a,b,h,s,int:real;
i,n:integer;

function F(x:real):real;
begin
F:=sqr(ln(x))/2;
end;
begin
a:=1;
b:=4;
n:=52;
h:=(b-a)/n;
s:=0;
for i:=1 to n-1 do
begin
s:=s+(F(a)+F(b))/2;
end;
int:=(b-a)/n*s;
writeln('Int=',int:3:3);
readln;
end.

[BratStigmat]

для трапеций:
program inttrap;
function F(x:real):real;
begin
F:=sqr(ln(x))/2;
end;
var
a,b,n,h,si,x:real;
i:integer;
begin

n:=52;
a:=1;
b:=4;
h:=(b-a)/n;
i:=0;
si:=0;
x:=a;
repeat
si:=si+(F(x)+F(x+h))/2*h;
x:=x+h;
i:=i+1;
until i>=n;
writeln ('Int=',si:3:3);
readln;
end.

[BratStigmat]

и наконец для симпсона:

program Simpson;
function F(x:real):real;
begin
F:=sqr(ln(x))/x;
end;
var
a,b,h,x:real;
n,i:integer;
int:real;

begin
a:=1;
b:=4;
n:=52;
h:=(b-a)/n;
int:=F(a)+F(b)+4*F(x+h);
for i:=1 to (n div 2)-1 do
begin
x:=a+2*h*i;
int:=int+2*F(x)*4*F(x+h);
end;
int:=h*int/3;
writeln ('Int=',int:2:3);
end.

[lu]

а ето что?

[BratStigmat]

это то что надо проверить и указать на ошибки, см. выше

[Данила]

есть такая замечательная вещь как дебаг...юзайте ее и сами смотрите в каком месте вычисления идут не туда и правьте

[Asix]

Наш главный программер вернулся =))
Ура =))

[Данила]

дада) я снова тут)