2820.sagews

# A few points [msa_x, msa_y, wavelength].
p1w1 = [-106, -102, 1]
p1w2 = [-106, -102, 5]
p2w1 = [-106, 102, 1]
p2w2 = [-106, 102, 5]
p2w3 = [-106, 102, 3]
# Form a matrix.
A = Matrix([p1w1, p1w2, p2w1]) 
# For previously defined points these are actual detector coordinates.
d1w1 = [-100, -102]
d1w2 = [-70, -102]
d2w1 = [-140, 102]
d2w2 = [-110, 102]
d2w3 = [-130, 102]
# Form answers for x and y.
Bx = vector([d1w1[0], d1w2[0], d2w1[0]])
By = vector([d1w1[1], d1w2[1], d2w1[1]])
# Find solutions.
Rx = A.solve_right(Bx)
Ry = A.solve_right(By)
# Show solutions.
Rx, Ry

((255/212, -10/51, 15/2), (0, 1, 0))

A

[-106 -102 1] [-106 -102 1] [-106 -102 5] [-106 -102 5] [-106 102 3] [-106 102 3]

Bx, By

(-100, -102, -70, -102, -130, 102)

Rx, Ry

\newcommand{\Bold}[1]{\mathbf{#1}}\left(\left(\frac{255}{212},\,-\frac{10}{51},\,\frac{15}{2}\right), \left(0,\,1,\,0\right)\right)