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)