A = (0.3, 2.88)
B = (1.18, 2.94)
C = (1.56, 2.44)
D = (1.86, 1.66)
E = (2.34, 1.8)
F = (2.98, 1.5)
G = (3.46, 1)
H = (3.98, 0.56)
a = FitLineY({A, B, C, D, E, F, G, H})
I: (A + B + C + D + E + F + G + H) / 8
b = Variance({x(A), x(B), x(C), x(D), x(E), x(F), x(G), x(H)})
c = Variance({y(A), y(B), y(C), y(D), y(E), y(F), y(G), y(H)})
d = Covariance({A, B, C, D, E, F, G, H})
matrice1: {{b, d}, {d, c}}
M: matrice1^(-1)
J = ApplyMatrix(M, (1, 0))
K = ApplyMatrix(M, (0, 1))
e: ((x - x(I))^(2) * x(J)) + ((((x - x(I)) * (y - y(I))) * 2) * y(J)) + ((y - y(I))^(2) * y(K)) = 4
f = FirstAxis(e)