1
/* determine latitude from authalic latitude */
3
static const char SCCSID[]="@(#)pj_auth.c 4.3 93/06/12 GIE REL";
6
# define P00 .33333333333333333333
7
# define P01 .17222222222222222222
8
# define P02 .10257936507936507936
9
# define P10 .06388888888888888888
10
# define P11 .06640211640211640211
11
# define P20 .01641501294219154443
14
pj_authset(double es) {
17
if (APA = (double *)pj_malloc(APA_SIZE * sizeof(double))) {
30
pj_authlat(double beta, double *APA) {
32
return(beta + APA[0] * sin(t) + APA[1] * sin(t+t) + APA[2] * sin(t+t+t));