<pre>&lt;Cfinclude template=&quot;udflib.cfm&quot; /&gt;

&lt;cfoutput&gt;
&lt;html&gt;
&lt;head&gt;
&lt;title&gt;Earth/Sun&lt;/title&gt;
&lt;style type=&quot;text/css&quot;&gt;
pre {
	border: 1px dashed ##999;
	background-color: ##f7f7f8;
	margin: 10px 0px;
	padding: 5px 1%;
	clear: both;
}
div.row {
	clear: both;
	padding: 0px;
	margin: 0px;
}
pre.half {
	clear: none;
	float: left;
	width: 47%;
}
pre.half:last-child {
	float: right;
}
pre.third {
	clear: none;
	float: left;
	width: 30%;
}
pre.third:first-child {
	float: left;
	margin-right: 1.5%;
	width: 30%;
}
pre.third:last-child {
	float: right;
	margin-left: 1.5%;
	width: 30%;
}
table.pre th,
table.pre td {
	font-family: monospace;
	white-space: pre;
}
&lt;/style&gt;
&lt;link rel=&quot;stylesheet&quot; href=&quot;/s/table_data.css&quot; type=&quot;text/css&quot; /&gt;
&lt;/head&gt;
&lt;body&gt;
&lt;cfscript&gt;
startTime = getTickCount();
attributes.restart = structKeyExists(attributes, &quot;restart&quot;) AND isBoolean(attributes.restart) AND attributes.restart;
attributes.iterate = structKeyExists(attributes, &quot;iterate&quot;) AND isBoolean(attributes.iterate) AND attributes.iterate;

variables.my.dsn = &quot;sun_earth&quot;;
CHART_INTERVAL = 60; // s (don't change this)
BUILD_INTERVAL = 3 * 60 * 60; // s

Me = 5.974 * 10 ^ 24; // kg
Ms = 1.989 * 10 ^ 30; // kg
r = 1.496 * 10 ^ 11; // m
G = 6.673 * 10 ^ -11; // m^3 * k^-1 * s^-2 == m^3 / (kg * s^2)
&lt;/cfscript&gt;
&lt;cffunction name=&quot;cache&quot;&gt;
	&lt;cfargument name=&quot;t&quot; /&gt;
	&lt;cfargument name=&quot;Det&quot; /&gt;
	&lt;cfargument name=&quot;Dst&quot; /&gt;
	&lt;cfset var Vet = V_e(t) /&gt;
	&lt;cfset var Vst = V_s(t) /&gt;
	&lt;cfset var Aet = A_e(t) /&gt;
	&lt;cfset var Ast = A_s(t) /&gt;
	&lt;cfset var Dt = Dst + Det /&gt;
	&lt;cfset var rt = r - Dt /&gt;
	&lt;cfset var Ft = F(Me, Ms, rt) /&gt;
	&lt;cfset var AetF = a(Ft, Me) /&gt;
	&lt;cfset var AstF = a(Ft, Ms) /&gt;
	&lt;cfquery datasource=&quot;#variables.my.dsn#&quot; name=&quot;&quot;&gt;
		insert delayed into point
			(t, Dt, rt, Ft, AetF, AstF,
			 Det, Dst, Vet, Vst, Aet, Ast
			)
		values
			(&lt;cfqueryparam cfsqltype=&quot;cf_sql_integer&quot; value=&quot;#t#&quot; /&gt;,
			 &lt;cfqueryparam cfsqltype=&quot;cf_sql_double&quot; value=&quot;#Dt#&quot; /&gt;,
			 &lt;cfqueryparam cfsqltype=&quot;cf_sql_double&quot; value=&quot;#rt#&quot; /&gt;,
			 &lt;cfqueryparam cfsqltype=&quot;cf_sql_double&quot; value=&quot;#Ft#&quot; /&gt;,
			 &lt;cfqueryparam cfsqltype=&quot;cf_sql_double&quot; value=&quot;#AetF#&quot; /&gt;,
			 &lt;cfqueryparam cfsqltype=&quot;cf_sql_double&quot; value=&quot;#AstF#&quot; /&gt;,
			 &lt;cfqueryparam cfsqltype=&quot;cf_sql_double&quot; value=&quot;#Det#&quot; /&gt;,
			 &lt;cfqueryparam cfsqltype=&quot;cf_sql_double&quot; value=&quot;#Dst#&quot; /&gt;,
			 &lt;cfqueryparam cfsqltype=&quot;cf_sql_double&quot; value=&quot;#Vet#&quot; /&gt;,
			 &lt;cfqueryparam cfsqltype=&quot;cf_sql_double&quot; value=&quot;#Vst#&quot; /&gt;,
			 &lt;cfif t EQ 0&gt;
			 	&lt;!--- we can't actually derive from distance these at t = 0, so we'll use the force-derived ones. ---&gt;  
				&lt;cfqueryparam cfsqltype=&quot;cf_sql_double&quot; value=&quot;#AetF#&quot; /&gt;,
				&lt;cfqueryparam cfsqltype=&quot;cf_sql_double&quot; value=&quot;#AstF#&quot; /&gt;
			 &lt;cfelse&gt;
				&lt;cfqueryparam cfsqltype=&quot;cf_sql_double&quot; value=&quot;#Aet#&quot; /&gt;,
				&lt;cfqueryparam cfsqltype=&quot;cf_sql_double&quot; value=&quot;#Ast#&quot; /&gt;
			 &lt;/cfif&gt;
			)
	&lt;/cfquery&gt;
