Conversion between Equatorial and Ecliptic Coordinates

[alpha,delta] = eclip2equat(lambda,beta) returns the right ascension and declination in the equatorial coordinate system corresponding to the longitudinal and latitudinal coordinates in the ecliptic coordinate system (lambda and beta).

[lambda,beta] = equat2eclip(alpha,delta) returns the longitudinal and latitudinal coordinates in the ecliptic coordinate system corresponding to coordinates of right ascension and declination (alpha and delta) in the equatorial coordinate system.

All inputs and outputs are in units of radians. longitudinal coordinates will be in the range of [-pi,pi] while latitudinal coordinates will be in the range [-pi/2, pi/2]. Inputs may be vectors but must be of equal size;