Для языка free pascal: вычислить значение суммы ряда y(x) в произвольной точке х , если y(x) = sinx + sin^2(x) + sin^3(x) + + sin^n(x) величину х вводить выполнить вычисления для 5,10,15 и 20 слагаемых

var x,y,p:real; n,i:integer;
begin
write('x = ');
readln(x);
write('n = ');
readln(n);
p:=1;
for i:=1 to n do
begin
p:=p*sin(x);
y:=y+p;
end;
writeln('y = ',y:6:4);
end.

Пример:
x = 1.5
n = 5
y = 4.9626