Операторы цикла (Pascal)

[Stacy31]

Задание: Выбрать алгоритм, составить блок-схему и программу для вычис-
ления в точках xi = a + i×h, i = 0,1,2…,n, h = (b – a)/n промежутка [a,b]
наибольшего и среднего значений функции y = f(x), указанной в варианте задания. Для выполнения задания воспользоваться оператором цикла
FOR.

 program U;
uses crt;
var p,v,x,y,a,b,h,la,max:real;
    g,d,n:word;
begin
 write('Vvedite chislo razbieniy(n): ');
     read(n);
     a:=1;
     b:=2;
     h:=(b-a)/n;
 la:=0;
 max:=0;
 for d:=0 to n do
    begin
    x:=a+d*h;
    writeln('x=',x:4:2);

        p:=exp(3*ln(x))/((x+1)*(x+2));
        v:=arctan((1-x)/sqrt(1-(1-x)*(1-x))/exp((1/3)*ln(1-ln(x))));
        y:=p+v;
    writeln('y=',y:4:2);
    la:=la+y;
    writeln('la=',la:4:2);
    if y>max then
      begin
      max:=y;
      writeln('max=',max:4:2);
       end;
       end;
    la:=la/n;
    writeln('srednee arif= ',la:7:5);
    writeln('max= ',max:7:5);
     readln;
     readln;

end.

Выдает ошибку деление на 0. Подскажите, что исправить?

[type_Oleg]

В формулах что-то наворочено.
Зачем городить такую жуть, как exp(3*ln(x)) , что нельзя написать простоx*x*x ? В крайнем случае, есть функция Power для всяких степеней.
И там почему-то арктангенс вместо арксинуса.

[Stacy31]

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

В формулах что-то наворочено...
И там почему-то арктангенс вместо арксинуса.

В Pascal же арксинус числа равен arcsin X = Arctan(X/Sqrt(l - X*X));
Вот поэтому там арктангенс. Или это не правильно?

[type_Oleg]

В Pascal есть и функция ArcSin в модуле Math
Поэтому можно так

// ...
uses crt, Math;
// и т д. 
 y:=x*x*x/(x+1)*(x+2)+ArcSin(1-x)/Power(1-Ln(x),1/3);
// и т .д.

А переменные p, v выбросить

[Stilet]

program Project1;

var p,v,x,y,a,b,h,la,max:real;
    g,d,n:word;

    function asin(x:real):real;
    var d:real;
    begin
      d:=sqrt(1-x*x);
      if d<>0 then
       asin:=arctan(x/d)
      else
       asin:=0;
    end;

begin
 write('Vvedite chislo razbieniy(n): ');
     n:=10;
     a:=1;
     b:=2;
     h:=(b-a)/n;
 la:=0;
 max:=0;
 for d:=0 to n do
    begin
    x:=a+d*h;
    writeln('x=':5,x:6:2);

        p:=exp(3*ln(x))/((x+1)*(x+2));
        v:=asin(1-x)/exp((1/3)*ln(1-ln(x)));
        y:=p+v;
    write('y=':5,y:6:2);
    la:=la+y;
    writeln('la=':10,la:4:2);
    if y>max then
      begin
      max:=y;
      writeln('max=',max:4:2);
       end;
       end;
    la:=la/n;
    writeln('srednee arif= ',la:7:5);
    writeln('max= ',max:7:5);
     readln;


end.

Годится?

[Stacy31]

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

В Pascal есть и функция ArcSin в модуле Math
Поэтому можно так

// ...
uses crt, Math;
// и т д. 
 y:=x*x*x/(x+1)*(x+2)+ArcSin(1-x)/Power(1-Ln(x),1/3);
// и т .д.

А переменные p, v выбросить

Спасибо большое. Все получилось! Урашечки!)))))))