#include <u.h>#include <libc.h>#include <stdio.h>#ifdef PLOT#include PLOT#else#include "iplot.h"#endif #ifndef PI#define PI	3.1415926535897932384626433832795028841971693993751#endif #define TWOPI (2*PI)#define RAD (PI/180)double	hypot(double, double);	/* sqrt(a*a+b*b) */double	tan(double);		/* not in K&R library */ #define ECC .08227185422	/* eccentricity of earth */#define EC2 .006768657997 #define FUZZ .0001#define UNUSED 0.0		/* a dummy float parameter */ struct coord {	float l;	/* lat or lon in radians*/	float s;	/* sin */	float c;	/* cos */};struct place {	struct coord nlat;	struct coord wlon;}; typedef int (*proj)(struct place *, float *, float *); struct index {		/* index of known projections */	char *name;	/* name of projection */	proj (*prog)(float, float);			/* pointer to projection function */	int npar;	/* number of params */	int (*cut)(struct place *, struct place *, float *);			/* function that handles cuts--eg longitude 180 */	int poles;	/*1 S pole is a line, 2 N pole is, 3 both*/	int spheroid;	/* poles must be at 90 deg if nonzero */};  proj	aitoff(void);proj	albers(float, float);int	Xazequalarea(struct place *, float *, float *);proj	azequalarea(void);int	Xazequidistant(struct place *, float *, float *);proj	azequidistant(void);proj	bicentric(float);proj	bonne(float);proj	conic(float);proj	cylequalarea(float);int	Xcylindrical(struct place *, float *, float *);proj	cylindrical(void);proj	elliptic(float);proj	fisheye(float);proj	gall(float);proj	globular(void);proj	gnomonic(void);int	guycut(struct place *, struct place *, float *);int	Xguyou(struct place *, float *, float *);proj	guyou(void);proj	harrison(float, float);int	hexcut(struct place *, struct place *, float *);proj	hex(void);proj	homing(float);proj	lagrange(void);proj	lambert(float, float);proj	laue(void);proj	loxodromic(float);	/* not in library */proj	mecca(float);proj	mercator(void);proj	mollweide(void);proj	newyorker(float);proj	ortelius(float, float);	/* not in library */int	Xorthographic(struct place *place, float *x, float *y);proj	orthographic(void);proj	perspective(float);int	Xpolyconic(struct place *, float *, float *);proj	polyconic(void);proj	rectangular(float);proj	simpleconic(float, float);int	Xsinusoidal(struct place *, float *, float *);proj	sinusoidal(void);proj	sp_albers(float, float);proj	sp_mercator(void);proj	square(void);int	Xstereographic(struct place *, float *, float *);proj	stereographic(void);int	Xtetra(struct place *, float *, float *);int	tetracut(struct place *, struct place *, float *);proj	tetra(void);proj	trapezoidal(float, float);proj	vandergrinten(void);proj	wreath(float, float);	/* not in library */ void	findxy(double, double *, double *);void	albscale(float, float, float, float);void	invalb(float, float, float *, float *); void	cdiv(double, double, double, double, double *, double *);void	cmul(double, double, double, double, double *, double *);void	csq(double, double, double *, double *);void	csqrt(double, double, double *, double *);void	ccubrt(double, double, double *, double *);double	cubrt(double);int	elco2(double, double, double, double, double, float *, float *);void	cdiv2(double, double, double, double, double *, double *);void	csqr(double, double, double *, double *); void	orient(float, float, float);void	latlon(float, float, struct place *);void	deg2rad(float, struct coord *);void	sincos(struct coord *);void	normalize(struct place *);void	invert(struct place *);void	norm(struct place *, struct place *, struct coord *);void	printp(struct place *);void	copyplace(struct place *, struct place *); int	picut(struct place *, struct place *, float *);int	ckcut(struct place *, struct place *, float);float	reduce(float); void	getsyms(char *);int	putsym(struct place *, char *, double, int);void	filerror(char *s, char *f);void	error(char *s);int	doproj(struct place *, int *, int *);int	cpoint(int, int, int);int	plotpt(struct place *, int);int	nocut(struct place *, struct place *, float *); extern int (*projection)(struct place *, float *, float *);