quadartic_function(16)

[x1,x2]=f(1,0,25)
function [x1,x2]=f(a,b,c)
    D=b*b-4*a*c;
    %%No real solutions
    if(D<0)
        x1=(-b+sqrt(-D)*i)/(2*a);
        x2=(-b-sqrt(-D)*i)/(2*a);
        %%x=-10:0.1:10;
    elseif(D==0)
        x1=-b/(2*a);
        x2=x1;
        %%x=-10+x1:0.1:10+x1;
    else
        x1=(-b+sqrt(D))/(2*a);
        x2=(-b-sqrt(D))/(2*a);
        %%x=-10+x2:0.1:10+x1;
    end
    x=-10:0.01:10;
    y=a*x.*x+b*x+c;
    plot(x,y);
    axis([-10 10 -20 20])
end