clear all

tt = 8;
a  = 1;
b  = vpa(62.10,tt);
c  = 1;

b2 = vpa(b^2,tt);
ac = vpa(4 * vpa(a*c,tt), tt);
d1 = vpa(b2-ac,tt);
d1 = vpa(sqrt(d1),tt);
x1 = vpa(vpa(-b+d1,tt) / vpa(2*a,tt), tt);
x2 = vpa(vpa(-b-d1,tt) / vpa(2*a,tt), tt);

tt = 4;
a  = 1;
b  = vpa(62.10,tt);
c  = 1;

b2  = vpa(b^2,tt);
ac  = vpa(4 * vpa(a*c,tt), tt);
d1  = vpa(b2-ac,tt);
d1  = vpa(sqrt(d1),tt);
x10 = vpa(vpa(-b+d1,tt) / vpa(2*a,tt), tt);
x20 = vpa(vpa(-b-d1,tt) / vpa(2*a,tt), tt);

x11 = -vpa(vpa(2*c,tt) / vpa(b+d1,tt), tt);