Simulate the Solar System FOR REAL

{"name":"Simulate the Solar System","description":"Experience the wonders of gravity in this interactive 2D simulation! Create your own universe and watch as celestial bodies attract each other based on mass and distance.","photo":"https://res.cloudinary.com/dzv7hjwfn/image/upload/lzam3qxjwjkomxvifq5r","view_count":26,"fork_count":0,"expires_at":null,"created_at":"2025-08-06T08:51:39.723122Z","updated_at":"2025-08-06T08:51:39.723123Z","user_id":246,"status":"DELETED","parent_id":2106,"id":2108,"user":{"photo":null,"email":"kristen@recess.gg","recess_id":"c4ad84f6-113d-480b-bc29-80549518b3ee","created_at":"2024-10-31T20:27:47.342142Z","reset_token":null,"stripe_customer_id":null,"seats":null,"username":"chicothedog","id":246,"password_hash":"$2b$12$QgZsBBZ/7CFyOhVkfVhTNeoe4jAY.pIJOieTXeg0HWS23ZcCpiAeO","updated_at":"2024-10-31T20:27:47.342148Z","reset_token_expiry":null,"subscription_status":null},"code":{"script.js":"/*\\n==============================================\\nSIMULATE THE SOLAR SYSTEM: the ~fancy~ version\\n==============================================\\n\\nFirst off, credit where credit is due. Took a lot from:\\nhttp://www.planetaryorbits.com/tutorial-javascript-orbit-simulation.html\\n\\nBelow you will find my adaptation of code from the above link.\\n\\nI want to admit my mistakes, show how I corrected them, and give you guys a chance to learn from them.\\n\\nMy main mistake was thinking that we could build an accurate simulation with minimal celestial data.\\n\\nAll the physics we've done so far, from velocity and acceleration to force and mass to trigonometry to momentum, is 100% legit.\\n\\nIn principle, we could dial in exact positions and velocity values and create stable orbits.\\n\\nBut the problem is with the precision and accuracy necessary to accomplish this.\\nThe solar system is a big place, and being a pixel off on our little screens is like being thousands of miles off in space.\\nOn top of that, we have to convert between physical units — kilometers and days — to computer units — pixels and frames.\\nAny tiny error there (even due to rounding) can shatter the delicate gravitational balance of the solar system.\\n\\nThe bottom line: \\nOur physics so far will work for many purposes, including games, art, and simulations, all with pretty light code.\\nBut to achieve our true goal of Simulating the Actual real-deal Solar System, we need a little more data and a little more precision.\\n\\nI tracked down some helpful data and some code that accomplishes this, and did what I could to make it readable and editable.\\n\\nI also added the asteroids (which are not realistic, and are really more like comets than asteroids, but are fun nonetheless...).\\n\\nIn any case, behold! Admire the rhythm of the cosmos.\\n\\nAnd play around! Add a soundtrack. Crank up gravity. Program a controllable spaceship that can bump into the planets.\\nThe only limit is your imagination.\\n*/\\n\\n\\n// First up, the basics. Variables to hold the images of the planets.\\nlet sunImg, mercuryImg, venusImg, earthImg, marsImg, jupiterImg, saturnImg, uranusImg, neptuneImg;\\n\\n// Load 'em up.\\nfunction preload() {\\n sunImg = loadImage(\"@asset('sun')\");\\n mercuryImg = loadImage(\"@asset('mercury')\");\\n venusImg = loadImage(\"@asset('venus')\");\\n earthImg = loadImage(\"@asset('earth')\");\\n marsImg = loadImage(\"@asset('mars')\");\\n jupiterImg = loadImage(\"@asset('jupiter')\");\\n saturnImg = loadImage(\"@asset('saturn')\");\\n uranusImg = loadImage(\"@asset('uranus')\");\\n neptuneImg = loadImage(\"@asset('neptune')\");\\n}\\n\\n// An empty list that will hold asteroids.\\nlet asteroids = [];\\n\\n// For simplicity, I am setting all asteroids to the same mass.\\n// You can change that!\\nconst asteroidMass = 3.0e14;\\n\\n// A basic setup function. \"T6\" is the name of the main simulation.\\nfunction setup() {\\n createCanvas(windowWidth, windowHeight);\\n frameRate(60);\\n T6.contextbackground = createGraphics(windowWidth, windowHeight);\\n T6.contextforeground = createGraphics(windowWidth, windowHeight);\\n init_T6();\\n}\\n\\n// Click to make an asteroid.\\nfunction mousePressed() {\\n addAsteroid();\\n}\\n\\n// Most everything happens in T6. That means draw() is small and simple — nice!\\nfunction draw() {\\n run_T6();\\n image(T6.contextbackground, 0, 0);\\n image(T6.contextforeground, 0, 0);\\n drawAsteroids();\\n}\\n\\n/*\\n OK now for the big cheese. T6 is gonna hold a big list of properties.\\n Most will start as 'null' and get filled in during the actual simulation.\\n Then we will write some functions that work with T6 (which stores data) to make stuff *move* — that's physics!\\n*/\\nvar T6 = {\\n contextbackground: null,\\n contextforeground: null,\\n width: window.innerWidth,\\n height: window.innerHeight,\\n fps: 60,\\n now: null,\\n then: Date.now(),\\n interval: null,\\n delta: null,\\n nowTPS: null,\\n thenTPS: Date.now(),\\n avgTPSCount: 0,\\n TPSCount: 0,\\n deltaTPS: 0,\\n nowFPS: null,\\n thenFPS: Date.now(),\\n avgFPSCount: 0,\\n FPSCount: 0,\\n deltaFPS: 0,\\n julianCenturyInJulianDays: 36525,\\n julianEpochJ2000: 2451545.0,\\n julianDate: null,\\n current: null,\\n newDate: null,\\n DAY: null,\\n MONTH: null,\\n YEAR: null,\\n // This looks crazy, I know. These are numbers from NASA. Google unfamiliar words if you want to learn more.\\n // From the original data:\\n // At January 2000: semimajor Axis, Eccentricity, Inclination, mean longitude (L), longitude of perihelion, longitude of ascending node\\n planetElements: [\\n [0.38709927, 0.20563593, 7.00497902, 252.25032350, 77.45779628, 48.33076593],\\n [0.72333566, 0.00677672, 3.39467605, 181.97909950, 131.60246718, 76.67984255],\\n [1.00000261, 0.01671123, -0.00001531, 100.46457166, 102.93768193, 0.0],\\n [1.52371034, 0.09339410, 1.84969142, -4.55343205, -23.94362959, 49.55953891],\\n [5.20288700, 0.04838624, 1.30439695, 34.39644051, 14.72847983, 100.47390909],\\n [9.53667594, 0.05386179, 2.48599187, 49.95424423, 92.59887831, 113.66242448],\\n [19.18916464, 0.04725744, 0.77263783, 313.23810451, 170.95427630, 74.01692503],\\n [30.06992276, 0.00859048, 1.77004347, -55.12002969, 44.96476227, 131.78422574]\\n ],\\n\\n // rates: A, E, I, mean longitude (L), longitude of perihelion, longitude of ascending node\\n planetRates: [\\n [0.00000037, 0.00001906, -0.00594749, 149472.67411175, 0.16047689, -0.1253408],\\n [0.00000390, -0.00004107, -0.00078890, 58517.81538729, 0.00268329, -0.27769418],\\n [0.00000562, -0.00004392, -0.01294668, 35999.37244981, 0.32327364, 0.0],\\n [0.00001847, 0.00007882, -0.00813131, 19140.30268499, 0.44441088, -0.29257343],\\n [-0.00011607, -0.00013253, -0.00183714, 3034.74612775, 0.21252668, 0.20469106],\\n [-0.00125060, -0.00050991, 0.00193609, 1222.49362201, -0.41897216, -0.28867794],\\n [-0.00196176, -0.00004397, -0.00242939, 428.48202785, 0.40805281, 0.04240589],\\n [0.00026291, 0.00005105, 0.00035372, 218.45945325, -0.32241464, -0.00508664]\\n ],\\n orbitalElements: null,\\n xMercury: null,\\n yMercury: null,\\n xVenus: null,\\n yVenus: null,\\n xEarth: null,\\n yEarth: null,\\n xMars: null,\\n yMars: null,\\n xJupiter: null,\\n yJupiter: null,\\n xSaturn: null,\\n ySaturn: null,\\n xUranus: null,\\n yUranus: null,\\n xNeptune: null,\\n yNeptune: null,\\n // We scale planets' images a bit so we can see them better. Otherwise, big is too big or small is too small.\\n scale: 150,\\n jupiterScaleDivider: 2.5,\\n saturnScaleDivider: 3.5,\\n uranusScaleDivider: 6.2,\\n neptuneScaleDivider: 8.7\\n}\\n\\n// This function sets up the simulation.\\nfunction init_T6() {\\n // One smart thing the original author of this simulation did was to split up the background and foreground into two canvases.\\n // That way, the static elements (the Sun and background) don't have to be rendered over and over, which improves performance.\\n T6.contextbackground = createGraphics(windowWidth, windowHeight);\\n T6.contextforeground = createGraphics(windowWidth, windowHeight);\\n\\n T6.interval = 1000 / T6.fps;\\n // This is all timing stuff. Our simulation begins with real data from January 2000, so we want to reflect that.\\n T6.current = new Date(Date.UTC(2000, 0, 2, 0, 0, 0));\\n T6.DAY = T6.current.getDate();\\n T6.MONTH = T6.current.getMonth() + 1;\\n T6.YEAR = T6.current.getFullYear();\\n T6.julianDate = getJulianDate_T6(T6.YEAR, T6.MONTH, T6.DAY);\\n T6.T = (T6.julianDate - T6.julianEpochJ2000) / T6.julianCenturyInJulianDays;\\n\\n // This gets the ball rolling.\\n renderBackground_T6();\\n run_T6();\\n}\\n\\n/*\\n This function runs the simulation.\\n It's not exactly readable.\\n The main idea is a bit different than what we have done before.\\n Instead of calculating forces and updating particles (planets, balls) one \"step\" at a time...\\n ... it uses orbital data (how long, wide, squished...) and calculates the next step along each orbit.\\n So it is almost as if the planets are just moving along their own tracks.\\n That's not as satisfying as if you calculate every force at every moment, but it lets us stay precise and accurate.\\n It's a tradeoff!\\n Adjust some numbers in the matrices above to get different orbits (but make sure to save the original delicate values!)\\n You gotta do what you gotta do.\\n\\n Learn more here: http://www.stargazing.net/kepler/ellipse.html\\n\\n (The asteroid system, a few functions below, is more like what we've done so far.)\\n*/\\nfunction run_T6() {\\n\\n // Read this chunk carefully to understand how time works in the simulation.\\n T6.now = Date.now();\\n T6.delta = T6.now - T6.then;\\n\\n if (T6.delta > T6.interval) {\\n T6.then = T6.now - (T6.delta % T6.interval);\\n\\n T6.FPSCount++;\\n T6.nowFPS = Date.now();\\n T6.deltaFPS = T6.nowFPS - T6.thenFPS;\\n if (T6.deltaFPS > 1000) {\\n T6.avgFPSCount = T6.FPSCount;\\n T6.FPSCount = 0;\\n T6.thenFPS = T6.nowFPS - (T6.deltaFPS % 1000);\\n }\\n\\n updateDate_T6(1);\\n\\n // Here's a whiff of our old thinking: X and Y components, separated!\\n T6.orbitalElements = plotPlanet_T6(T6.T, 0); T6.xMercury = T6.orbitalElements[0]; T6.yMercury = T6.orbitalElements[1];\\n T6.orbitalElements = plotPlanet_T6(T6.T, 1); T6.xVenus = T6.orbitalElements[0]; T6.yVenus = T6.orbitalElements[1];\\n T6.orbitalElements = plotPlanet_T6(T6.T, 2); T6.xEarth = T6.orbitalElements[0]; T6.yEarth = T6.orbitalElements[1];\\n T6.orbitalElements = plotPlanet_T6(T6.T, 3); T6.xMars = T6.orbitalElements[0]; T6.yMars = T6.orbitalElements[1];\\n T6.orbitalElements = plotPlanet_T6(T6.T, 4); T6.xJupiter = T6.orbitalElements[0]; T6.yJupiter = T6.orbitalElements[1];\\n T6.orbitalElements = plotPlanet_T6(T6.T, 5); T6.xSaturn = T6.orbitalElements[0]; T6.ySaturn = T6.orbitalElements[1];\\n T6.orbitalElements = plotPlanet_T6(T6.T, 6); T6.xUranus = T6.orbitalElements[0]; T6.yUranus = T6.orbitalElements[1];\\n T6.orbitalElements = plotPlanet_T6(T6.T, 7); T6.xNeptune = T6.orbitalElements[0]; T6.yNeptune = T6.orbitalElements[1];\\n\\n updateAsteroids();\\n }\\n\\n // TPS == 'ticks per second'\\n T6.TPSCount++;\\n T6.nowTPS = Date.now();\\n T6.deltaTPS = T6.nowTPS - T6.thenTPS;\\n if (T6.deltaTPS > 1000) {\\n T6.avgTPSCount = T6.TPSCount;\\n T6.TPSCount = 0;\\n T6.thenTPS = T6.nowTPS - (T6.deltaTPS % 1000);\\n }\\n\\n renderForeground_T6();\\n}\\n\\n// static elements\\nfunction renderBackground_T6() {\\n T6.contextbackground.background(15, 15, 50);\\n\\n if (sunImg) {\\n T6.contextbackground.image(sunImg, T6.width / 2 - 25, T6.height / 2 - 25, 50, 50);\\n } else {\\n console.error(\"Sun image not loaded\");\\n }\\n}\\n\\n// dynamic elements\\nfunction renderForeground_T6() {\\n T6.contextforeground.clear();\\n\\n // Error handling for images\\n if (!mercuryImg || !venusImg || !earthImg || !marsImg || !jupiterImg || !saturnImg || !uranusImg || !neptuneImg) {\\n console.error(\"One or more planet images not loaded\");\\n return;\\n }\\n\\n // Helper function that makes sure Saturn especially(!) isn't squished.\\n function drawImageWithAspectRatio(img, x, y, desiredWidth, desiredHeight) {\\n const aspectRatio = img.width / img.height;\\n let width = desiredWidth;\\n let height = desiredHeight;\\n\\n // Adjust width and height to preserve aspect ratio\\n if (width / height > aspectRatio) {\\n width = height * aspectRatio;\\n } else {\\n height = width / aspectRatio;\\n }\\n\\n T6.contextforeground.image(img, x - width / 2, y - height / 2, width, height);\\n }\\n\\n // Draw ALL the pictures.\\n drawImageWithAspectRatio(mercuryImg, (T6.width / 2) + (T6.xMercury * T6.scale), (T6.height / 2) - (T6.yMercury * T6.scale), 10, 10);\\n drawImageWithAspectRatio(venusImg, (T6.width / 2) + (T6.xVenus * T6.scale), (T6.height / 2) - (T6.yVenus * T6.scale), 20, 20);\\n drawImageWithAspectRatio(earthImg, (T6.width / 2) + (T6.xEarth * T6.scale), (T6.height / 2) - (T6.yEarth * T6.scale), 20, 20);\\n drawImageWithAspectRatio(marsImg, (T6.width / 2) + (T6.xMars * T6.scale), (T6.height / 2) - (T6.yMars * T6.scale), 14, 14);\\n drawImageWithAspectRatio(jupiterImg, (T6.width / 2) + (T6.xJupiter * T6.scale / T6.jupiterScaleDivider), (T6.height / 2) - (T6.yJupiter * T6.scale / T6.jupiterScaleDivider), 40, 40);\\n drawImageWithAspectRatio(saturnImg, (T6.width / 2) + (T6.xSaturn * T6.scale / T6.saturnScaleDivider), (T6.height / 2) - (T6.ySaturn * T6.scale / T6.saturnScaleDivider), 36, 36);\\n drawImageWithAspectRatio(uranusImg, (T6.width / 2) + (T6.xUranus * T6.scale / T6.uranusScaleDivider), (T6.height / 2) - (T6.yUranus * T6.scale / T6.uranusScaleDivider), 28, 28);\\n drawImageWithAspectRatio(neptuneImg, (T6.width / 2) + (T6.xNeptune * T6.scale / T6.neptuneScaleDivider), (T6.height / 2) - (T6.yNeptune * T6.scale / T6.neptuneScaleDivider), 28, 28);\\n\\n T6.contextforeground.fill(105, 105, 105);\\n T6.contextforeground.textSize(20);\\n T6.contextforeground.textAlign(LEFT);\\n T6.contextforeground.text(\"FPS = \" + T6.avgFPSCount, 10, 30);\\n T6.contextforeground.textAlign(LEFT);\\n T6.contextforeground.text(\"GREGORIAN DATE = \" + T6.YEAR + \".\" + T6.MONTH + \".\" + T6.DAY, 10, T6.height - 10);\\n}\\n\\n// Date stuff. If you care, look up the Julian versus Gregorian calendar.\\nfunction getJulianDate_T6(Year, Month, Day) {\\n var inputDate = new Date(Year, Month, Math.floor(Day));\\n var switchDate = new Date(\"1582\", \"10\", \"15\");\\n\\n var isGregorianDate = inputDate >= switchDate;\\n\\n if (Year < 0) {\\n Year++;\\n }\\n\\n if (Month == 1 || Month == 2) {\\n Year = Year - 1;\\n Month = Month + 12;\\n }\\n\\n var A = Math.floor(Year / 100);\\n var B = 2 - A + Math.floor(A / 4);\\n\\n if (!isGregorianDate) { B = 0; }\\n\\n return ((Math.floor(365.25 * Year)) + (Math.floor(30.6001 * (Month + 1))) + Day + 1720994.5 + B);\\n}\\n\\nfunction updateDate_T6(increment) {\\n T6.newDate = new Date(T6.current.getFullYear(), T6.current.getMonth(), T6.current.getDate() + increment);\\n T6.current = T6.newDate;\\n\\n newdd = T6.newDate.getDate();\\n newmm = T6.newDate.getMonth() + 1;\\n newyyyy = T6.newDate.getFullYear();\\n\\n T6.YEAR = newyyyy;\\n T6.MONTH = newmm;\\n T6.DAY = newdd;\\n\\n T6.julianDate = getJulianDate_T6(T6.YEAR, T6.MONTH, T6.DAY);\\n T6.T = (T6.julianDate - T6.julianEpochJ2000) / T6.julianCenturyInJulianDays;\\n}\\n\\n// This is the main function that controls each individual planet. This is the real physics!\\nfunction plotPlanet_T6(TGen, planetNumber) {\\n // a = semimajor axis, ec = eccentricity, E = eccentric anomaly (angle)\\n\\n // The recipe is the same: take the starting state of a planet, add the rate at which that changes times the time step, then update the state.\\n // THIS is actually familiar to you: it's the same as how velocity and acceleration work!\\n // This is calculus, baby! You can do it with anything that changes *continuously* — not just velocity.\\n aGen = T6.planetElements[planetNumber][0] + (T6.planetRates[planetNumber][0] * TGen);\\n eGen = T6.planetElements[planetNumber][1] + (T6.planetRates[planetNumber][1] * TGen);\\n iGen = T6.planetElements[planetNumber][2] + (T6.planetRates[planetNumber][2] * TGen);\\n // The % (modulus) symbol wraps values around a number. In this case, 360, for angles.\\n iGen = iGen % 360;\\n WGen = T6.planetElements[planetNumber][5] + (T6.planetRates[planetNumber][5] * TGen);\\n WGen = WGen % 360;\\n wGen = T6.planetElements[planetNumber][4] + (T6.planetRates[planetNumber][4] * TGen);\\n wGen = wGen % 360;\\n if (wGen < 0) { wGen = 360 + wGen; }\\n LGen = T6.planetElements[planetNumber][3] + (T6.planetRates[planetNumber][3] * TGen);\\n LGen = LGen % 360;\\n if (LGen < 0) { LGen = 360 + LGen; }\\n\\n MGen = LGen - (wGen);\\n if (MGen < 0) { MGen = 360 + MGen; }\\n\\n EGen = EccAnom_T6(eGen, MGen, 6);\\n\\n trueAnomalyArgGen = (Math.sqrt((1 + eGen) / (1 - eGen))) * (Math.tan(toRadians_T6(EGen) / 2));\\n\\n K = Math.PI / 180.0;\\n if (trueAnomalyArgGen < 0) {\\n nGen = 2 * (Math.atan(trueAnomalyArgGen) / K + 180);\\n }\\n else {\\n nGen = 2 * (Math.atan(trueAnomalyArgGen) / K)\\n }\\n\\n rGen = aGen * (1 - (eGen * (Math.cos(toRadians_T6(EGen)))));\\n\\n xGen = rGen * (Math.cos(toRadians_T6(WGen)) * Math.cos(toRadians_T6(nGen + wGen - WGen)) - Math.sin(toRadians_T6(WGen)) * Math.sin(toRadians_T6(nGen + wGen - WGen)) * Math.cos(toRadians_T6(iGen)));\\n yGen = rGen * (Math.sin(toRadians_T6(WGen)) * Math.cos(toRadians_T6(nGen + wGen - WGen)) + Math.cos(toRadians_T6(WGen)) * Math.sin(toRadians_T6(nGen + wGen - WGen)) * Math.cos(toRadians_T6(iGen)));\\n zGen = rGen * (Math.sin(toRadians_T6(nGen + wGen - WGen)) * Math.sin(toRadians_T6(iGen)));\\n\\n // At the end of the day, all we really care about is where the planet is: its X and its Y!\\n return [xGen, yGen];\\n}\\n\\n// The guy I took the above from took the below from: http://www.jgiesen.de/kepler/kepler.html\\n// Thanks to Juergen Giesen.\\n// The point of this (mathematical-as-heck!) function is to compute E, the *Eccentric Anomaly* of an orbit.\\n// That is a number that describes deviation from a standard elliptical (oval-shaped) path.\\n// The planets tilt and wobble in space, and we need to account for that to dial in our orbits.\\n// Lucky, some ancient dweebs figured out the math for us! Now we can just put it into code.\\nfunction EccAnom_T6(ec, m, dp) {\\n // ec = eccentricity, m = mean anomaly,\\n // dp = number of decimal places\\n\\n var pi = Math.PI, K = pi / 180.0;\\n var maxIter = 30, i = 0;\\n var delta = Math.pow(10, -dp);\\n var E, F;\\n\\n m = m / 360.0;\\n m = 2.0 * pi * (m - Math.floor(m));\\n\\n // Notice this funky style with no brackets.\\n if (ec < 0.8) E = m; else E = pi;\\n\\n F = E - ec * Math.sin(m) - m;\\n\\n // If nothing else, notice sine and cosine in here!\\n while ((Math.abs(F) > delta) && (i < maxIter)) {\\n E = E - F / (1.0 - ec * Math.cos(E));\\n F = E - ec * Math.sin(E) - m;\\n i = i + 1;\\n }\\n\\n E = E / K;\\n\\n return Math.round(E * Math.pow(10, dp)) / Math.pow(10, dp);\\n}\\n\\n// You can just use radians() to do this, but some people like to write their own helper functions.\\nfunction toRadians_T6(deg) {\\n return deg * (Math.PI / 180);\\n}\\n\\nfunction round_T6(value, decimals) {\\n return Number(Math.round(value + 'e' + decimals) + 'e-' + decimals);\\n}\\n\\n/*\\n All this asteroid stuff is my code, added to the simulation that draws on others' work.\\n The idea is not to be realisitc but rather to show how some of the physics ideas from before still apply.\\n Asteroids launch semi-randomly towards the sun from a click and the sun's gravity pulls them in.\\n They can even bounce off planets!\\n*/\\nfunction addAsteroid() {\\n const x = mouseX;\\n const y = mouseY;\\n const sunPos = createVector(width / 2, height / 2);\\n\\n // This makes a vector from the asteroid that points towards the Sun.\\n let direction = p5.Vector.sub(sunPos, createVector(x, y));\\n\\n // Then rotate that vector by up to 60º so that it doesn't fly directly into the Sun.\\n const angleOffset = radians(random(-60, 60));\\n direction.rotate(angleOffset);\\n\\n // We give each new asteroid a random velocity between a min and a max.\\n const minVelocity = 0;\\n const maxVelocity = 20;\\n const velocityMagnitude = random(minVelocity, maxVelocity);\\n\\n const velocity = direction.setMag(velocityMagnitude);\\n\\n // This is the asteroid object, that gets added to the list.\\n // All you need is position, velocity, mass, and a trail for a nice visual effect.\\n asteroids.push({\\n position: createVector(x, y),\\n velocity: velocity,\\n mass: asteroidMass,\\n trail: [] // The trail will be an array of points.\\n });\\n console.log(`Asteroid added at position: (${x}, ${y}) with velocity magnitude: ${velocityMagnitude}`);\\n}\\n\\n// Once an asteroid is born, it needs to move!\\nfunction updateAsteroids() {\\n const G = 6.67430e-11; // A God-given number representing the strength of gravity.\\n\\n for (let asteroid of asteroids) {\\n // Apply gravity from the Sun\\n let sunPos = createVector(width / 2, height / 2);\\n let force = p5.Vector.sub(sunPos, asteroid.position);\\n // We calculate the square of the distance to the Sun because gravity depends on distance squared, not just distance.\\n let distanceSq = max(force.magSq(), 1000000); // Don't want this number to get too big.\\n // This is a force, and therefore acceleration, calculation.\\n let strength = (G * asteroid.mass * 1.989e30) / distanceSq; // Sun's mass: 1.989e30 kg\\n force.setMag(strength * 1e-29); // We move the decimal point 29 times to scale down to the pixel level rather than astronomical units.\\n // The key moment: add to the velocity. That's what acceleration (force) does.\\n asteroid.velocity.add(force);\\n\\n // We do the same thing with each planet (doesn't effect the asteroid that much tbh... could comment this out.)\\n let planets = [\\n { x: T6.xMercury * T6.scale, y: T6.yMercury * T6.scale, mass: 3.3011e23 },\\n { x: T6.xVenus * T6.scale, y: T6.yVenus * T6.scale, mass: 4.8675e24 },\\n { x: T6.xEarth * T6.scale, y: T6.yEarth * T6.scale, mass: 5.97237e24 },\\n { x: T6.xMars * T6.scale, y: T6.yMars * T6.scale, mass: 6.4171e23 },\\n { x: T6.xJupiter * T6.scale / T6.jupiterScaleDivider, y: T6.yJupiter * T6.scale / T6.jupiterScaleDivider, mass: 1.8982e27 },\\n { x: T6.xSaturn * T6.scale / T6.saturnScaleDivider, y: T6.ySaturn * T6.scale / T6.saturnScaleDivider, mass: 5.6834e26 },\\n { x: T6.xUranus * T6.scale / T6.uranusScaleDivider, y: T6.yUranus * T6.scale / T6.uranusScaleDivider, mass: 8.6810e25 },\\n { x: T6.xNeptune * T6.scale / T6.neptuneScaleDivider, y: T6.yNeptune * T6.scale / T6.neptuneScaleDivider, mass: 1.02413e26 }\\n ];\\n\\n for (let planet of planets) {\\n let planetPos = createVector(width / 2 + planet.x, height / 2 - planet.y);\\n let force = p5.Vector.sub(planetPos, asteroid.position);\\n let distanceSq = max(force.magSq(), 1000000);\\n let strength = (G * asteroid.mass * planet.mass) / distanceSq;\\n force.setMag(strength * 1e-29);\\n asteroid.velocity.add(force);\\n\\n // Check for collision with a Planet\\n if (asteroid.position.dist(planetPos) < 20) {\\n // Conservation of momentum!\\n const planetVelocity = createVector(0, 0); // Planets are so big relative to asteroids, it's like they aren't even moving (this keeps our math simple).\\n // This is a version of the momentum calculatino from a week ago.\\n const newVelocity = p5.Vector.add(\\n p5.Vector.mult(asteroid.velocity, asteroid.mass),\\n p5.Vector.mult(planetVelocity, planet.mass)\\n ).div(asteroid.mass + planet.mass);\\n asteroid.velocity = newVelocity;\\n //asteroid.toRemove = true; // Optional: remove asteroid after collision\\n break;\\n }\\n }\\n\\n // Add velocity to position. Back to day one!\\n asteroid.position.add(asteroid.velocity);\\n\\n // The trail is the last 50 points where the asteroid was for each frame.\\n asteroid.trail.push(asteroid.position.copy());\\n if (asteroid.trail.length > 50) {\\n asteroid.trail.shift(); // Remove the oldest point if the trail is longer than 50\\n }\\n\\n // Check for collision with the Sun: remove the asteroid.\\n if (asteroid.position.dist(sunPos) < 20) { // Sun's radius in pixels\\n asteroid.toRemove = true;\\n continue;\\n }\\n }\\n\\n // Remove asteroids that are off the canvas or marked for removal\\n asteroids = asteroids.filter(asteroid => {\\n const isOffCanvas = asteroid.position.x < -200 || asteroid.position.x > width + 200 || asteroid.position.y < -200 || asteroid.position.y > height + 200;\\n if (isOffCanvas) {\\n //console.log(`Asteroid left the canvas at position (${asteroid.position.x}, ${asteroid.position.y})`);\\n }\\n return !asteroid.toRemove && !asteroid.isOffCanvas;\\n });\\n}\\n\\nfunction drawAsteroids() {\\n // Draw asteroids and their trails\\n for (let asteroid of asteroids) {\\n // Draw trail\\n for (let i = 0; i < asteroid.trail.length; i++) {\\n let trailPos = asteroid.trail[i];\\n fill(255, 255, 255, map(i, 0, asteroid.trail.length, 0, 255)); // Fade trail points\\n noStroke();\\n ellipse(trailPos.x, trailPos.y, 2, 2); // Draw each point in the trail as a small circle\\n }\\n\\n // Draw asteroid\\n fill(255);\\n noStroke();\\n ellipse(asteroid.position.x, asteroid.position.y, 10, 10); // Draw asteroid as a circle\\n }\\n}\\n\\n///////////////////////////////////////////////////////////////////////////////////////////\\n\\n/*\\n You may want to extend this project. Go for it!\\n Add stars to the background.\\n Code up an astronaut or a spaceship that can fly around.\\n Cook up a speculative new planet!\\n \\n Anything you want. Have fun!\\n*/\\n\\n///////////////////////////////////////////////////////////////////////////////////////////\\n\\n/*\\n We have come so far. Give yourselves a hand!\\n A standard physics course will last a whole year in high school or a semester in college, and mostly has you solving math problems by hand.\\n That can be fun indeed. Math is fun.\\n But it's a whole new world, baby!\\n We can *code*.\\n\\n What I want you to come away with:\\n - Some physics *vocabulary* so you can describe what you see and feel.\\n - Some knowledge of the kinds of mathematical and computation patterns that make physics work.\\n - The confidence that you can build terrific projects fast with help from others:\\n Scientists from the past, coders from the present, and your peers of course.\\n - Most of all, curiosity to learn more!\\n \\n Where to go next in physics? I can see a few paths:\\n - Put physics into a game project! Make characters and objects move in interesting ways. That's a whole career's worth of mastery.\\n - Make some art — creative coding — using physics and math ideas.\\n - Study basic physics the \"old school\" way. Read a book, do some problems, watch some lectures. Challenge yourself to code those ideas as a way to master them.\\n - Learn some more math to support that journey. Algebra and Geometry are the father and mother of physics. Those ingredients combine into Calculus, the most powerful physics tool of all.\\n - Explore other topics in phsyics:\\n - Work and Energy\\n - Electricity and Magnetism\\n - Light and optics\\n - Other forces: tension, friction, springs, ...\\n - Nuclear physics\\n There's a lot to learn!\\n\\n My request: drop me a line when you do something cool. support@jippity.pro\\n\\n Proud of your work; can't wait to see where it takes you.\\n*/","style.css":"body {\\n color:#c32de05e\\n}","index.html":"\\n\\n\\n\\n \\n \\n Simulate the Solar System FOR REAL\\n \\n\\n\\n\\n \\n

\\n\\n \\n

\\n \\n

\\n\\n

\\n \\n \\n\\n\\n"},"chat_history":[],"assets":[{"project_id":2108,"creator":null,"cloudinary_id":"lzam3qxjwjkomxvifq5r","deleted":false,"id":11590,"created_at":"2024-07-01T14:44:04.367267Z","name":"SS","resource_type":"image","tags":null},{"project_id":2108,"creator":null,"cloudinary_id":"v9monl3gqcd7jdhzewxo","deleted":false,"id":11589,"created_at":"2024-07-01T14:06:12.769161Z","name":"mercury","resource_type":"image","tags":null},{"project_id":2108,"creator":null,"cloudinary_id":"i6exysrb54inn0dbzz5e","deleted":false,"id":11588,"created_at":"2024-07-01T14:06:12.768166Z","name":"venus","resource_type":"image","tags":null},{"project_id":2108,"creator":null,"cloudinary_id":"rsoiws15wjy50gbdsiiv","deleted":false,"id":11587,"created_at":"2024-07-01T14:06:12.766884Z","name":"earth","resource_type":"image","tags":null},{"project_id":2108,"creator":null,"cloudinary_id":"dyzkeopx9xzgqdzwhhun","deleted":false,"id":11586,"created_at":"2024-07-01T14:06:12.765645Z","name":"mars","resource_type":"image","tags":null},{"project_id":2108,"creator":null,"cloudinary_id":"knztmajsviml6ydjrjbx","deleted":false,"id":11585,"created_at":"2024-07-01T14:06:12.764266Z","name":"jupiter","resource_type":"image","tags":null},{"project_id":2108,"creator":null,"cloudinary_id":"cp7elsfrletlgcrpc2mw","deleted":false,"id":11584,"created_at":"2024-07-01T14:06:12.763034Z","name":"saturn","resource_type":"image","tags":null},{"project_id":2108,"creator":null,"cloudinary_id":"zkpxbh9u1z2gqqd50rku","deleted":false,"id":11583,"created_at":"2024-07-01T14:06:12.761742Z","name":"uranus","resource_type":"image","tags":null},{"project_id":2108,"creator":null,"cloudinary_id":"zun3dpkcjgdgd4jwqwc1","deleted":false,"id":11582,"created_at":"2024-07-01T14:06:12.760171Z","name":"neptune","resource_type":"image","tags":null},{"project_id":2108,"creator":null,"cloudinary_id":"m2fnyrzanyad18imd8nr","deleted":false,"id":11581,"created_at":"2024-07-01T14:06:12.758238Z","name":"sun","resource_type":"image","tags":null}],"messages_sent":23,"comments":[],"steps":[]}