gravitational response of a ball

import%20marimo%0A%0A__generated_with%20%3D%20%220.14.15%22%0Aapp%20%3D%20marimo.App(width%3D%22medium%22)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_()%3A%0A%20%20%20%20import%20marimo%20as%20mo%0A%20%20%20%20import%20array%20as%20ar%0A%20%20%20%20import%20matplotlib.pylab%20as%20plt%0A%20%20%20%20import%20math%0A%20%20%20%20from%20collections%20import%20namedtuple%2C%20UserList%0A%20%20%20%20import%20plotly.graph_objects%20as%20go%0A%20%20%20%20from%20plotly.subplots%20import%20make_subplots%0A%20%20%20%20return%20UserList%2C%20ar%2C%20go%2C%20math%2C%20mo%2C%20namedtuple%2C%20plt%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.hstack(%5B%0A%0A%20%20%20%20%20%20%20%20mo.vstack(%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20rf%22%22%22%0A%20%20%20%20%23%20Gravitational%20Responses%20of%202%20%22imaginary%22%20perfect%20sphere%20objects%20below%20surface%0A%0A%20%20%20%20A%20ball%20with%20density%20of%20%5C(%5Crho%5C)%20located%20at%20%5C(x%2Cy%2Cz%5C)%20below%20the%20surface%20will%20induce%20a%20vector%20field%0A%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(r%22%22%22%23%20%24%5Cvec%7Bg%7D(F_1)%20%3D%20-%5Cgamma%20%5Cfrac%7B%5Cfrac%7B4%7D%7B3%7D%20%5Cpi%20a%5E2%20%5Crho%7D%7BR%5E2%7D%20%5Chat%7Br%7D%24%22%22%22).center()%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(r%22%22%22**Where%3A**%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20-%20%5C(%5Cvec%7Bg%7D(F_1)%5C)%20%3D%20Gravitational%20field%20at%20point%20%5C(F_1%5C)%20due%20to%20a%20mass%20distribution%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20-%20%5C(%5Cgamma%5C)%20%3D%20Gravitational%20constant%20(also%20written%20as%20%5C(G%5C))%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20-%20%5C(a%5C)%20%3D%20Radius%20of%20the%20ball%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20-%20%5C(%5Crho%5C)%20%3D%20Mass%20density%20of%20the%20object%20(mass%20per%20unit%20volume%20or%20surface%20area%20depending%20on%20context)%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20-%20%5C(R%5C)%20%3D%20Distance%20from%20the%20center%20of%20the%20mass%20distribution%20to%20the%20point%20%5C(F_1%5C)%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20-%20%5C(%5Chat%7Br%7D%5C)%20%3D%20Unit%20vector%20pointing%20from%20the%20mass%20distribution%20to%20the%20point%20%5C(F_1%5C)%20(direction%20of%20the%20field)%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20-%20%5C(%5Cfrac%7B4%7D%7B3%7D%20%5Cpi%20a%5E2%5C)%20%3D%20Volume%20of%20a%20ball%0A%20%20%20%20%20%20%20%20%20%20%20%20%22%22%22)%5D%2C%20align%3D%22center%22)%2C%0A%20%20%20%20%20%20%20%20mo.image(%0A%20%20%20%20%20%20%20%20%20%20%20%20src%3Dstr(mo.notebook_location()%20%2F%20%22public%2Fgravity.jpeg%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20width%3D%22800px%22%2C%20height%3D%22320px%22%2C%20rounded%3DTrue%0A%20%20%20%20%20%20%20%20).center()%2C%0A%20%20%20%20%20%20%20%20mo.md(%22%5Cn---%22)%0A%20%20%20%20%5D%2C%20justify%3D%22start%22%2C%20align%3D%22center%22).center()%0A%0A%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20r%22%22%22%0A%20%20%20%20---%0A%0A%20%20%20%20%23%20The%20Exciting%20Part%0A%20%20%20%20%22%22%22%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(math)%3A%0A%20%20%20%20GAMMA%20%3D%206.6E-11%0A%20%20%20%20SI2MG%20%3D%201.0E5%0A%20%20%20%20PI%20%3D%20math.pi%0A%20%20%20%20KM2M%20%3D%201e3%20%23%20konversi%20satuan%0A%20%20%20%20IERROR%20%3D%200%0A%20%20%20%20NX%20%3D%20200%0A%20%20%20%20NY%20%3D%20200%0A%0A%20%20%20%20LEFT%20%3D%20-50%0A%20%20%20%20BOTTOM%20%3D%20-50%0A%20%20%20%20RIGTH%20%3D%2050%0A%20%20%20%20TOP%20%3D%2050%0A%20%20%20%20return%20BOTTOM%2C%20GAMMA%2C%20KM2M%2C%20LEFT%2C%20NX%2C%20NY%2C%20PI%2C%20RIGTH%2C%20SI2MG%2C%20TOP%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(namedtuple)%3A%0A%20%20%20%20Ball%20%3D%20namedtuple(%22Ball%22%2C%20%5B%22r%22%2C%20%22rho%22%2C%20%22x%22%2C%20%22y%22%2C%20%22z%22%5D%2C%20defaults%3D%5B10.0%2C%2090.0%2C%200%2C%200%2C%200%5D)%0A%20%20%20%20return%20(Ball%2C)%0A%0A%0A%40app.cell%0Adef%20_(BOTTOM%2C%20LEFT%2C%20NX%2C%20NY%2C%20RIGTH%2C%20TOP%2C%20linspace%2C%20meshgrid)%3A%0A%20%20%20%20x%20%3D%20linspace(LEFT%2C%20RIGTH%2C%20NX)%0A%20%20%20%20y%20%3D%20linspace(BOTTOM%2C%20TOP%2C%20NY)%0A%20%20%20%20xx%2C%20yy%20%3D%20meshgrid(x%2C%20y)%0A%20%20%20%20return%20xx%2C%20yy%0A%0A%0A%40app.cell%0Adef%20_(%0A%20%20%20%20Ball%2C%0A%20%20%20%20NX%2C%0A%20%20%20%20NY%2C%0A%20%20%20%20density_ball_a%2C%0A%20%20%20%20density_ball_b%2C%0A%20%20%20%20myArr%2C%0A%20%20%20%20radius_ball_a%2C%0A%20%20%20%20radius_ball_b%2C%0A%20%20%20%20slider_posx_ball_a%2C%0A%20%20%20%20slider_posx_ball_b%2C%0A%20%20%20%20slider_posy_ball_a%2C%0A%20%20%20%20slider_posy_ball_b%2C%0A%20%20%20%20slider_posz_ball_a%2C%0A%20%20%20%20slider_posz_ball_b%2C%0A%20%20%20%20zeros%2C%0A)%3A%0A%20%20%20%20ball_a%20%3D%20Ball(r%3Dradius_ball_a.value%2C%20rho%3Ddensity_ball_a.value%2C%20x%3Dslider_posx_ball_a.value%2C%20y%3Dslider_posy_ball_a.value%2C%20z%3Dslider_posz_ball_a.value)%0A%20%20%20%20ball_b%20%3D%20Ball(r%3Dradius_ball_b.value%2C%20rho%3Ddensity_ball_b.value%2C%20x%3Dslider_posx_ball_b.value%2C%20y%3Dslider_posy_ball_b.value%2C%20z%3Dslider_posz_ball_b.value)%0A%20%20%20%20G1%20%20%3DmyArr(zeros((NX%2CNY)))%0A%20%20%20%20G1_a%3DmyArr(zeros((NX%2CNY)))%0A%20%20%20%20G1_b%3DmyArr(zeros((NX%2CNY)))%0A%20%20%20%20return%20G1%2C%20G1_a%2C%20G1_b%2C%20ball_a%2C%20ball_b%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(UserList)%3A%0A%20%20%20%20class%20myArr(UserList)%3A%0A%20%20%20%20%20%20%20%20ndim%20%3D%202%0A%20%20%20%20return%20(myArr%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(G1%2C%20G1_a%2C%20G1_b%2C%20NX%2C%20NY%2C%20ball_a%2C%20ball_b%2C%20bola%2C%20math%2C%20xx%2C%20yy)%3A%0A%20%20%20%20MIN%20%3D%20math.inf%0A%20%20%20%20MAX%20%3D%20-math.inf%0A%0A%20%20%20%20MIN_a%20%3D%20math.inf%0A%20%20%20%20MAX_a%20%3D%20-math.inf%0A%0A%20%20%20%20MIN_b%20%3D%20math.inf%0A%20%20%20%20MAX_b%20%3D%20-math.inf%0A%0A%20%20%20%20for%20i%20in%20range(NX)%3A%0A%20%20%20%20%20%20%20%20for%20j%20in%20range(NY)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20gx1%2C%20gy1%2C%20gz1%20%3D%20bola(ball_a%2C%20xx%5Bi%5D%5Bj%5D%2C%20yy%5Bi%5D%5Bj%5D%2C%20-2)%0A%20%20%20%20%20%20%20%20%20%20%20%20gx2%2C%20gy2%2C%20gz2%20%3D%20bola(ball_b%2C%20xx%5Bi%5D%5Bj%5D%2C%20yy%5Bi%5D%5Bj%5D%2C%20-2)%0A%20%20%20%20%20%20%20%20%20%20%20%20G1%5Bi%5D%5Bj%5D%3Dgz1%2Bgz2%0A%20%20%20%20%20%20%20%20%20%20%20%20G1_a%5Bi%5D%5Bj%5D%3Dgz1%0A%20%20%20%20%20%20%20%20%20%20%20%20G1_b%5Bi%5D%5Bj%5D%3Dgz2%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20MIN%20%3E%20G1%5Bi%5D%5Bj%5D%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20MIN%20%3D%20G1%5Bi%5D%5Bj%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20MAX%20%3C%20G1%5Bi%5D%5Bj%5D%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20MAX%20%3D%20G1%5Bi%5D%5Bj%5D%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20MIN_a%20%3E%20G1_a%5Bi%5D%5Bj%5D%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20MIN_a%20%3D%20G1_a%5Bi%5D%5Bj%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20MAX_a%20%3C%20G1_a%5Bi%5D%5Bj%5D%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20MAX_a%20%3D%20G1_a%5Bi%5D%5Bj%5D%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20MIN_b%20%3E%20G1_b%5Bi%5D%5Bj%5D%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20MIN_b%20%3D%20G1_b%5Bi%5D%5Bj%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20MAX_b%20%3C%20G1_b%5Bi%5D%5Bj%5D%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20MAX_b%20%3D%20G1_b%5Bi%5D%5Bj%5D%0A%0A%20%20%20%20return%20MAX%2C%20MAX_a%2C%20MAX_b%2C%20MIN%2C%20MIN_a%2C%20MIN_b%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(input_ball_a%2C%20input_ball_b%2C%20mo)%3A%0A%20%20%20%20mo.hstack(%5Binput_ball_a%2C%20input_ball_b%5D)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(G1%2C%20G1_a%2C%20G1_b%2C%20MAX%2C%20MAX_a%2C%20MAX_b%2C%20MIN%2C%20MIN_a%2C%20MIN_b%2C%20mo%2C%20plt%2C%20xx%2C%20yy)%3A%0A%20%20%20%20fig%2C%20axes%20%3D%20plt.subplot_mosaic(%5B%5B%22a%22%2C%20%22a%2Bb%22%2C%20%22b%22%5D%5D%2C%20figsize%3D(15%2C5)%2C%20sharey%3DTrue%2C%20layout%3D%22tight%22)%0A%20%20%20%20fig.suptitle(%22Gravitational%20Field%22%2C%20fontsize%3D15)%0A%0A%20%20%20%20contour_a_b%20%3D%20axes%5B%22a%2Bb%22%5D.contourf(xx%2C%20yy%2C%20G1%2C%20cmap%3D%22jet%22%2C%20vmin%3DMIN%2C%20vmax%3DMAX)%0A%20%20%20%20axes%5B%22a%2Bb%22%5D.set_title(%22Combine%22%2C%20fontsize%3D13%2C%20loc%3D%22left%22)%0A%20%20%20%20axes%5B%22a%2Bb%22%5D.grid(alpha%3D0.25)%0A%20%20%20%20cbar_a_b%20%3D%20fig.colorbar(contour_a_b%2C%20ax%3Daxes%5B%22a%2Bb%22%5D)%0A%20%20%20%20cbar_a_b.set_label('Gravity%20Response')%0A%0A%20%20%20%20contour_a%20%3D%20axes%5B%22a%22%5D.contourf(xx%2C%20yy%2C%20G1_a%2C%20cmap%3D%22coolwarm%22%2C%20vmin%3DMIN_a%2C%20vmax%3DMAX_a)%0A%20%20%20%20axes%5B%22a%22%5D.set_title(%22Ball%20A%22%2C%20fontsize%3D13%2C%20loc%3D%22left%22)%0A%20%20%20%20axes%5B%22a%22%5D.grid(alpha%3D0.25)%0A%20%20%20%20cbar_a%20%3D%20fig.colorbar(contour_a%2C%20ax%3Daxes%5B%22a%22%5D)%0A%20%20%20%20cbar_a.set_label('Gravity%20Response')%0A%0A%0A%20%20%20%20contour_b%20%3D%20axes%5B%22b%22%5D.contourf(xx%2C%20yy%2C%20G1_b%2C%20cmap%3D%22coolwarm%22%2C%20vmin%3DMIN_b%2C%20vmax%3DMAX_b)%0A%20%20%20%20axes%5B%22b%22%5D.set_title(%22Ball%20B%22%2C%20fontsize%3D13%2C%20loc%3D%22left%22)%0A%20%20%20%20axes%5B%22b%22%5D.grid(alpha%3D0.25)%0A%20%20%20%20cbar_b%20%3D%20fig.colorbar(contour_b%2C%20ax%3Daxes%5B%22b%22%5D)%0A%20%20%20%20cbar_b.set_label('Gravity%20Response')%0A%0A%20%20%20%20for%20_%2C%20_ax%20in%20axes.items()%3A%0A%20%20%20%20%20%20%20%20_ax.set(ylabel%3D%22Y%22%2C%20xlabel%3D%22X%22)%0A%0A%20%20%20%20mo.as_html(plt.gca()).center()%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(%0A%20%20%20%20density_ball_a%2C%0A%20%20%20%20mo%2C%0A%20%20%20%20radius_ball_a%2C%0A%20%20%20%20slider_posx_ball_a%2C%0A%20%20%20%20slider_posy_ball_a%2C%0A%20%20%20%20slider_posz_ball_a%2C%0A)%3A%0A%20%20%20%20input_ball_a%20%3D%20mo.vstack(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(f%22%22%22%23%23%20Ball%20A%22%22%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(f%22%22%22%0A%20%20%20%20-%20the%20larger%20Z%20the%20deeper%20the%20ball%20below%20the%20surface%0A%20%20%20%20%22%22%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.hstack(%5Bslider_posx_ball_a%2C%20slider_posy_ball_a%2C%20slider_posz_ball_a%5D%2C%20justify%3D%22start%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.hstack(%5Bradius_ball_a%2C%20density_ball_a%5D%2C%20justify%3D%22start%22)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%23.callout()%0A%20%20%20%20return%20(input_ball_a%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(%0A%20%20%20%20density_ball_b%2C%0A%20%20%20%20mo%2C%0A%20%20%20%20radius_ball_b%2C%0A%20%20%20%20slider_posx_ball_b%2C%0A%20%20%20%20slider_posy_ball_b%2C%0A%20%20%20%20slider_posz_ball_b%2C%0A)%3A%0A%20%20%20%20input_ball_b%20%3D%20mo.vstack(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(f%22%22%22%23%23%20Ball%20B%22%22%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(f%22%22%22%0A%20%20%20%20-%20the%20larger%20Z%20the%20deeper%20the%20ball%20below%20the%20surface%0A%20%20%20%20%22%22%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.hstack(%5Bslider_posx_ball_b%2C%20slider_posy_ball_b%2C%20slider_posz_ball_b%5D%2C%20justify%3D%22start%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.hstack(%5Bradius_ball_b%2C%20density_ball_b%5D%2C%20justify%3D%22start%22)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%23.callout()%0A%20%20%20%20return%20(input_ball_b%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(ball_a%2C%20ball_b%2C%20mo)%3A%0A%20%20%20%20stats%20%3D%20mo.hstack(%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.stat(label%3D%22Ball%20A%22%2C%20value%3Dstr(ball_a)%2C%20bordered%3DTrue)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.stat(label%3D%22Ball%20B%22%2C%20value%3Dstr(ball_b)%2C%20bordered%3DTrue)%0A%20%20%20%20%20%20%20%20%5D%2C%20justify%3D%22center%22)%0A%0A%20%20%20%20stats.center()%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(BOTTOM%2C%20LEFT%2C%20RIGTH%2C%20TOP%2C%20mo)%3A%0A%20%20%20%20slider_posx_ball_a%20%3D%20mo.ui.slider(label%3D%22%60X%60%22%2C%20value%3D%2020%2C%20start%3DLEFT%2C%20stop%3DRIGTH%2C%20step%3D1)%0A%20%20%20%20slider_posy_ball_a%20%3D%20mo.ui.slider(label%3D%22%60Y%60%22%2C%20value%3D-20%2Cstart%3DBOTTOM%2C%20stop%3DTOP%2C%20step%3D1)%0A%20%20%20%20slider_posz_ball_a%20%3D%20mo.ui.slider(label%3D%22%60Z%60%22%2C%20value%3D%20%202%2Cstart%3D0.5%2C%20stop%3D10_000%2C%20step%3D0.5)%0A%20%20%20%20radius_ball_a%20%3D%20mo.ui.number(label%3D%22Radius%22%2C%20start%20%3D%201%2C%20stop%20%3D%2050%2C%20step%3D0.5%2C%20value%3D12)%0A%20%20%20%20density_ball_a%20%3D%20mo.ui.number(label%3Dr%22(%24%5Crho%24)%22%2C%20start%20%3D%201%2C%20stop%20%3D%2050_000%2C%20step%3D0.5%2C%20value%3D10000.10)%0A%0A%0A%20%20%20%20slider_posx_ball_b%20%3D%20mo.ui.slider(label%3D%22%60X%60%22%2C%20value%3D-20%2C%20start%3DLEFT%2C%20stop%3DRIGTH%2C%20step%3D1)%0A%20%20%20%20slider_posy_ball_b%20%3D%20mo.ui.slider(label%3D%22%60Y%60%22%2C%20value%3D20%2Cstart%3DBOTTOM%2C%20stop%3DTOP%2C%20step%3D1)%0A%20%20%20%20slider_posz_ball_b%20%3D%20mo.ui.slider(label%3D%22%60Z%60%22%2C%20value%3D%20%202%2Cstart%3D0.5%2C%20stop%3D10_000%2C%20step%3D0.5)%0A%20%20%20%20radius_ball_b%20%3D%20mo.ui.number(label%3D%22Radius%22%2C%20start%20%3D%201%2C%20stop%20%3D%2050%2C%20step%3D0.5%2C%20value%3D20)%0A%20%20%20%20density_ball_b%20%3D%20mo.ui.number(label%3Dr%22(%24%5Crho%24)%22%2C%20start%20%3D%201%2C%20stop%20%3D%2050_000%2C%20step%3D0.5%2C%20value%3D12000.10)%0A%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20density_ball_a%2C%0A%20%20%20%20%20%20%20%20density_ball_b%2C%0A%20%20%20%20%20%20%20%20radius_ball_a%2C%0A%20%20%20%20%20%20%20%20radius_ball_b%2C%0A%20%20%20%20%20%20%20%20slider_posx_ball_a%2C%0A%20%20%20%20%20%20%20%20slider_posx_ball_b%2C%0A%20%20%20%20%20%20%20%20slider_posy_ball_a%2C%0A%20%20%20%20%20%20%20%20slider_posy_ball_b%2C%0A%20%20%20%20%20%20%20%20slider_posz_ball_a%2C%0A%20%20%20%20%20%20%20%20slider_posz_ball_b%2C%0A%20%20%20%20)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(G1%2C%20G1_a%2C%20G1_b%2C%20go%2C%20mo%2C%20xx%2C%20yy)%3A%0A%20%20%20%20_xx%20%3D%20%5Bx.tolist()%20for%20x%20in%20xx%5D%0A%20%20%20%20_yy%20%3D%20%5By.tolist()%20for%20y%20in%20yy%5D%0A%0A%20%20%20%20_G1_a%20%3D%20%5Bx.tolist()%20for%20x%20in%20G1_a%5D%0A%20%20%20%20_G1%20%3D%20%5Bx.tolist()%20for%20x%20in%20G1%5D%0A%20%20%20%20_G1_b%20%3D%20%5Bx.tolist()%20for%20x%20in%20G1_b%5D%0A%0A%0A%20%20%20%20%23%20Create%203D%20surface%20plot%0A%20%20%20%20fig3d%20%3D%20go.Figure(data%3D%5Bgo.Surface(%0A%20%20%20%20%20%20%20%20x%3D_xx%2C%20%20%23%20assuming%20xx%20is%20your%20meshgrid%20x-coordinates%0A%20%20%20%20%20%20%20%20y%3D_yy%2C%20%20%23%20assuming%20yy%20is%20your%20meshgrid%20y-coordinates%0A%20%20%20%20%20%20%20%20z%3D_G1%2C%0A%20%20%20%20%20%20%20%20colorscale%3D'Jet'%2C%0A%20%20%20%20%20%20%20%20showscale%3DTrue%2C%0A%20%20%20%20%20%20%20%20colorbar%3Ddict(%0A%20%20%20%20%20%20%20%20%20%20%20%20x%3D0.85%2C%20%20%23%20position%20of%20colorbar%0A%20%20%20%20%20%20%20%20%20%20%20%20len%3D0.75%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20)%5D)%0A%0A%20%20%20%20%23%20Update%20layout%20with%20camera%20view%20and%20labels%0A%20%20%20%20fig3d.update_layout(%0A%20%20%20%20%20%20%20%20autosize%3DTrue%2C%0A%20%20%20%20%20%20%20%20width%3D1200%2C%0A%20%20%20%20%20%20%20%20height%3D800%2C%0A%20%20%20%20%20%20%20%20scene%3Ddict(%0A%20%20%20%20%20%20%20%20%20%20%20%20xaxis_title%3D'X'%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20yaxis_title%3D'Y'%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20zaxis_title%3D'Gravity'%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20camera%3Ddict(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20center%3Ddict(x%3D0%2C%20y%3D0%2C%20z%3D0)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20eye%3Ddict(x%3D1.25%2C%20y%3D1.25%2C%20z%3D1.25)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20up%3Ddict(x%3D0%2C%20y%3D0%2C%20z%3D1)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20aspectratio%3Ddict(x%3D1%2C%20y%3D1%2C%20z%3D1)%20%20%23%20similar%20to%20%5C%22auto%5C%22%20aspect%2C%0A%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20margin%3Ddict(l%3D0%2C%20r%3D0%2C%20b%3D0%2C%20t%3D0)%0A%20%20%20%20)%0A%0A%20%20%20%20fig3d.update_scenes(%0A%20%20%20%20%20%20%20%20xaxis_showbackground%3DFalse%2C%0A%20%20%20%20%20%20%20%20yaxis_showbackground%3DFalse%2C%0A%20%20%20%20%20%20%20%20zaxis_showbackground%3DTrue%0A%20%20%20%20)%0A%0A%20%20%20%20mo.as_html(fig3d).center()%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(input_ball_a%2C%20input_ball_b%2C%20mo)%3A%0A%20%20%20%20mo.hstack(%5Binput_ball_a%2C%20input_ball_b%5D)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(Ball%2C%20GAMMA%2C%20KM2M%2C%20PI%2C%20SI2MG%2C%20math)%3A%0A%20%20%20%20def%20bola(ball%3A%20Ball%2C%20xp%2C%20yp%2C%20zp)%3A%0A%20%20%20%20%20%20%20%20rx%20%3D%20xp%20-%20ball.x%0A%20%20%20%20%20%20%20%20ry%20%3D%20yp%20-%20ball.y%0A%20%20%20%20%20%20%20%20rz%20%3D%20zp%20-%20ball.z%0A%0A%20%20%20%20%20%20%20%20%23%20menghitung%20jarak%20menggunakan%20persamaan%20Pythagoras%0A%20%20%20%20%20%20%20%20r%20%3D%20math.sqrt(math.pow(rx%2C2)%2Bmath.pow(ry%2C2)%2Bmath.pow(rz%2C2))%0A%0A%20%20%20%20%20%20%20%20%23%20apabila%20jarak%3D0%20maka%20dibuat%20fungsi%20error%0A%20%20%20%20%20%20%20%20if%20r%20%3D%3D%200%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20raise%20ValueError(%22Observer%20%60observer%60%20and%20Object%20%60ball%60%20should%20be%20in%20the%20same%20point%22)%0A%0A%20%20%20%20%20%20%20%20%23jarak%20pangkat%203%0A%20%20%20%20%20%20%20%20r3%20%3D%20math.pow(r%2C3)%0A%0A%20%20%20%20%20%20%20%20%23%20menghitung%20total%20mass%0A%20%20%20%20%20%20%20%20tmass%20%3D%204.0%20*%20PI%20*%20ball.rho%20*%20(math.pow(ball.r%2C2))%2F3%0A%0A%20%20%20%20%20%20%20%20%23%20menghitung%20nilai%20percepatan%20gravitasi%20pada%20masing-masing%20komponen%0A%20%20%20%20%20%20%20%20gx%20%3D%20-GAMMA%20*%20tmass%20*%20rx%2Fr3%0A%20%20%20%20%20%20%20%20gy%20%3D%20-GAMMA%20*%20tmass%20*%20ry%2Fr3%0A%20%20%20%20%20%20%20%20gz%20%3D%20-GAMMA%20*%20tmass%20*%20rz%2Fr3%0A%0A%20%20%20%20%20%20%20%20%23%20konversi%20satuan%0A%20%20%20%20%20%20%20%20gx%20%3D%20gx%20*%20SI2MG%20*%20KM2M%0A%20%20%20%20%20%20%20%20gy%20%3D%20gy%20*%20SI2MG%20*%20KM2M%0A%20%20%20%20%20%20%20%20gz%20%3D%20gz%20*%20SI2MG%20*%20KM2M%0A%0A%20%20%20%20%20%20%20%20return%20gx%2C%20gy%2C%20gz%0A%20%20%20%20return%20(bola%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(ar)%3A%0A%20%20%20%20def%20linspace(lower%2C%20upper%2C%20n)%20-%3E%20ar.array%3A%0A%20%20%20%20%20%20%20%20_step%20%3D%20(upper-lower)%2Fn%0A%20%20%20%20%20%20%20%20_out%20%3D%20ar.array('d'%2C%20%5Blower%2B(i*_step)%20for%20i%20in%20range(n)%5D)%0A%20%20%20%20%20%20%20%20return%20_out%0A%0A%20%20%20%20def%20meshgrid(*xi%3A%20list%7Car.array)%20-%3E%20list%5Bar.array%5D%3A%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20Reimplementation%20of%20numpy.meshgrid%20using%20Python%20built-in%20arrays%0A%0A%20%20%20%20%20%20%20%20Parameters%3A%0A%20%20%20%20%20%20%20%20*xi%20%3A%20sequence%20of%201-D%20arrays%20(lists)%0A%20%20%20%20%20%20%20%20%20%20%20%20Representing%20the%20coordinates%20of%20a%20grid.%0A%0A%20%20%20%20%20%20%20%20Returns%3A%0A%20%20%20%20%20%20%20%20list%20of%20array.array%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20if%20not%20xi%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20%5B%5D%0A%0A%20%20%20%20%20%20%20%20xi%20%3D%20%5Blist(x)%20for%20x%20in%20xi%5D%0A%20%20%20%20%20%20%20%20ndim%20%3D%20len(xi)%0A%0A%20%20%20%20%20%20%20%20%23%20Initialize%20output%0A%20%20%20%20%20%20%20%20output%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20for%20i%2C%20x%20in%20enumerate(xi)%3A%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Create%20the%20broadcasted%20array%0A%20%20%20%20%20%20%20%20%20%20%20%20broadcasted%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20ndim%20%3D%3D%201%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20broadcasted%20%3D%20%5Bar.array('d'%2C%20%5Bval%5D)%20for%20val%20in%20x%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20i%20%3D%3D%200%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%23%20For%20x%20in%20xy%20mode%2C%20repeat%20rows%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20for%20val%20in%20x%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20broadcasted.append(ar.array('d'%2C%20%5Bval%5D*len(xi%5B1%5D)))%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20elif%20i%20%3D%3D%201%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%23%20For%20y%20in%20xy%20mode%2C%20repeat%20columns%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20broadcasted%20%3D%20%5Bar.array('d'%2C%20%5Bval%20for%20val%20in%20x%5D)%20for%20_%20in%20xi%5B0%5D%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20pass%0A%20%20%20%20%20%20%20%20%20%20%20%20output.append(broadcasted)%0A%0A%20%20%20%20%20%20%20%20return%20output%0A%0A%20%20%20%20def%20zeros(dim%3A%20tuple%5Bint%2C%20int%5D)%20-%3E%20list%5Bar.array%5D%3A%0A%20%20%20%20%20%20%20%20x%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20(nx%2C%20ny)%20%3D%20dim%0A%20%20%20%20%20%20%20%20for%20_%20in%20range(nx)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20x.append(ar.array('d'%2C%20%5B0.0%5D*ny))%0A%20%20%20%20%20%20%20%20return%20x%0A%20%20%20%20return%20linspace%2C%20meshgrid%2C%20zeros%0A%0A%0Aif%20__name__%20%3D%3D%20%22__main__%22%3A%0A%20%20%20%20app.run()%0A