Untitled

disp('Vanishing line:');
disp(vanishingLine);

disp('Circumference matrix:');
disp(Circumference);



Circumference_norm = Circumference / norm(Circumference, 'fro');
disp('Cicornferenza normalizzata:');
disp(Circumference_norm);

vanishingLine_norm = vanishingLine / norm(vanishingLine);
disp('Vanishing line normalizzata:');
disp(vanishingLine_norm);

zr_norm = (-xr*vanishingLine_norm(1) - yr*vanishingLine_norm(2)) / vanishingLine_norm(3);
X_norm = [xr, yr, zr_norm];
new_eq_norm = simplify(X_norm * Circumference_norm * X_norm');
[sol_xr_norm,sol_yr_norm] = solve(new_eq_norm == 0, [xr, yr]);

disp('Soluzione con matrice normalizzata:');
disp(sol_xr_norm);
disp(sol_yr_norm);