ВЫЧИСЛЕНИЕ ОПРЕДЕЛЕННЫХ ИНТЕГРАЛОВ (метод трапеций)

[maliyusha]

Здравствуйте,программисты! Вы не могли бы мне помочь с решением следующей задачи: даны подынтегральная функция (x ⋅ sin x2) и пределы интегрирования (0; 1),нужно вычислить определенный интеграл методом трапеций с заданной точностью вычислений Eps. Алгоритм решения знаю,решала многократно,но программа все равно не работает. Пожалуйста,напишите решение с комментариями,заранее спасибо
P.S.: необходимо к пятнице (26 декабря).

[eyo]

Выложи свой вариант решения, сча посмотрю что не так.

[maliyusha]

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

Выложи свой вариант решения, сча посмотрю что не так.

А как это сделать? Пыталась через вложения не выходит!

[maliyusha]

Program Integral;
const a=0;b=1;
var
dx,c,c1,s,s3,s2,e:real;
n,i:integer;
function f(x:real):real;
begin
f:=x*sin(sqr(x));
end;
begin
writeln('vvedite n');
readln(n) ;
writeln('vvedite tochnost');
readln(e);
s2:=0;
repeat
dx:=(b-a)/n;
for i:=2 to n do begin
c:=a+(i-1)*dx;
c1:=a+i*dx;
s3:=(f(c)+f(c1))*0.5*dx;
s2:=s2+s3;
n:=2*n;
end;
until abs(s3-s2)<e;
writeln(n,' ',s2:7:4);
readln;
end.

[maliyusha]

Получилось

[smileman]

maliyusha,странно,запустил ваш листинг,но ни каких ошибок не вышло,или программа не работает так,как вам нужно?

[maliyusha]

Понимаете,после ввода точности программа должна выводить ответ ('S2= '). Но у меня не выводит. Как исправить я,к сожалению, не знаю.

[smileman]

Program Integral;
const a=0;b=1;
var
dx,c,c1,s,s3,s2,e:real;
n,i:integer;
function f(x:real):real;
begin
f:=x*sin(sqr(x));
end;
begin
writeln('vvedite n');
readln(n) ;
writeln('vvedite tochnost');
readln(e);
s2:=0;
repeat
dx:=(b-a)/n;
for i:=2 to n do begin
c:=a+(i-1)*dx;
c1:=a+i*dx;
s3:=(f(c)+f(c1))*0.5*dx;
s2:=s2+s3;
n:=2*n;
end;
until abs(s3-s2)<e;
writeln('n=',n,' ','s2= ',s2:7:4);
readln;
end.

[maliyusha]

Нет,все равно курсор после точности останавливается и мигает,но на клавиши клавиатуры не реагирует. Приходится каждый раз закрывать программу. А у вас программа высчитывает окончательное значение?

[smileman]

vvedite n
3
vvedite tochnost
2
n= 2 s2= 0.2420

вот что у меня выводит,это не то?