Y=x^3*tg^2*(x+1)^2+sqrt(1/x) решение для pascal abc

Var x, y: real;
begin readln(x);
          if (x<>0) then
          begin y:=x*x*x*sqr(tg(x))*sqr(x+1)+sqrt(1/x);
                    writeln('y=', y);
          end
          else writeln('нет решения');  {или, если вы деление на 0, то пишите у=∞}
readln;
end.