JavaScript Gravity Simulator: Check One Mutual Attraction Step

HTML for this example

Select the code and use your copy shortcut. Changes in this preview are temporary.

<!doctype html>
<html lang="en-US"><head><meta charset="utf-8"><meta name="viewport" content="width=device-width,initial-scale=1"><title>JavaScript Gravity Simulator: Check One Mutual Attraction Step</title><style>:root{--bg:#eaf3f7;--paper:#fff;--ink:#16354a;--accent:#12526f;--tint:#dcedf5;--line:#b5ccd7;--muted:#436171;--button:#fff}*{box-sizing:border-box}body{margin:0;background:var(--bg);color:var(--ink);font:17px/1.55 system-ui,sans-serif}main{max-width:1000px;margin:auto;padding:22px clamp(18px,4vw,42px) 38px}header{border-bottom:2px solid var(--accent);padding-bottom:14px;margin-bottom:18px}h1{font-size:clamp(26px,4vw,36px);line-height:1.15;margin:4px 0 10px}h2{font-size:22px;line-height:1.3;margin:0 0 14px}h3{font-size:19px}p{margin:0 0 14px}.eyebrow{font-size:16px;font-weight:750;color:var(--accent)}.lede{max-width:800px}.panel{background:var(--paper);border:1px solid var(--line);border-radius:9px;padding:20px;margin:18px 0}.workbench{border-top:4px solid var(--accent)}.grid{display:grid;grid-template-columns:repeat(2,minmax(0,1fr));gap:16px}.grid>*{min-width:0}label{display:block;font-size:16px;font-weight:650;margin:12px 0 6px}input,textarea,select{font:inherit;color:var(--ink);background:var(--paper);border:1px solid #7b8893;border-radius:5px;padding:10px;width:100%;min-width:0}textarea{resize:vertical;line-height:1.45}input[type=checkbox],input[type=radio]{width:21px;height:21px;vertical-align:middle;margin-right:9px}button{font:650 16px/1.4 system-ui,sans-serif;padding:11px 15px;min-height:46px;background:var(--accent);border:2px solid var(--accent);color:var(--button);border-radius:5px;cursor:pointer}button:disabled{opacity:.55;cursor:default}.buttons{display:flex;gap:9px;flex-wrap:wrap;margin:14px 0}.buttons>*{flex:1 1 auto}button:focus-visible,input:focus-visible,textarea:focus-visible,select:focus-visible,summary:focus-visible{outline:3px solid #c05b00;outline-offset:3px}.result{font:600 17px/1.6 ui-monospace,monospace;white-space:pre-wrap;overflow-wrap:anywhere;background:var(--tint);padding:16px;border-left:4px solid var(--accent);min-height:65px;margin:16px 0}.sample{padding:14px 16px;border-left:3px solid var(--accent);background:var(--tint);margin:12px 0}.muted{color:var(--muted)}pre,code{font:16px/1.6 ui-monospace,monospace;white-space:pre-wrap;overflow-wrap:anywhere}pre{margin:12px 0}table{border-collapse:collapse;width:100%;table-layout:fixed;font-size:16px}th,td{text-align:left;vertical-align:top;padding:10px 8px;border-bottom:1px solid var(--line);overflow-wrap:anywhere}th{background:var(--tint)}svg{max-width:100%;height:auto;display:block}details{margin:20px 0}summary{font-weight:700;cursor:pointer}footer{margin-top:28px;border-top:1px solid var(--line);padding-top:18px;font-size:16px}a{color:var(--accent)}.reading{font-family:Georgia,serif}.codeview{border-radius:0;box-shadow:7px 7px 0 var(--tint)}.token{display:inline-block;padding:8px 12px;margin:4px;border:1px solid var(--line);border-radius:4px;background:var(--paper)}.swatch{min-height:30px}.split{display:grid;grid-template-columns:1.1fr 1fr;gap:22px}ol{padding-left:25px}li{margin:9px 0}@media(max-width:600px){main{padding:16px}.grid,.split{grid-template-columns:1fr}.panel{padding:16px;margin:14px 0}h1{font-size:27px}header{padding-bottom:10px}.buttons{gap:8px}.result{padding:13px}th,td{padding:8px 5px}}</style></head><body><main><header><div class="eyebrow">Original practice fixture · local workbench</div><h1>Equal force, unequal acceleration</h1><p class="lede">Compare equal and opposite forces with unequal accelerations in one step. Settings: G = 1; h = 0.1; masses A = 1 and B = 2; positions A = −1 and B = 1; initial velocities = 0. Expected force magnitudes: 0.5; accelerations: 0.5 and −0.25; velocities: 0.05 and −0.025; positions: −0.995 and 0.9975; momentum: 0; center of mass: 1/3.</p></header><pre id="result" class="result" role="status" aria-live="polite" aria-atomic="true">Ready: A at −1 mass 1; B at 1, mass 2; velocities 0.</pre><div class="buttons"><button type="button" id="step">Calculate one step</button><button type="button" id="reset">Reset</button></div><label for="ma">Mass A</label><input id="ma" value="1" ><label for="mb">Mass B</label><input id="mb" value="2" ><label for="xa">Initial x A</label><input id="xa" value="-1" ><label for="xb">Initial x B</label><input id="xb" value="1" ><label for="g">Model G</label><input id="g" value="1" ><label for="h">Step h</label><input id="h" value=".1" ><footer><p>The dimensionless semi-implicit Euler exercise is not an orbital predictor, collision solver, or physical validation.</p><p>Use a prepared HTML explanation in <a href="https://nos-workspace.com/">NOS</a> when reviewing this mutual-attraction step calculator. Eligible static instructions can be revised in visual editing mode; dynamic outputs follow their controls.</p></footer></main><script>
(()=>{'use strict';const $=s=>document.querySelector(s);const out=s=>{$('#result').textContent=s;};const val=id=>$('#'+id).value;const num=id=>{const s=val(id).trim();return s===''?NaN:Number(s);};const on=(id,fn)=>$('#'+id).addEventListener('click',fn);const watch=ids=>{for(const id of ids)for(const event of ['input','change'])$('#'+id).addEventListener(event,()=>out('Inputs changed. Run the action to update this result.'));};const fmt=(n,d=4)=>Math.abs(n)<1e-10?(0).toFixed(d):n.toFixed(d);
let done=false;const ids=['ma','mb','xa','xb','g','h'];for(const id of ids)$('#'+id).addEventListener('input',()=>{done=false;out('Inputs changed. Calculate a fresh step.');});on('step',()=>{if(done){out('Single step already calculated. Change inputs or Reset.');return;}const [ma,mb,xa,xb,g,h]=ids.map(num);if(![ma,mb,xa,xb,g,h].every(Number.isFinite)||Math.min(ma,mb,g,h)<=0){out('Invalid: finite positions and positive masses, G, h required.');return;}if(xa===xb){out('Held: coincident centers.');return;}const d=xb-xa,F=g*ma*mb/(d*d)*Math.sign(d),aa=F/ma,ab=-F/mb,va=h*aa,vb=h*ab,na=xa+h*va,nb=xb+h*vb,momentum=ma*va+mb*vb,center=(ma*na+mb*nb)/(ma+mb);if(![F,aa,ab,va,vb,na,nb,momentum,center].every(Number.isFinite)){out('Held: nonfinite computed state.');return;}done=true;out('Force '+fmt(F)+' / '+fmt(-F)+' | Acceleration '+fmt(aa)+' / '+fmt(ab)+' | Velocity '+fmt(va)+' / '+fmt(vb)+' | Position '+fmt(na)+' / '+fmt(nb)+' | Momentum '+fmt(momentum)+' | Center '+fmt(center,6));});on('reset',()=>{ids.forEach((id,i)=>$('#'+id).value=['1','2','-1','1','1','.1'][i]);done=false;out('Ready: both velocities zero.');});
})();
</script></body></html>