Задача на итерационные циклы, вроде написал, а ответы дикие... Опции

[TheKnyazz]

Нужна помощь в решении задачи на итерациооные циклы
в общем вот мой код

 program lab4;
uses crt;
const maxiter = 500;
var
   x1,x2,dX,E:real; {x1 - Xnach, x2- Xkon, dX- shag, E -tochnost}
   n:integer; {schetchik vicheslennih elementov ryada}


function getY(x,E:real):real;
         var y,ch:real; {y-znachenie funcii arcctg(x), ch - element ryada}
         begin
             ch:= 1/x; {pervii element ryada}
             y:= -pi/2;
             n:=1;

             while (abs(ch)>E) do
               begin
                    y:=y+ch;
                    inc(n);
                    ch:=ch*(-1)*(2*n+1) /((2*n+3)*x*x); {vichislaem ocherednoi element   ryada}
                    {if (n mod 2) <> 0 then ch:=-ch;}
                    
                     {podschitivaem kol-vo vicheslennih elementov ryada}

                    if (n>maxIter) then
                       begin
                            write('ryad rashoditsa!');
                            break;
                       end;
               end;
             getY:=y;
         end;

begin
 clrscr;
 writeln('vvedite x1,x2,dX,E:');
 readln(x1,x2,dX,E);
 if (x1>=(-pi/2)) and (x1<=(pi/2)) and (x2>=(-pi/2)) and (x2<=(pi/2)) then {sootvetvuet oblasti dop. znachenii ? stroim tablicu}
    begin

      writeln('arcctg(x)------------------------------');
      writeln('|X                |Y                |N'); {shapka tablici}

      while (x1<=x2) do
         begin
              writeln('|',x1,'|',getY(x1,e),'|',n); {postroenie tablici znachenii ot x1 do x2 s shagom dX}
              x1:=x1+dX;
         end;
      writeln('---------------------------------------');

    end
 else
     writeln('nevernie dannie x1<=x2,x=[-pi/2;pi/2]!!');
begin
writeln('nazmite enter');
readln;
end;
end.

помогите пожалуйста. вот условие задачи.

[volvo]

1) отсутствует проверка на X = 0... Если уж ты пишешь, что X - из интервала -pi/2 .. +pi/2, то будь добр обрабатывать все значения, которые этот интервал попадают... При вводе 0 программа вылетает с делением на 0, что уже нехорошо...

2)

                    ch:=ch*(-1)*(2*n+1) /((2*n+3)*x*x); {vichislaem ocherednoi element   ryada}

неправда... На картинке представлено нечто другое... Что ты здесь делаешь?
Изначально: ch = 1/x
n = 1 -> ch = -3 / (x*x*x)
n = 2 -> ch = 3*5 / (x³*x*x) = 15/x⁵
(красным выделены предыдущие значения ch)

И так далее... А что должно быть?

Сообщение отредактировано: volvo - 26.10.2008 19:15