a=input('a of ax**2 +bx+c ');   %s6dn.m
b=input('b of ax**2 +bx+c ');
c=input('c of ax**2 +bx+c ');
d=b*b-4*a*c ;
if (d==0) ;
  fprintf('%f x**2 + %f x + %f : Double real root  Delta=%f\n',a,b,c,d) ;
else
  if  (d<0) ;
    fprintf('%f x**2 + %f x + %f : No real root      Delta=%f\n',a,b,c,d) ;
  else
    fprintf('%f x**2 + %f x + %f : Two real roots    Delta=%f\n',a,b,c,d) ;
  end
end
