r. goldstein   july 1977
fills the cons array as follows
cons(1)=alpha0 in radians
cons(2)=sin(alpha0)
cons(3)=cos(alpha0)
cons(4)=delta0 in radians
cons(5)=sin(delta0)
cons(6)=cos(delta0)
cons(7)=n in radians
cons(8)=sin(n)
cons(9)=cos(n)
cons(10)=j in radians
cons(11)=sin(j)
cons(12)=cos(j)
cons(13)=i0 in radians
cons(14)=sin(i0)
cons(15)=cos(i0)
cons(16)=psi0 in radians
cons(17)=sin(psi0)
cons(18)=cos(psi0)
cons(19)=phi0 in radians
cons(20)=sin(phi0)
cons(21)=cos(phi0)
cons(22)=alpha0-n
cons(23)=sin(alpha0-n)
cons(24)=cos(alpha0-n)
cons(25)=omeg3
cons(26)=p
cons(27)=mu
cons(28)=-pdot/(2*p0)