Convertion de projection en ellipsoïdiques

[micmonay]

Bonjours

Je cherche à créer un programme de conversion des coordonnées suisse en coordonnées ellipsoïdiques. J'ai suivit la docs de ce lien http://geomatics.ladetto.ch/swiss_projection_fr.pdf en appliquant les formules mais je n'arrive pas au même résultat pour S impossible, je tourne la formule dans tous les sens mais je trouve pas.j'ai essayer d'écrire la formule en changent les variables par rapport aux données plus haut mais j'obtiens pas le même résultât. Peut être qu'il y a quelque chose que j'ai mal interpréter, par ex. S= ln tan ... pour moi s'est ln(tan(...)).

Merci d'avance.

PS : Je ne sais pas si je suis au bonne endroits dite moi si c'est une erreur