I’d expect the output to be in radians. All the Degree multiplier does is convert the specified value form degrees to radians. Try typing Degree//N to get its numeric value. If you want to retain the units in the output, you could try wrapping the angular variables, x, in Quantity[x,”AngularDegrees”] - tip, you can generate a nicely formatted version of this Quantity expression by typing ctrl= followed by “x Angular Degrees”. This might return a Quantity result (in this case either radians or angular degrees). I’m not sure whether it would have the desired outcome, but might be worth a go to see what you get.
Sorry, haven’t been near my computer for a while. Having checked, I don’t think my suggestion works, I wasn’t sure that it would - see attached notebook. But the value of phi in your output is definitely in radians.
Thank you for this very nice demonstration.
I realised that, in my original function, some of the quantities are in Radians and some are in Degree. Therefore, my results were not consistent.