Quadratic Analytical Solution

X^4 + a3*X^3 + a2*X^2 + a1*X + a0 = 0

 

%%%% begin MATLAB code %%%%
a3 = 0;
a2 = -62;
a1 = 1;
a0 = 920;

T1 = -a3/4;
T2 = a2^2 - 3*a3*a1 + 12*a0;
T3 = (2*a2^3 - 9*a3*a2*a1 + 27*a1^2 + 27*a3^2*a0 - 72*a2*a0)/2;
T4 = (-a3^3 + 4*a3*a2 - 8*a1)/32;
T5 = (3*a3^2 - 8*a2)/48;

R1 = sqrt(T3^2 - T2^3);
R2 = (T3 + R1)^(1/3);
R3 = (1/12)*(T2/R2 + R2);
R4 = sqrt(T5 + R3);
R5 = 2*T5 - R3;
R6 = T4/R4;

x1 = T1 - R4 - sqrt(R5 - R6);
x2 = T1 - R4 + sqrt(R5 - R6);
x3 = T1 + R4 - sqrt(R5 + R6);
x4 = T1 + R4 + sqrt(R5 + R6);

y1 = 31-x1^2;
y2 = 31-x2^2;
y3 = 31-x3^2;
y4 = 31-x4^2;

x1, y1, x2, y2, x3, y3, x4, y4
%%%% end MATLAB code %%%