matplotlib Code

Brought to you by: cjgohlke, dsdale, efiring, heeres, and 8 others

[r8871]: / trunk / toolkits / basemap / src / pj_qsfn.c Maximize Restore History

17 lines (14 with data), 337 Bytes

/* determine small q */
#include <math.h>
#include <projects.h>

# define EPSILON 1.0e-7
	double
pj_qsfn(double sinphi, double e, double one_es) {
	double con;

	if (e >= EPSILON) {
		con = e * sinphi;
		return (one_es * (sinphi / (1. - con * con) -
		   (.5 / e) * log ((1. - con) / (1. + con))));
	} else
		return (sinphi + sinphi);
}