vg

Lignedeterre = Line((0, 0), (0, 1))
MUR: x = 0
A = (163.10995647187792, 148.5147305950029, 88.7810762092646)
B = (542.1445491879206, 642.3718317036137, 344.9271009060767)
h = Segment(A, B)
m = OrthogonalLine(B, SOL)
n = OrthogonalLine(A, SOL)
A_1 = Intersect(n, SOL)
B_1 = Intersect(m, SOL)
A1B1 = Segment(A_1, B_1)
PlanProjVertAB = Plane(B, A, A_1)
i = OrthogonalLine(A, MUR)
f = OrthogonalLine(B, MUR)
B_2 = Intersect(f, MUR)
A_2 = Intersect(i, MUR)
g = Segment(A_2, B_2)
PlanProjDeboutAB = Plane(B, A, A_2)
p = Plane(A, SOL)
j = IntersectPath(PlanProjVertAB, p)
C = Intersect(m, j)
α = Angle(B, A, C)
q = Plane(A, MUR)
k = IntersectPath(PlanProjDeboutAB, q)
D = Intersect(f, k)
β = Angle(B, A, D)