&lt;/cffunction&gt;
&lt;cfset D_e_cache = {} /&gt;
&lt;cfset D_s_cache = {} /&gt;
&lt;cfif attributes.restart&gt;
	&lt;cfquery datasource=&quot;#variables.my.dsn#&quot; name=&quot;&quot;&gt;
		delete from point
	&lt;/cfquery&gt;
	&lt;cfset cache(0, 0, 0) /&gt;
	&lt;cfset lastT = 0 /&gt;
&lt;cfelse&gt;
	&lt;cfquery datasource=&quot;#variables.my.dsn#&quot; name=&quot;lastT&quot;&gt;
		select max(t) t
		from point
	&lt;/cfquery&gt;
	&lt;cfset lastT = lastT.t /&gt;&lt;!--- s ---&gt;
&lt;/cfif&gt;

&lt;!--- locked! ---&gt;
&lt;cfset attributes.t = 5486400 /&gt;&lt;!--- 63.5 days / ~88.8% ---&gt;
&lt;cfset attributes.showCharts = true /&gt;

&lt;cfset maxT = lastT + BUILD_INTERVAL /&gt;&lt;!--- s ---&gt;
&lt;cfif NOT structKeyExists(attributes, &quot;t&quot;)&gt;
	&lt;cfif structKeyExists(attributes, &quot;d&quot;) AND isNumeric(attributes.d) AND attributes.d GTE 0 AND attributes.d LTE BUILD_INTERVAL&gt;
		&lt;cfset attributes.t = lastT + attributes.d /&gt;
	&lt;cfelse&gt;
		&lt;cfset attributes.t = maxT /&gt;
	&lt;/cfif&gt;
