it acually is an assignment not really “homework” and i didnt post my “effort” cause all the wording is in a foreign language and it would complicate things. all the results are correct, just in the wrong order. im posting the order i want to achieve in the screenshot below. I hope you can help.

c = 152.34
omega=0
phi=0
kappa=0
import math
r11= math.cos(phi)*math.cos(kappa)
r12= math.sin(omega)*math.sin(phi)*math.cos(kappa) + math.cos(omega)*math.sin(kappa)
r13= -math.cos(omega)*math.sin(phi)*math.cos(kappa)+math.sin(omega)*math.sin(kappa)
r21=-math.cos(phi)*math.sin(kappa)
r22=-math.sin(omega)*math.sin(phi)*math.sin(kappa) + math.cos(omega)*math.cos(kappa)
r23=math.cos(omega)*math.sin(phi)*math.sin(kappa) + math.sin(omega)*math.cos(kappa)
r31=math.sin(phi)
r32=-math.sin(omega)*math.cos(phi)
r33= math.cos(omega)*math.cos(phi)
Xo= (6050.75+6777.80+6450.20+6102.55+6653.40)/5
Yo= (12424.50+12450.45+12050.10+11407.65+11475.80)/5
Zo=944.343
X= [6050.75, 6777.80, 6450.20, 6102.55,6653.40]
Y= [12424.50,12450.45,12050.10,11407.65,11475.80]
Z= [201.70,257.6,169.2,181.3,205.8]
A=
B=
C=
D=
E=
for i in range(5):
A1=r11*(X[i]-Xo) +r12*(Y[i]-Yo) +r13*(Z[i]-Zo)
A.append(A1)
A2=r21*(X[i]-Xo) +r22*(Y[i]-Yo) +r23*(Z[i]-Zo)
B.append(A2)
P=r31*(X[i]-Xo)+r32*(Y[i]-Yo)+r33*(Z[i]-Zo)
C.append(P)
x=-c*(A1/P)
D.append(x)
y=-c*(A2/P)
E.append(y)
print("Α1 ",A)
print(“Α2 :”,B)
print(“Π:”,C)
print("x ",D)
print(" y ",E)
F=
for i in range(5):
a1=(c/C[i]**2)*(-A[i]*r31+C[i]*r11)
F.append(a1)
G=
for i in range(5):
a2=(c/C[i]**2)*(-A[i]*r32+C[i]*r12)
G.append(abs(a2))
H=
for i in range(5):
a3=(c/C[i]**2)*(-A[i]*r33+C[i]*r13)
H.append(a3)
I=
for i in range(5):
a4=(-c/C[i]2)(C[i](-r13(Y[i]-Yo)+r12(Z[i]-Zo))-A[i](-r33(Y[i]-Yo)+r32*(Z[i]-Zo)))/636620
I.append(a4)
f11=-math.sin(phi)*math.cos(kappa)
f12=math.sin(omega)*math.cos(phi)*math.cos(kappa)
f13=-math.cos(omega)*math.cos(phi)*math.cos(kappa)
f21=math.sin(phi)*math.sin(kappa)
f22=-math.sin(omega)*math.cos(phi)*math.sin(kappa)
f23=math.cos(omega)*math.cos(phi)*math.sin(kappa)
f31=math.cos(phi)
f32=math.sin(omega)*math.sin(phi)
f33=-math.cos(omega)*math.sin(omega)
J=
for i in range(5):
a5=(-c/C[i]2)(C[i](f11(X[i]-Xo)+f12(Y[i]-Yo)+f13*(Z[i]-Zo))-A[i](f31(X[i]-Xo)+f32*(Y[i]-Yo)+f33*(Z[i]-Zo)))/636620
J.append(a5)
K=
for i in range(5):
a6=((-c/C[i])(r21(X[i]-Xo)+r22*(Y[i]-Yo)+r23*(Z[i]-Zo)))/636620
K.append(a6)
L=
for i in range(5):
b1=(c/C[i]**2)*(-B[i]*r31+C[i]*r21)
L.append(b1)
M=
for i in range(5):
b2=(c/C[i]**2)*(-B[i]*r32+C[i]*r22)
M.append(b2)
N=
for i in range(5):
b3=(c/C[i]**2)*(-B[i]*r33+C[i]*r23)
N.append(b3)
O=
for i in range(5):
b4=(-c/C[i]2)(C[i](-r23(Y[i]-Yo)+r22(Z[i]-Zo))-B[i](-r33(Y[i]-Yo)+r32*(Z[i]-Zo)))/636620
O.append(b4)
P=
for i in range(5):
b5=(-c/C[i]2)(C[i](f21(X[i]-Xo)+f22(Y[i]-Yo)+f23*(Z[i]-Zo))-B[i](f31(X[i]-Xo)+f32*(Y[i]-Yo)+f33*(Z[i]-Zo)))/636620
P.append(b5)
Q=
for i in range(5):
b6=((-c/C[i])(-r11(X[i]-Xo)-r12*(Y[i]-Yo)-r13*(Z[i]-Zo)))/636620
Q.append(b6)
print(F)
print(G)
print(H)
print(I)
print(J)
print(K)
print(L)
print(M)
print(N)
print(O)
print(P)
print(Q)
import numpy as np
xa=[-73.206, 65.578, -2.616,-70.080, 34.277]
xa = np.array(xa)
xa
ya=[80.153, 94.446, 10.861,-106.349,-103.873]
dl1=xa[0]-D[0]
dl2=ya[0]-E[0]
dl3=xa[1]-D[1]
dl4=ya[1]-E[1]
dl5=xa[2]-D[2]
dl6=ya[2]-E[2]
dl7=xa[3]-D[3]
dl8=ya[3]-E[3]
dl9=xa[4]-D[4]
dl10=ya[4]-E[4]
dl=[dl1,dl2,dl3,dl4,dl5,dl6,dl7,dl8,dl9,dl10]
print(" δλ :",dl)
import array
#Α
Ai=[F,G,H,I,J,K,L,M,N,O,P,Q]
np.set_printoptions(suppress=True)
array6 = np.array(Ai)
array6=np.reshape(Ai,(10,6))
print(“A=”,arrayFromAi)