Среда turbo pascal. нужен алгоритм решения неполного квадратного уравнения y=ax^2+c

var
  a, c: real;
begin
  readln(a, c);
  a := sqrt(-c / a); //находим  
  writeln('x1 = ', a);
  writeln('x2 = ', -a);
end.