Image:Gravity well plot.svg

[edit] Summary

Graph of gravity well. Above: Distortion of space-time. Below: Gravitational potential function.

Asymptote source code:

import graph3;

typedef real fun(real);

fun Phi(real R, real GM_R, real shift=0) {
  return new real(real r) {
    real r_R = abs(r) / R;
    if (r_R >= 1)
      return -GM_R / r_R + shift;
    else
      return 0.5 GM_R * (r_R^2 - 3) + shift;
  };
}

typedef real fun2(pair);
fun2 revolve(fun f) { return new real(pair z) { return f(abs(z)); }; }

size(300,0);
currentprojection=perspective(0,20,4);
add(surface(revolve(Phi(2,7)), nsub=4, (-12,-12), (12,12), nx=20, light=O));
draw(graph(Phi(2,7,-10),-24,24));

[edit] Licensing

I, the copyright holder of this work, hereby release it into the public domain. This applies worldwide. In case this is not legally possible: I grant anyone the right to use this work for any purpose, without any conditions, unless such conditions are required by law. Afrikaans | Alemannisch | Aragonés | العربية | Asturianu | Български | Català | Česky | Cymraeg | Dansk | Deutsch | Eʋegbe | Ελληνικά | English | Español | Esperanto | Euskara | Estremeñu | فارسی | Français | Galego | 한국어 | हिन्दी | Hrvatski | Ido | Bahasa Indonesia | Íslenska | Italiano | עברית | Kurdî / كوردی | Latina | Lietuvių | Latviešu | Magyar | Македонски | Bahasa Melayu | Nederlands | ‪Norsk (bokmål)‬ | ‪Norsk (nynorsk)‬ | 日本語 | Polski | Português | Ripoarisch | Română | Русский | Shqip | Slovenčina | Slovenščina | Српски / Srpski | Svenska | ไทย | Tagalog | Türkçe | Українська | Tiếng Việt | Walon | ‪中文(简体)‬ | ‪中文(繁體)‬ | zh-yue-hant | +/-

• Gravitational well