Euler

cp = vector( [0.5,0.7,19.3] ) p1 = vector( [0.36292265447125238,15.739176957226150,19.290000794048506] ) p2 = vector( [-12.391110865165327,-6.8956218037985124,19.312672977586356] ) p3 = vector( [13.588187132562275,-6.6235568961199212,19.267327264346871] ) side1 = p1 - p2 side2 = p1 - p3 # Calc W w = side1.cross_product(side2) w = w/w.norm() # Calc U u = cp - p1 u = u/u.norm() # Calc V v = u.cross_product(w) v = v/v.norm() show( [u,v,w] ) show( asin( w[0] ).n() * 180 / 3.141592 ) show( atan2( w[1], w[2] ).n() ) show( atan2( v[0], u[0] ).n() ) 

\newcommand{\Bold}[1]{\mathbf{#1}}\left[\left(0.00911430339304340,\,-0.999958242852610,\,0.000664849442334250\right), \left(-0.999956951360168,\,-0.00911313117424658,\,0.00174535574594054\right), \left(0.001745261930295941,\,0.00001824711806909706,\,0.9999984768627587\right)\right]
\newcommand{\Bold}[1]{\mathbf{#1}}0.0999962143181992
\newcommand{\Bold}[1]{\mathbf{#1}}0.0000182471458599793
\newcommand{\Bold}[1]{\mathbf{#1}}-1.56168188342327

\newcommand{\Bold}[1]{\mathbf{#1}}\left[\left(0.00911430339304340,\,-0.999958242852610,\,0.000664849442334250\right), \left(-0.999956951360168,\,-0.00911313117424658,\,0.00174535574594054\right), \left(0.001745261930295941,\,0.00001824711806909706,\,0.9999984768627587\right)\right]
\newcommand{\Bold}[1]{\mathbf{#1}}0.0999962143181992
\newcommand{\Bold}[1]{\mathbf{#1}}0.0000182471458599793
\newcommand{\Bold}[1]{\mathbf{#1}}-1.56168188342327