Incorrect angle detected between two planes

I want to calculate the angle between 2 planes, Reference plane and Plane1. I've attached the NX3 and MX3 matrices containing the 3D co-ordinates of both planes. When I feed this to the function plane_fit.m (by Kevin Mattheus Moerman), I get the coefficients:

reference_plane_coeff: [-0.13766204 -0.070385590 130.69409]

Plane1_coeff: [0.0044337390 -0.0013548643 95.890228]

Next, I find the intersection of both planes, separately on the XZ plane and get a line equation; ref_line_XZ and plane1_line_XZ respectively. I make the second coefficient 0.

Aref = reference_plane_coeff(1);
Cref = reference_plane_coeff(3);
ref_line_XZ = [Aref Cref];

Arun = Plane1_coeff(1);
Crun = Plane1_coeff(3);
plane1_line_XZ = [Arun Crun];

angle_XZ = acosd(  dot(ref_line_XZ,plane1_line_XZ ) / (norm(ref_line_XZ) * norm(plane1_line_XZ ))    )

I get the angle_XZ value as 0.0685 degrees.

When I plot these planes on a graph and view it, the angle seems to be much more than 0.0012 degrees. I'm talking about the angle made by the two lines after intersection of both planes with the XZ plane.

What am I doing wrong?

Also, when I tried to find angle between its normals, using:

angle_beta_deg = acosd(  dot(reference_plane_coeff,Plane1_coeff) / (norm(reference_plane_coeff) * norm(Plane1_coeff))    )

I got the angle as 0.0713.

On visual inspection of both planes' plots and manually calculating from the plot, angle_XZ should be around 9 degrees.

plane_fit.m (by Kevin Mattheus Moerman): http://in.mathworks.com/matlabcentral/fileexchange/22042-plane-fit/content/plane_fit.m