Thank you. I don't know how to use MathML, so I had to convert your equation to an Excel-friendly format. I finally got δ = asin(cosθ/(sqrt(a^2+b^2)))+((atan(b/a))) and plugging in the values of a, b and θ, I got δ = 27.1°. Using a slightly more accurate value of θ (60.19°), I got δ = 27.13°. The declination of the North Galactic Pole in Stellarium is 27°07'41.7" which converts to 27.13°.

Does b = sin 192.8583° sin 23.439289° mean (sin 192.8583°*sin 23.439289°)? That's where I'm getting messed up trying to duplicate your calculations.

Well, that's brilliant if you don't mind me saying so. I've been looking for this calculation for several years.