Вычисление арккосинуса рядом тейлора

[isla]

Ребят, прошу помощи, всегда старался обойтись своими навыками(и вроде удавалось), но сейчас как-то странно, знаю в чем ошибка, но исправить не могу.

Задание:

вот мой код:

program laba3;
Const maxiter =2;
var arccos,rez, x,xnach, eps,xk,dx: real;
c,y,a: real;
n:integer;
done:boolean;
begin
writeln ('x, xk, dx, eps - ?');
readln (x, xk, dx, eps);

  done:=true;


writeln ('|    x   | |     y    |');
while abs(x)<1 do begin


     c:=1;
     y:=0;
     n:=0;



   while abs(c)>eps do begin
      c:=c*( (sqr(2*n+1)*sqr(x)*n) / ((n+1)*(2*n+3)*(2*n-1)) ) ;
      writeln ('', c);
      y:=y+c;
            inc(n);

          if n>maxiter then begin
          writeln ('ooops');
           done:=false;
            break
            end;
         end;



     
              if done then begin

                rez:=pi/2-(x+y);
                writeln (x:6:5,'     ',rez:6:4);
                end;



y:=0;
x:=x+dx;
end;

writeln ('|________| |__________|');
readln;
end.

________
Код нужно оформлять по правилам:
тегом [CODE]..[/СODE] (это кнопочка с решёточкой #)
Не забывайте об этом!
Модератор.

[Serge_Bliznykov]

1) а что за странная одинокая точка с запятой посредине программы?

2)

но сейчас как-то странно, знаю в чем ошибка, но исправить не могу.

ну так расскажите, в чём ошибка, что не работает, тогда больше шансов получить помощь...

[isla]

Прошу прошения у модератора, не нашёл кнопочку.

Проанализорав код, я понял, что у меня что-то с входными данными (c и y); x^3 как и все остальные нечетные степени не получаться никогда, как исправить?

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

Спасибо заранее

Убрал точку с запятой. Насчет того, что делал
1)Реккурентной формулой вывел множитель, если нужно могу рассписать рассчет.
2)Остальной алгоритм в учебнике, его выписывать не буду, думаю он и так понятен)

[Sweta]

Напишите исходную формулу.

[isla]

arccosX=Pi/2-(x+сумматор(по N) ((2n-1)*x^(2n+1))/(2n*(2n+1)))

[Sweta]

Похоже что Вы неправильно поняли формулу.
Подставляем n=1, получим ((1)*x^(3))/(2*3)
Подставляем n=2, получим (3)*x^(5))/(4*5)
Подставляем n=3, получим (5)*x^(7))/(6*7) , что не соответствует Вашему примеру
(у меня очень размытая картинка)
Очевидно в Вашей формуле (((2n-1)!)*x^(2n+1))/(2*n!*(2n+1))) и чтобы получить нечетную степень первое значение д.б. Х

[isla]

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

Похоже что Вы неправильно поняли формулу.
Подставляем n=1, получим ((1)*x^(3))/(2*3)
Подставляем n=2, получим (3)*x^(5))/(4*5)
Подставляем n=3, получим (5)*x^(7))/(6*7) , что не соответствует Вашему примеру
(у меня очень размытая картинка)
Очевидно в Вашей формуле (((2n-1)!)*x^(2n+1))/(2*n!*(2n+1))) и чтобы получить нечетную степень первое значение д.б. Х

Мммм, я тоже в Вики видел формулу с факториалами, но в моем учебнике без факториалов, и нужно сделать по данной формуле. могу сфотографировать еще раз

[Sweta]

Вопрос не в факториалах, а как Вы поняли запись формулы.
Если не нравятся факториалы, то можно рассмотреть как
arccosX=Pi/2-(x+сумматор(по N)((x^(2n+1))*Произведение(по N)(2n-1))/((2n+1)*(Произведение(по N)2n))

Множитель получится ((x^2)*((2n-1)^2))/(2*n*(2*n+1)), при начальном c=x

[Z1000000]

program arccos;
var
 x,y,sumodd,sumeven,e: real;
 farccosstep,farccos : real;
 n : Integer;
begin
write('Vvedite x=');
readln(x);
write('Vvedite tochnost e=');
readln(e);
n := 1; // øàã
sumeven := 1;
sumodd := 1;
farccos := Pi/2 - x;
farccosstep := x;
while ( farccosstep > e ) do
 begin
 sumeven := sumeven * (2*n -1);
 sumodd := sumodd * (2*n);
 farccosstep := (sumeven * Power(x,2*n+1)) / (sumodd * (2*n+1));
 inc(n);
 writeln(farccosstep);
 farccos := farccos - farccosstep;
 end;
writeln(farccos);
end.

[isla]

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

program arccos;
var
 x,y,sumodd,sumeven,e: real;
 farccosstep,farccos : real;
 n : Integer;
begin
write('Vvedite x=');
readln(x);
write('Vvedite tochnost e=');
readln(e);
n := 1; // øàã
sumeven := 1;
sumodd := 1;
farccos := Pi/2 - x;
farccosstep := x;
while ( farccosstep > e ) do
 begin
 sumeven := sumeven * (2*n -1);
 sumodd := sumodd * (2*n);
 farccosstep := (sumeven * Power(x,2*n+1)) / (sumodd * (2*n+1));
 inc(n);
 writeln(farccosstep);
 farccos := farccos - farccosstep;
 end;
writeln(farccos);
end.

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