&lt;cfelseif attributes.t GT maxT&gt;
	&lt;p&gt;Your supplied 't' value was greater than the maximum allowed (#numberFormat(maxT, ',')#), so it was lowered to that.
	&lt;/p&gt;
	&lt;cfset attributes.t = maxT /&gt;
&lt;/cfif&gt;
&lt;cfquery datasource=&quot;#variables.my.dsn#&quot; name=&quot;get&quot;&gt;
	(
		select t, Det, Dst
		from point
		where t &lt;= &lt;cfqueryparam cfsqltype=&quot;cf_sql_integer&quot; value=&quot;#attributes.t#&quot; /&gt;
		order by t desc
		limit 10
	) union distinct (
		select t, Det, Dst
		from point
		order by t desc
		limit 10
	)
&lt;/cfquery&gt;
&lt;cfloop query=&quot;get&quot;&gt;
	&lt;cfset D_e_cache[t] = Det /&gt;
	&lt;cfset D_s_cache[t] = Dst /&gt;
&lt;/cfloop&gt;

&lt;cfscript&gt;
function F(m1, m2, r) {
	return G * m1 * m2 / r ^ 2; // kg * m  / s^2
}
function a(F, m) {
	return F / m; // m / s^2
}
D_e_localCache = {
	0 = 0
};
D_s_localCache = duplicate(D_e_localCache);
structAppend(D_e_localCache, D_e_cache);
structAppend(D_s_localCache, D_s_cache);
function V_e(t) {
	return V_x(D_e, t);
}
function V_s(t) {
	return V_x(D_s, t);
}
function V_x(d, t) {
	if (t LT 1) {
		return 0;
	}
	return d(t) - d(t - 1); // m / s
}
function A_e(t) {
	return A_x(V_e, t);
}
function A_s(t) {
	return A_x(V_s, t);
}
function A_x(v, t) {
	if (t LT 2) {
		return 0;
	}
	return (v(t) - v(t - 1)); // m / s^2
}
function D_e(t) {
	return D_x(t, Me, D_e_cache, D_e_localCache);
}
function D_s(t) {
	return D_x(t, Ms, D_s_cache, D_s_localCache);
}
function D_x(t, M, cache, localCache) {
	if (t LTE 0) {
		return 0;
	} else if (structKeyExists(localCache, t)) {
		return localCache[t];
	} else {
		build_D_x(t);
		return localCache[t];
	}
}
function maxLte(s, v) {
	var result = 0;
	for (i in s) {
		if (i GT result AND i LTE v) {
			result = i;
		}
	}
	return result;
}
function build_D_x(t) {
	var re = 0;
	var rs = 0;
	var n = 1 + min(
		maxLte(D_e_localCache, t),
		maxLte(D_s_localCache, t)
	);
	var lF = 0;
	var De1 = 0;
	var Ds1 = 0;
	writeOutput(&quot;&lt;pre&gt;build_D_x(#t#) w/ n = #n#&lt;/pre&gt;&quot;);
	for (; n LTE t; n += 1) {
		De1 = D_e(n - 1);
		Ds1 = D_s(n - 1);
		lF = F(Me, Ms, r - De1 - Ds1);
		re = De1 + V_e(n - 1) + a(lF, Me); // m
		rs = Ds1 + V_s(n - 1) + a(lF, Ms) ; // m
		D_e_localCache[n] = re;
		D_s_localCache[n] = rs;
		if (n MOD CHART_INTERVAL EQ 0) {
			// cache the magic intervals
			cache(n - 1, D_e_localCache[n - 1], D_s_localCache[n - 1]);
			cache(n, re, rs);
		}
	}
	cache(t - 1, D_e_localCache[t - 1], D_s_localCache[t - 1]);
	cache(t, D_e_localCache[t], D_s_localCache[t]);
}
function log2(n) {
	// log2(x) = log2(10) * log10(x) = log10(x) / log10(2)
	return log10(n) / log10(2);
}
function n(n) {
	if (n LT 10000 AND n GTE 1000) {
		return numberFormat(n, &quot;9999&quot;);
	} else if (n LT 1000 AND n GTE 100) {
		return numberFormat(n, &quot;999.0&quot;);
	} else if (n LT 100 AND n GTE 10) {
		return numberFormat(n, &quot;99.00&quot;);
	} else if (n LT 10 AND n GTE 1) {
		return numberFormat(n, &quot;0.000&quot;);
	} else if (n LT 1 AND n GTE 0.1) {
		return numberFormat(n, &quot;0.0000&quot;);
	} else if (n LT 0.1 AND n GTE 0.01) {
		return numberFormat(n, &quot;0.00000&quot;);
	} else if (n LT 0.01 AND n GTE 0.001) {
		return numberFormat(n, &quot;0.000000&quot;);
	}
	return REReplace(scientificFormat(n, 4), 'E(-?[0-9]+)', ' x 10&lt;sup&gt;\1&lt;/sup&gt;');
}
function scientificFormat(n,m) {
    return createObject(&quot;java&quot;,&quot;java.text.DecimalFormat&quot;).init(&quot;0.#repeatString(&quot;0&quot;, m - 1)#E0&quot;).format(createObject(&quot;java&quot;,&quot;java.math.BigDecimal&quot;).init(n));
}
function printLookup(c) {
	var a = listToArray(structKeyList(c));
	var i = 0;
	var l = arrayLen(a);
	var ml = 0;
	arraySort(a, 'numeric');
	ml = len(a[l]);
	for (i = 1; i LTE l; i++) {
		writeOutput(&quot;#rJustify(a[i], ml)#=#c[a[i]]##chr(10)#&quot;);
	}
}
&lt;/cfscript&gt;

&lt;p&gt;What if the earth suddenly stopped orbiting the sun?  The gravitational pull would
suddenly no longer have the earth's inertia to counteract it, and the earth would
plummet into the sun.  But how fast?  How long would it take?
&lt;/p&gt;

&lt;p&gt;What if both bodies were point masses (i.e., had no volume), there was no
other motion, and were no other bodies around?  Now THAT, I can approximate.  :)
&lt;/p&gt;

&lt;p&gt;These are the results of that model, evaluated each second, starting at &lt;code style=&quot;white-space:nowrap;&quot;&gt;t = 0&lt;/code&gt;: 
&lt;/p&gt;
&lt;div class=&quot;row&quot;&gt;
&lt;pre class=&quot;third&quot;&gt;
Names:
------
t = time
r = Earth-Sun separation (orbital radius)
F = gravitational force
a = acceleration due to gravity
D = distance covered
V = velocity (based on &amp;##916;distance)
A = acceleration (based on &amp;##916;velocity)
&lt;/pre&gt;
&lt;pre class=&quot;third&quot;&gt;
Formulae:
---------

F(M&lt;sub&gt;1&lt;/sub&gt;, M&lt;sub&gt;2&lt;/sub&gt;, r) = G * M&lt;sub&gt;1&lt;/sub&gt; * M&lt;sub&gt;2&lt;/sub&gt; / r&lt;sup&gt;2&lt;/sup&gt;
a(M, r) = F(M, M&lt;sub&gt;other&lt;/sub&gt;, r) / M
D(t) = D(t - 1) + V(t - 1) + A(t - 1)
V(t) = D(t) - D(t - 1)
A(t) = V(t) - V(t - 1)
&lt;/pre&gt;
&lt;cfset F0 = F(Me, ms, r) /&gt;
&lt;pre class=&quot;third&quot;&gt;
Constants:
----------
M&lt;sub&gt;e&lt;/sub&gt; = Mass&lt;sub&gt;earth&lt;/sub&gt; = #n(Me)# kg
M&lt;sub&gt;s&lt;/sub&gt; = Mass&lt;sub&gt;sun&lt;/sub&gt; = #n(Ms)# kg
r(0) = Earth's mean orbital radius = #n(r)# m
G = Gravitational constant = #n(G)# m&lt;sup&gt;3&lt;/sup&gt; / (kg * s&lt;sup&gt;2&lt;/sup&gt;)
&lt;/pre&gt;
&lt;/div&gt;

&lt;cfloop list=&quot;0,#attributes.t#,#lastT#&quot; index=&quot;t&quot;&gt;
	&lt;cfset Det = D_e(t) /&gt;
	&lt;cfset Dst = D_s(t) /&gt;
	&lt;cfset Vet = V_e(t) /&gt;
	&lt;cfset Vst = V_s(t) /&gt;
	&lt;cfset Aet = A_e(t) /&gt;
	&lt;cfset Ast = A_s(t) /&gt;
	&lt;cfset Dt = Dst + Det /&gt;
	&lt;cfset rt = r - Dt /&gt;
	&lt;cfset Ft = F(Me, Ms, rt) /&gt;
	&lt;cfset AetF = a(Ft, Me) /&gt;
	&lt;cfset AstF = a(Ft, Ms) /&gt;
	&lt;div class=&quot;row&quot;&gt;
	&lt;cfset s = 't = #numberFormat(t, &quot;,&quot;)# s (#round(t / 86400 * 100) / 100# days)' /&gt;
	&lt;pre class=&quot;third&quot;&gt;
#s##chr(10)##repeatString(&quot;-&quot;, len(s))#
D(t) = #n(Dt)# m (#n(Dt / r * 100)#%)
r(t) = #n(rt)# m&lt;cfif Vet + Vst Gt 0&gt; (&amp;lt; #decimalTimespanFormat(rt / (Vet + Vst))#)&lt;/cfif&gt;
F(t) = #n(Ft)# kg * m / s&lt;sup&gt;2&lt;/sup&gt;
	&lt;/pre&gt;
	&lt;pre class=&quot;third&quot;&gt;
a&lt;sub&gt;e&lt;/sub&gt;(t) = #n(AetF)# m / s&lt;sup&gt;2&lt;/sup&gt;
D&lt;sub&gt;e&lt;/sub&gt;(t) = #n(Det)# m
V&lt;sub&gt;e&lt;/sub&gt;(t) = #n(Vet)# m / s
A&lt;sub&gt;e&lt;/sub&gt;(t) = #n(Aet)# m / s&lt;sup&gt;2&lt;/sup&gt;
	&lt;/pre&gt;
	&lt;pre class=&quot;third&quot;&gt;
a&lt;sub&gt;s&lt;/sub&gt;(t) = #n(AstF)# m / s&lt;sup&gt;2&lt;/sup&gt;
D&lt;sub&gt;s&lt;/sub&gt;(t) = #n(Dst)# m
V&lt;sub&gt;s&lt;/sub&gt;(t) = #n(Vst)# m / s
A&lt;sub&gt;s&lt;/sub&gt;(t) = #n(Ast)# m / s&lt;sup&gt;2&lt;/sup&gt;
	&lt;/pre&gt;
	&lt;/div&gt;
&lt;/cfloop&gt;
&lt;pre&gt;#numberFormat(getTickCount() - startTime, &quot;,&quot;)# ms&lt;/pre&gt;
&lt;cfif attributes.iterate&gt;
	&lt;cfif rt GTE 0.1 * r&gt;
		&lt;meta http-equiv=&quot;refresh&quot; content=&quot;1&quot; /&gt;
	&lt;cfelse&gt;
		&lt;cfset attributes.showCharts = true /&gt;
	&lt;/cfif&gt;
&lt;/cfif&gt;
&lt;cfset startTime = getTickCount() /&gt;

&lt;cfif NOT structKeyExists(attributes, &quot;showCharts&quot;) OR NOT isBoolean(attributes.showCharts) OR NOT attributes.showCharts&gt;
	&lt;p&gt;Add '&lt;code&gt;&lt;a href=&quot;?showCharts=true&quot;&gt;?showCharts=true&lt;/a&gt;&lt;/code&gt;' to see some charts of the progress
	&lt;/p&gt;
	&lt;cfexit method=&quot;exittemplate&quot; /&gt;
&lt;/cfif&gt;
&lt;div id=&quot;loading&quot;&gt;Loading charts...&lt;/div&gt;
&lt;cfflush /&gt;
&lt;script type=&quot;text/javascript&quot;&gt;try { var el = document.getElementById(&quot;loading&quot;);el.parentNode.removeChild(el); } catch (e) { /* oh well */ }&lt;/script&gt;

&lt;!--- lets build some charts! ---&gt;
&lt;style type=&quot;text/css&quot;&gt;
.charts img {
	border: 1px solid ##ddd;
	margin: 5px;
	padding: 5px;
}
&lt;/style&gt;
&lt;cfset width = 600 /&gt;
&lt;cfset height = 200 /&gt;
&lt;cfset encoding = &quot;e&quot; /&gt;
&lt;cfset Ce = &quot;0099ff&quot; /&gt;
&lt;cfset Cs = &quot;ffee00&quot; /&gt;
&lt;cfset CF = &quot;cc6600&quot; /&gt;
&lt;cfquery datasource=&quot;#variables.my.dsn#&quot; name=&quot;get&quot;&gt;
	select max(t) maxT
	from point
	where t &lt;= &lt;cfqueryparam cfsqltype=&quot;cf_sql_integer&quot; value=&quot;#attributes.t#&quot; /&gt;
&lt;/cfquery&gt;
&lt;cfif get.recordCount EQ 0 OR NOT isNumeric(get.maxT)&gt;
	&lt;p&gt;There isn't enough data to render charts yet...
	&lt;/p&gt;
	&lt;cfexit method=&quot;exittemplate&quot; /&gt;
&lt;/cfif&gt;
&lt;p&gt;Charts are plotted through #numberFormat(attributes.t, &quot;,&quot;)# sec (#attributes.t / 86400# days).  Gravity goes insane after that.  Be careful of y-axis units; they're inconsistent (for different scales), and some are unusual (e.g., gigameters).
&lt;/p&gt;
&lt;cfset interval = ceiling(get.maxT / 40 / CHART_INTERVAL) * CHART_INTERVAL /&gt;
&lt;cfquery datasource=&quot;#variables.my.dsn#&quot; name=&quot;get&quot;&gt;
	select t as tt, t / 86400 as t, -- hours
		Dt / 1000 / 1000 / 1000 as Dt, -- gigameters
		rt / 1000 / 1000 / 1000 as rt, -- gigameters
		(#r# - Det) / 1000 / 1000 / 1000 as Pet, -- gigameters
		Dst / 1000 / 1000 / 1000 as Pst, -- gigameters
		Ft / 1000 / 1000 / 1000 / 1000 / 1000 / 1000 / 1000 Ft, -- zettanewtons
		Det / 1000 / 1000 / 1000 as Det, -- gigameters
		Dst / 1000 as Dst, -- kilometers
		Vet / 1000 as Vet, -- kilometers / second
		Vst * 100 as Vst, -- centimeters / second
		Aet as Aet, -- meters / second
		Ast * 1000 * 1000 as Ast -- micrometers / second
	from point
	where (t &lt; &lt;cfqueryparam cfsqltype=&quot;cf_sql_integer&quot; value=&quot;#get.maxT#&quot; /&gt;
			and t MOD &lt;cfqueryparam cfsqltype=&quot;cf_sql_integer&quot; value=&quot;#interval#&quot; /&gt; = 0
		)
		or t = &lt;cfqueryparam cfsqltype=&quot;cf_sql_integer&quot; value=&quot;#get.maxT#&quot; /&gt;
	order by t
&lt;/cfquery&gt;
&lt;cffunction name=&quot;chart&quot;&gt;
	&lt;cfargument name=&quot;labels&quot; /&gt;
	&lt;cfargument name=&quot;series&quot; /&gt;
	&lt;cfargument name=&quot;sharedAxis&quot; default=&quot;false&quot; /&gt;
	&lt;cfargument name=&quot;colors&quot; default=&quot;#Ce#,#Cs#&quot; /&gt;
	&lt;cfset var i = 0 /&gt;
	&lt;cfset var temp = &quot;&quot; /&gt;
	&lt;cfif isSimpleValue(labels)&gt;
		&lt;cfset labels = listToArray(labels, &quot;|&quot;) /&gt;
	&lt;/cfif&gt;
	&lt;cfif isSimpleValue(series)&gt;
		&lt;cfset series = listToArray(series, &quot;,&quot;) /&gt;
	&lt;/cfif&gt;
	&lt;cfif isSimpleValue(colors)&gt;
		&lt;cfset colors = listToArray(colors) /&gt;
	&lt;/cfif&gt;
	&lt;cfloop from=&quot;1&quot; to=&quot;#arrayLen(series)#&quot; index=&quot;i&quot;&gt;
		&lt;cfset temp = series[i] /&gt;
		&lt;cfset series[i] = {
			name = temp,
			label = labels[i],
			color = colors[(i - 1) MOD arrayLen(colors) + 1]
		} /&gt;
	&lt;/cfloop&gt;
	&lt;cfif sharedAxis&gt;
		&lt;cfset temp = {
			min = 0,
			max = 0
		} /&gt;
		&lt;cfloop from=&quot;1&quot; to=&quot;#arrayLen(series)#&quot; index=&quot;i&quot;&gt;
			&lt;cfset temp.min = min(temp.min, arrayMin(valueArray(&quot;get.#series[i].name#&quot;))) /&gt;
			&lt;cfset temp.max = max(temp.max, arrayMax(valueArray(&quot;get.#series[i].name#&quot;))) /&gt;
		&lt;/cfloop&gt;
		&lt;cfset temp = getAxisBounds(temp.min, temp.max) /&gt;
		&lt;cfloop from=&quot;1&quot; to=&quot;#arrayLen(series)#&quot; index=&quot;i&quot;&gt;
			&lt;cfset series[i].bounds = temp /&gt;
		&lt;/cfloop&gt;
	&lt;cfelse&gt;
		&lt;cfloop from=&quot;1&quot; to=&quot;#arrayLen(series)#&quot; index=&quot;i&quot;&gt;
			&lt;cfset series[i].bounds = getAxisBounds(arrayMin(valueArray(&quot;get.#series[i].name#&quot;)), arrayMax(valueArray(&quot;get.#series[i].name#&quot;))) /&gt;
		&lt;/cfloop&gt;
	&lt;/cfif&gt;
	&lt;cfsavecontent variable=&quot;chartUrl&quot;&gt;
		http://chart.apis.google.com/chart
		?cht=lc
		&amp;chs=#width#x#height#
		&amp;chd=#encoding#:
			&lt;cfloop from=&quot;1&quot; to=&quot;#arrayLen(series)#&quot; index=&quot;i&quot;&gt;
				&lt;cfif i GT 1&gt;,&lt;/cfif&gt;
				#chartEncode(valueList(&quot;get.#series[i].name#&quot;), series[i].bounds.lower, series[i].bounds.upper, encoding)#
			&lt;/cfloop&gt;
		&amp;chco=&lt;cfloop from=&quot;1&quot; to=&quot;#arrayLen(series)#&quot; index=&quot;i&quot;&gt;&lt;cfif i GT 1&gt;,&lt;/cfif&gt;#series[i].color#&lt;/cfloop&gt;
		&amp;chm=&lt;cfloop from=&quot;1&quot; to=&quot;#arrayLen(series)#&quot; index=&quot;i&quot;&gt;&lt;cfif i GT 1&gt;|&lt;/cfif&gt;d,#series[i].color#,#i - 1#,-1,7&lt;/cfloop&gt;
		&amp;chdl=&lt;cfloop from=&quot;1&quot; to=&quot;#arrayLen(series)#&quot; index=&quot;i&quot;&gt;&lt;cfif i GT 1&gt;|&lt;/cfif&gt;#urlEncodedFormat(series[i].label)#&lt;/cfloop&gt;
		&amp;chdlp=t
		&amp;chxt=x,x&lt;cfif sharedAxis&gt;,r&lt;cfelse&gt;#repeatString(&quot;,r&quot;, arrayLen(series))#&lt;/cfif&gt;
		&amp;chxl=0:&lt;cfloop query=&quot;get&quot;&gt;|&lt;cfif currentRow MOD 4 EQ 1&gt;#numberFormat(t, &quot;.0&quot;)#&lt;/cfif&gt;&lt;/cfloop&gt;
			|1:|time+(days)
			&lt;cfloop from=&quot;1&quot; to=&quot;#arrayLen(series)#&quot; index=&quot;i&quot;&gt;
				|#1 + i#:&lt;cfloop from=&quot;#series[i].bounds.lower#&quot; to=&quot;#series[i].bounds.upper#&quot; index=&quot;i&quot; step=&quot;#series[i].bounds.interval * iif(series[i].bounds.intervalCount GT 7, 2, 1)#&quot;&gt;|#i#&lt;/cfloop&gt;
				&lt;cfif sharedAxis&gt;
					&lt;cfbreak /&gt;
				&lt;/cfif&gt;
			&lt;/cfloop&gt;
		&amp;chxp=1,50
		&lt;cfif arrayLen(series) GT 1 AND NOT sharedAxis&gt;
			&lt;!--- colored axis labels ---&gt;
			&amp;chxs=
				&lt;cfloop from=&quot;1&quot; to=&quot;#arrayLen(series)#&quot; index=&quot;i&quot;&gt;
					&lt;cfif i GT 1&gt;|&lt;/cfif&gt;
					#1 + i#,#series[i].color#
				&lt;/cfloop&gt;
		&lt;cfelse&gt;
			&amp;chg=20,#int(100 / (series[1].bounds.intervalCount / iif(series[1].bounds.intervalCount GT 7, 2, 1)) * 100) / 100#
		&lt;/cfif&gt;
	&lt;/cfsavecontent&gt;
	&lt;cfset chartUrl = REReplace(chartUrl, &quot;[[:space:]]+&quot;, &quot;&quot;, &quot;all&quot;) /&gt;
	&lt;!---
	&lt;pre&gt;#len(chartUrl)# chars#chr(10)##replace(chartUrl, &quot;&amp;&quot;, &quot;#chr(10)#  &amp;&quot;, &quot;all&quot;)#&lt;/pre&gt;
	---&gt;
	&lt;img src=&quot;#chartUrl#&quot; /&gt;
&lt;/cffunction&gt;
&lt;div class=&quot;charts&quot;&gt;
&lt;cfset chart(&quot;De(t) (Gm)|Ds(t) (km)&quot;, &quot;Det,Dst&quot;, true) /&gt;
&lt;cfset chart(&quot;Pe(t) (Gm)|Ps(t) (Gm)&quot;, &quot;Pet,Pst&quot;, true) /&gt;
&lt;cfset chart(&quot;Ve(t) (km / s)|Vs(t) (cm / s)&quot;, &quot;Vet,Vst&quot;, true) /&gt;
&lt;cfset chart(&quot;F(t) (ZN)&quot;, &quot;Ft&quot;, false, CF) /&gt;
&lt;cfset chart(&quot;Ae(t) (m / s^2)|As(t) (μm / s^2)&quot;, &quot;Aet,Ast&quot;, true) /&gt;
&lt;/div&gt;

&lt;pre&gt;#numberFormat(getTickCount() - startTime, &quot;,&quot;)# ms&lt;/pre&gt;

&lt;a href=&quot;sun_earth_code.cfm&quot;&gt;View the code (masochists only)&lt;/a&gt;
&lt;/body&gt;
&lt;/html&gt;
&lt;/cfoutput&gt;</pre>