Matlab - como programo la ecuacion cuadratica

%Ecuaicon cuadratica

p=1

while p==1

    a=input('Dame le valor de a:')

    b=input('Dame le valor de b:')

    c=input('Dame le valor de c:')

    s=(4*a*c)

    k=(b*b)

    if a<=0

    else

        if s>k

        else

        xu=(-b+( sqrt((k)-(s))))/(2*a)

        xd=(-b-( sqrt((k)-(s))))/(2*a)

        fprintf('X1 = %d',xu)

        fprintf('X2 = %d',xd)

        ecu=(

    xe=(a*(xu*xu)+b*xu+c)

%Script file: pregunta4.m

%Purpose:

%Record of revisions:

%     Date                 Programmer             Description of change

%  ==========        =======================      =====================

%  13/11/2014         GILMER CHALLCO MAMANI           

%Group: G18

%Define variables:

clear all;

clc;

global A B C X1 X2 X3 X4 D

disp('Ecuación de segundo grado A*X^2+B*X+C=0')

A=input('Ingrese el valor de A: ');

B=input('Ingrese el valor de B: ');

C=input('Ingrese el valor de C: ');

D=B^2-4*A*C;

if (D>=0)

    X1=(-B+D^0.5)/(2*A);

    X2=(-B-D^0.5)/(2*A);

    fprintf('La primera solución de la ecuación cuadratica es %2.2f\n',X1);

    fprintf('La segunda solución de la ecuación cuadratica es %2.2f\n',X2);

elseif D==0

    X1=(-B)/(2*A);

    X2=(-B)/(2*A);

    fprintf('La primera solución de la ecuación cuadratica es %2.2f\n',X1);

    fprintf('La segunda solución de la ecuación cuadratica es %2.2f\n',X2);

else

    X1=(-B)/(2*A);

    X2=(-B)/(2*A);

    X3=(abs(D))^0.5/(2*A);

    X4=-(abs(D))^0.5/(2*A);

    fprintf('La primera solución de la ecuación cuadratica es %2.2f + %2.2f i\n',X1,X3);

    fprintf('La segunda solución de la ecuación cuadratica es %2.2f + %2.2f i\n',X2,X4);

end

syms x

f=5*x^2+8*x+1

ezplot(f)