% symbolics
 
  clear
  syms a0 a1 a2 x

  % no numerical values to coefficients
  eqn1 = a2*x^2 + a1*x + a0;

  % MATLAB finds only symbolic solutions
  sols1 = solve(eqn1)
  pretty(sols1)

  % now give coefficients numerical values
  a0 =   4;
  a1 =  -4;
  a2 =   1;

  eqn2 = a2*x^2 + a1*x + a0;
  
  % MATLAB can find numerical roots now
  sols2 = solve(eqn2)

%------------------------------------------------------------------------------
% Output:
% sols1 =
% 
% [ 1/2/a2*(-a1+(a1^2-4*a2*a0)^(1/2))]
% [ 1/2/a2*(-a1-(a1^2-4*a2*a0)^(1/2))]
% 
% 
%                        [             2           1/2]
%                        [    -a1 + (a1  - 4 a2 a0)   ]
%                        [1/2 ------------------------]
%                        [               a2           ]
%                        [                            ]
%                        [             2           1/2]
%                        [    -a1 - (a1  - 4 a2 a0)   ]
%                        [1/2 ------------------------]
%                        [               a2           ]
% 
% sols2 =
% 
% [ 2]
% [ 2]
%------------------------------------------------------------------------------