a=input('a=');
b=input('b=');
c=input('c=');
x1=0;
x2=0;
y1=0;
y2=0;
d=b^2-4*a*c;
if d<0
    fprintf('the roots are not real')
else
    x1=(-b+sqrt(d))/(2*a);
    x2=(-b-sqrt(d))/(2*a);
    y1=2*c/(-b-sqrt(d));
    y2=2*c/(-b+sqrt(d));
    x1
    x2
    y1
    y2
end
