J'ai fait les calculs dans un tableur :
N ; p ; f ; d ; tmax
4 ; 6,4 ; 24 ; 45 ; =(19,3*B2+13*A2+0,1*C2)/(C2*COS(RADIANS(D2)))
Le tout est de bien penser à convertir la déclinaison en radians, unité d'angle utilisée par la fonction cosinus dans les tableurs