Вычисление определённого интеграла

[Gebbellz]

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

2*X-0.4*ln*(0.5*X+0.1);

Подинтегральная функция
x/0.5*X+0.1
Интервал интегрирования (a.b) - 3.5
количество разбиения частей n1.n2.n3 -50.500.1000


Можете помочь please?

[кусака]

http://vkontakte.ru/photo23922934_162968907

метод симпсона

[Gebbellz]

Сообщение от кусака

http://vkontakte.ru/photo23922934_162968907

метод симпсона

за метод симпсона спасибо огромное А первые два не знаешь:?)

uses crt;
var
a,b,h,f:real;
begin
crlscr;
writeln('a,b');readln(a);readln(b);
writeln('n');readln(n);
h;=(b-a)/n;
s:=0;
for i:=1 to n-1 do
s:=s+(a+h*i)/2*x-0.4*ln(0.5*x+0.1);
writeln(f:8:6);readkey;
end;end.

Помогите найти ошибку в первой методе(((

Что никто не знает?

[Alex_sim]

Переменные должны быть другого типа и в конце кода лучше будет если написать f:7:3:1

[Gebbellz]

Помогите ещё пожалуйста со вторым методом - прямоугольников...приминительно к моим данным

[Gebbellz]

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

uses crt;
var i, n: integer; a, b, h, x, s, c : real;
function f(x : real):real;
begin f:=(1/x)*sin(3.14*x/2); end;
begin
clrscr;
write('Введите нижний предел интегрирования '); readln(a);
write('Введите верхний предел интегрирования '); readln(b);
write('Введите количество отрезков '); readln(n);
h:=(b-a)/n; c:=a+h/2;
for i:=0 to n-1 do
begin x:=c+i*h; s:=s+f(x)*h; end;
write('Интеграл равен ',s:15:10); readln;
end.

метод трапеций

uses crt;
function f(x:real):real;
begin
f:=x/(0.5*x+0.1);
end;
function ff(x:real):real;
begin
ff:=2*x-0.4*ln(0.5*x+0.1);
end;
var
i,h:integer;
a,b,n,e,s,c,k:real;
begin
clrscr;
writeln ('a=');readln(a);writeln('b=');read ln(b);writeln('n=');readln(n);
h:=(b-a)/n;
e:=(f(a)+f(b))/2;
begin
s:=0
for i:=1 to n-1 do
s:=s+f(a+i*h);end;
c:=(s+e)/h;
writeln (c:5:10);
k:=ff(b)-ff(a);
writeln(c:5:10);
readkey;
end.

[Gebbellz]

метод симпсона...где ошибка помогите ((((

uses crt;
function f(x:real):real;
begin
f:=x/(0.5*x+0.1);end;
var
a,b,h,x:real;
i,n:integer;
begin
write ('a=');readln(a);
write('b=');readln(b);
writel('n=');readln(n);
if (n mod 2)>0 then begin n:=n+1
writeln ('chislo ne4etnoe zamenit` na n=',n);end;
h:=(b-a)/n;
s:=f(a)+f(b)+4*f(a+i*h);
for i:=1 to n-1 do
begin
x:=a+2*h*i;
s:=s+2*f(x)+4*f(x+h);end;
s:=h*s/3;
writeln('integral=',s);
end.