/*  * This file is part of the UCB release of Plan 9. It is subject to the license * terms in the LICENSE file found in the top-level directory of this * distribution and at http://akaros.cs.berkeley.edu/files/Plan9License. No * part of the UCB release of Plan 9, including this file, may be copied, * modified, propagated, or distributed except according to the terms contained * in the LICENSE file. */ #include <u.h>#include <libc.h>#include "map.h" static intXmollweide(struct place *place, double *x, double *y){	double z;	double w;	z = place->nlat.l;	if(fabs(z)<89.9*RAD)		do {	/*newton for 2z+sin2z=pi*sin(lat)*/			w = (2*z+sin(2*z)-PI*place->nlat.s)/(2+2*cos(2*z));			z -= w;		} while(fabs(w)>=.00001);	*y = sin(z);	*x = - (2/PI)*cos(z)*place->wlon.l;	return(1);} projmollweide(void){	return(Xmollweide);}