1  function __angy_angz(p1, p2) = 
2      let(v = p2 - p1) 
3      [
4          atan2(v.z, norm([v.x, v.y])), 
5          atan2(v.y, v.x)
6      ];