Интеграл

[Nobel]

народ помогите дописат код .. пожалуйста .
в кдое не работает второй цикл
[HTML]program Khan;
var a,b,dx,ST,SP,SS,k:real;
n,i,v,t:integer;
function f(x:real):real;
begin
f:=2*x+4;
end;
begin
writeln ('vvedite a,b,n');
readln (a,b,n);
SP:=0;
ST:=0;
SS:=0;
t:=0;
repeat
n:=n+1;
dx:=((b-a)/n);
For i:=1 to n do
begin
k:=a;
SP:=f(a)*dx+SP;
ST:=(((f(a)+f(a+dx))/2*dx)+ST);
SS:=((f(a)+4*f(a+dx/2)+f(a+dx))*dx)/6+SS ;
a:=k;
a:=a+dx;
end;
t:=t+1;
if t>800 then break;
until t>t+1
writeln ('SP=',SP:2:2,' ST=',ST:2:2,' SS=',SS:2:2);
readln;
end.

[McFrey]

что программа должна делать?выложи полный текст задания

[McFrey]

ну раз ты не хочешь выкладывать задание то за правильность не ручаюсь

program Khan;

var
  a, b, dx, ST, SP, SS: real;
  n, i, t: integer;

function f(x: real): real;
begin
  f := 2 * x + 4;
end;

begin
  writeln('vvedite a,b,n');
  readln(a, b, n);
  SP := 0;
  ST := 0;
  SS := 0;
  t := 0;
  repeat
    dx := ((b - a) / n);
    for i := 1 to n do
    begin
      SP := (f(a) * dx) + SP; 
      ST := ((f(a) + f(a + dx)) / 2 * dx) + ST; 
      SS := (((f(a) + 4 * f(a + dx / 2 * dx) + f(a + dx)) * dx) / 6 ) + SS;   
      a := a + dx;
    end;
    t := t + 1;
    n := n + 1;
    if t > 800 then break;
  until t > t + 1;
  writeln('SP=', SP:2:2, ' ST=', ST:2:2, ' SS=', SS:2:2);
  readln;
end.

и кстати "break" это отвратительно

[Nobel]

в интеграле надо найти оптимальное количество разбиений в котором все три методи равны... Праямоуголник, Триуголник и Симпсон..
т.е надо найти n

[McFrey]

как это оптимальное количество разбиений?...что означает "оптимальное" ? ведь чем больше количество разбиений, тем точнее результат. а в чем тогда должна выражаться оптимальность?

[Nobel]

условия зачаи у меня нету. Ну знаю что надо найти n в котором все 3 методы равны.

[McFrey]

интеграл можно вычислить этими 3 и еще несколькими методами, а n - число промежутков, на которые ты разбиваешь график функции для подсчета интеграла методом прямоугольников или треугольников или Симпсона...

[Nobel]

для все 3 метода . у меня в коде SP прямоуголник ST тиуголник SS симпсон

[McFrey]

я не могу понять....тебе нужно чтобы при определенном количестве разбиений n...результат интеграла для всех 3-х методов был равным?....тогда просто в своей программе n:=n+1 переставь перед t:=t+1 и подбери такое n при котором результат вычисления интеграла будет равным во всех 3-х случаях....например введи a=10,b=20,n=10000 и тогда значение интеграла будет примерно одинаковым во всех 3-х случаях.

[Nobel]

a=0 b=2
n должно найти программа ..
я вводил в ручную при n=800 все три методы =12 ..
но надо чтобы прграмма находыла n