La fórmula es correcta.
Tal vez te has equivocado en alguna operación.
Para verificar puedes también utilizar la formula
cos(90 - a) = cos(90 - d).cos(90 - l) + sin (90-d).sin(90-l).cos(P)
que es lo mismo que
sin(a)=sin(d).sin(l)+cos(d).cos(l).cos(P)
Depejando cos(P) obtienes
cos(P) = (sin(a) - sin(d).sin(l)) / ( cos(d).cos(l) )
Una vez calculado la derecha, haces el inverso de cos
para obtener el ángulo P en el polo.
