Hard to answer, because it's unclear (to me) what it is you're trying to compute. Your 2 given planes, with normals N1 and N2, intersect the XY plane in 2 lines. Is it the angle between these lines that you want? If so, that angle is,
Nxy=[0 0 1];
N1=N1/norm(N1);
N2=N2/norm(N2);
Angle_XY=acosd(dot(cross(N1,Nxy), cross(N2,Nxy)));
Same sort of thing for the XZ and YZ planes.
Best Answer