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

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