/*  * 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 "astro.h" voidsun(void){	double mven, merth, mmars, mjup, msat;	double dmoon, mmoon, gmoon;	double pturbb, pturbl, pturbr, lograd; 	ecc = .01675104 - 4.180e-5 * capt - 1.26e-7*capt2;	incl = 0;	node = 0;	argp = 281.220833 + .0000470684*eday + .000453*capt2		 + .000003*capt3;	mrad = 1;	anom = 358.475845 + .9856002670*eday - .000150*capt2		 - .000003*capt3;	motion = .9856473354; 	dmoon = 350.737681+12.1907491914*eday-.001436*capt2;	gmoon = 11.250889 + 13.2293504490*eday - .003212*capt2;	mmoon = 296.104608 + 13.0649924465*eday + 9.192e-3*capt2;	mven  = 212.448 + 1.602121635*eday;	merth = 358.476 + 0.985600267*eday;	mmars = 319.590 + .524024095*eday;	mjup = 225.269 + .083082362*eday;	msat  = 175.593 + .033450794*eday; 	dmoon = fmod(dmoon, 360.)*radian;	gmoon = fmod(gmoon, 360.)*radian;	mmoon = fmod(mmoon, 360.)*radian;	mven  *= radian;	merth *= radian;	mmars *= radian;	mjup *= radian;	msat  *= radian; 	icosadd(sunfp, suncp);	anom += cosadd(4, mmars, merth, mven, mjup)/3600.;	anom += sinadd(5, mmars, merth, mven, mjup, .07884*capt)/3600.; 	incl *= radian;	node *= radian;	argp *= radian;	anom = fmod(anom, 360.)*radian; /* *	computation of elliptic orbit */ 	lambda = anom + argp; 	pturbl = (6910.057 - 17.240*capt - 0.052*capt2)*sin(anom)		 + (72.338 - 0.361*capt) * sin(2.*anom)		 + (1.054 - 0.001*capt) * sin(3.*anom)		 + 0.018 * sin(4.*anom); 	lambda += pturbl*radsec; 	beta = 0.; 	lograd = (30.57e-6 - 0.15e-6*capt)		 - (7274.12e-6 - 18.14e-6*capt - 0.05e-6*capt2)*cos(anom)		 - (91.38e-6 - 0.46e-6*capt) * cos(2.*anom)		 - (1.45e-6 - 0.01e-6*capt) * cos(3.*anom)		 - 0.02e-6 * cos(4.*anom); 	pturbl = cosadd(5, mmars, merth, mven, mjup, msat);	pturbl += sinadd(3, dmoon, mmoon, merth) + .9;	pturbl *= radsec; 	pturbb =  cosadd(3, merth, mven, mjup);	pturbb += sinadd(3, gmoon, mmoon, dmoon);	pturbb *= radsec; 	pturbr =  cosadd(5, mmars, merth, mven, mjup, msat);	pturbr += cosadd(3, dmoon, mmoon, merth); 	lambda += pturbl;	if(lambda > pipi)		lambda -= pipi; 	beta += pturbb; 	lograd = (lograd+pturbr) * 2.30258509;	rad = 1 + lograd * (1 + lograd * (.5 + lograd/6)); 	motion *= radian*mrad*mrad/(rad*rad); 	semi = 961.182;	if(flags['o'])		semi = 959.63;	mag = -26.5;}