#include <u.h>#include <libc.h>#include "map.h" static struct coord stdpar;static double r0; staticXbonne(struct place *place, double *x, double *y){	double r, alpha;	r = r0 - place->nlat.l;	if(r<.001)		if(fabs(stdpar.c)<1e-10)			alpha = place->wlon.l;		else if(fabs(place->nlat.c)==0)			alpha = 0;		else 			alpha = place->wlon.l/(1+				stdpar.c*stdpar.c*stdpar.c/place->nlat.c/3);	else		alpha = place->wlon.l * place->nlat.c / r;	*x = - r*sin(alpha);	*y = - r*cos(alpha);	return(1);} projbonne(double par){	if(fabs(par*RAD) < .01)		return(Xsinusoidal);	deg2rad(par, &stdpar);	r0 = stdpar.c/stdpar.s + stdpar.l;	return(Xbonne);}