Matlab - porque cuando ejecuto este codigo no puedo poner variables negativas

a=input('¿Cuál es el valor de a?');

b=input('¿Cuál es el valor de b?');

c=input('¿Cuál es el valor de c?');

d=input('¿Cual es el valor de d?');

opc= "s"

x=d

while opc~="n"

    cal1=(((((a*x+b)*x)+c)*x)+d)

    cal2=((((a*3)*x+b*2)*x)+c)

    cal3=((6*a)*x)+(2*b)

    xn=x-((cal1*cal2)/((cal2*cal2)-(cal1*cal3)))

    x=xn

    if cal1==0  then opc="n"

    a1=a;

    b1=((a*xn)+b)

    c1=(((a*xn)+b)*xn+c)

    x2=(-b1+((b1*b1-4*a1*c1)^(1/2)))/(2*a1)

    x1=(-b1-((b1*b1-4*a1*c1)^(1/2)))/(2*a1)

disp (cal1, "F1=")

disp ('xn')

disp (xn)

disp ('x2')

disp (x2)

disp ('x1')

disp (x1)

liminf=xn-2

limsup=x2+2

    x=linspace(liminf,limsup);

    fun=a*(x.^3)+b*(x.^2)+(c*x)+d;

    y=evstr(fun);

    plot(x,y,'-r')

if cal1~=0 then opc="s"

end

end

end