peach

import%20marimo%0A%0A__generated_with%20%3D%20%220.23.16%22%0Aapp%20%3D%20marimo.App(width%3D%22full%22)%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20import%20marimo%20as%20mo%0A%20%20%20%20import%20matplotlib.pyplot%20as%20plt%0A%20%20%20%20import%20numpy%20as%20np%0A%0A%20%20%20%20import%20marimo_sandbox.jp_font%20as%20jp%0A%20%20%20%20from%20marimo_sandbox.data%20import%20load_titanic%0A%0A%20%20%20%20jp.setup_japanese_font()%0A%20%20%20%20return%20load_titanic%2C%20mo%2C%20np%2C%20plt%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23PEACH.py%EF%BC%88%E8%87%AA%E7%94%B1%E5%B8%B3%EF%BC%89%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20%E8%87%AA%E7%94%B1%E5%B8%B3%20%E2%80%94%20%E6%95%B0%E5%AD%A6%E3%81%AE%E7%BE%8E%E3%81%97%E3%81%95%E3%82%92%E5%91%B3%E3%82%8F%E3%81%86%0A%0A%20%20%20%20%E3%82%B5%E3%82%A4%E3%82%B3%E3%83%AD%E3%82%84%E5%AE%9D%E3%81%8F%E3%81%98%E3%81%AE%E9%81%8A%E3%81%B3%E3%82%92%E5%8D%92%E6%A5%AD%E3%81%97%E3%81%A6%E3%80%81%E7%A2%BA%E7%8E%87%E3%83%BB%E3%82%AB%E3%82%AA%E3%82%B9%E3%83%BB%E3%83%95%E3%83%A9%E3%82%AF%E3%82%BF%E3%83%AB%E3%83%BB%E6%95%B0%E8%AB%96%E3%81%8C%E6%8F%8F%E3%81%8F%0A%20%20%20%20%E3%80%8C%E8%87%AA%E7%84%B6%E7%95%8C%E3%81%AE%E8%83%8C%E5%BE%8C%E3%81%AB%E3%81%82%E3%82%8B%E7%A7%A9%E5%BA%8F%E3%80%8D%E3%82%92%E8%A6%B3%E5%AF%9F%E3%81%99%E3%82%8B%E3%83%8E%E3%83%BC%E3%83%88%E3%81%A7%E3%81%99%E3%80%82%0A%20%20%20%20%E3%81%A9%E3%82%8C%E3%82%82%E6%95%B0%E5%8D%83%E8%A1%8C%E3%81%AE%E8%A8%88%E7%AE%97%E3%81%8C%E4%B8%80%E7%9E%AC%E3%81%A7%E7%B5%82%E3%82%8F%E3%82%8B%E3%80%81%E3%82%B3%E3%83%B3%E3%83%94%E3%83%A5%E3%83%BC%E3%82%BF%E3%81%A0%E3%81%8B%E3%82%89%E3%81%93%E3%81%9D%E8%A6%8B%E3%81%88%E3%82%8B%E7%BE%8E%E3%81%97%E3%81%95%E3%81%A7%E3%81%99%E3%80%82%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20---%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%23%201.%20%E3%83%A2%E3%83%B3%E3%83%86%E3%82%AB%E3%83%AB%E3%83%AD%E6%B3%95%E3%81%A7%E5%86%86%E5%91%A8%E7%8E%87%E3%82%92%E6%8E%A8%E5%AE%9A%E3%81%99%E3%82%8B%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20n_samples%20%3D%20mo.ui.slider(%0A%20%20%20%20%20%20%20%201_000%2C%20200_000%2C%20step%3D1_000%2C%20value%3D50_000%2C%20label%3D%22%E3%82%B5%E3%83%B3%E3%83%97%E3%83%AB%E7%82%B9%E6%95%B0%22%0A%20%20%20%20)%0A%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20r%22%22%22%0A%20%20%20%20%20%20%20%20%E5%8D%98%E4%BD%8D%E6%AD%A3%E6%96%B9%E5%BD%A2%E3%81%AB%E3%83%A9%E3%83%B3%E3%83%80%E3%83%A0%E3%81%AB%E7%82%B9%E3%82%92%E6%89%93%E3%81%A1%E3%80%81%E5%8D%98%E4%BD%8D%E5%86%86(%E7%AC%AC1%E8%B1%A1%E9%99%90)%E3%81%AB%E5%85%A5%E3%82%8B%E5%89%B2%E5%90%88%E3%82%92%E6%95%B0%E3%81%88%E3%82%8C%E3%81%B0%0A%20%20%20%20%20%20%20%20%E5%86%86%E5%91%A8%E7%8E%87%E3%81%8C%E6%8E%A8%E5%AE%9A%E3%81%A7%E3%81%8D%E3%81%BE%E3%81%99%E3%80%82%E4%B9%B1%E6%95%B0%E3%81%A8%E3%81%84%E3%81%86%E3%80%8C%E5%81%B6%E7%84%B6%E3%80%8D%E3%82%92%E5%A4%A7%E9%87%8F%E3%81%AB%E9%87%8D%E3%81%AD%E3%82%8B%E3%81%A8%E3%80%81%E7%A2%BA%E7%8E%87%E3%81%8C%0A%20%20%20%20%20%20%20%20%E7%9C%9F%E7%90%86%E3%81%B8%E3%81%A8%E5%8F%8E%E6%9D%9F%E3%81%97%E3%81%A6%E3%81%84%E3%81%8F%E6%A7%98%E5%AD%90%E3%82%92%E5%AE%9F%E6%84%9F%E3%81%A7%E3%81%8D%E3%81%BE%E3%81%99%E3%80%82%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20)%0A%20%20%20%20n_samples%0A%20%20%20%20return%20(n_samples%2C)%0A%0A%0A%40app.cell%0Adef%20_(mo%2C%20n_samples%2C%20np%2C%20plt)%3A%0A%20%20%20%20rng%20%3D%20np.random.default_rng(2025)%0A%20%20%20%20pts%20%3D%20rng.random((n_samples.value%2C%202))%0A%20%20%20%20inside%20%3D%20pts%5B%3A%2C%200%5D%20**%202%20%2B%20pts%5B%3A%2C%201%5D%20**%202%20%3C%3D%201.0%0A%20%20%20%20pi_est%20%3D%20float(4%20*%20inside.mean())%0A%20%20%20%20error%20%3D%20abs(np.pi%20-%20pi_est)%0A%0A%20%20%20%20n%20%3D%20np.arange(1%2C%20n_samples.value%20%2B%201)%0A%20%20%20%20cum_pi%20%3D%204%20*%20np.cumsum(inside)%20%2F%20n%0A%0A%20%20%20%20fig_mc%2C%20axes_mc%20%3D%20plt.subplots(1%2C%202%2C%20figsize%3D(12%2C%204.5))%0A%20%20%20%20axes_mc%5B0%5D.scatter(%0A%20%20%20%20%20%20%20%20pts%5B~inside%2C%200%5D%2C%20pts%5B~inside%2C%201%5D%2C%20s%3D1%2C%20c%3D%22tomato%22%2C%20alpha%3D0.4%2C%20label%3D%22%E5%86%86%E3%81%AE%E5%A4%96%E5%81%B4%22%0A%20%20%20%20)%0A%20%20%20%20axes_mc%5B0%5D.scatter(%0A%20%20%20%20%20%20%20%20pts%5Binside%2C%200%5D%2C%20pts%5Binside%2C%201%5D%2C%20s%3D1%2C%20c%3D%22steelblue%22%2C%20alpha%3D0.4%2C%20label%3D%22%E5%86%86%E3%81%AE%E5%86%85%E5%81%B4%22%0A%20%20%20%20)%0A%20%20%20%20theta%20%3D%20np.linspace(0%2C%20np.pi%20%2F%202%2C%20200)%0A%20%20%20%20axes_mc%5B0%5D.plot(np.cos(theta)%2C%20np.sin(theta)%2C%20%22k-%22%2C%20linewidth%3D1.5)%0A%20%20%20%20axes_mc%5B0%5D.set_aspect(%22equal%22)%0A%20%20%20%20axes_mc%5B0%5D.set_title(f%22%E5%8D%98%E4%BD%8D%E5%86%86%E3%81%B8%E3%81%AE%E7%82%B9%E3%81%AE%E6%89%93%E7%82%B9%20%E2%80%94%20%E6%8E%A8%E5%AE%9A%E5%80%A4%20%CF%80%20%E2%89%88%20%7Bpi_est%3A.5f%7D%22)%0A%20%20%20%20axes_mc%5B0%5D.legend(loc%3D%22upper%20left%22)%0A%20%20%20%20axes_mc%5B1%5D.plot(n%2C%20cum_pi%2C%20lw%3D1%2C%20color%3D%22steelblue%22)%0A%20%20%20%20axes_mc%5B1%5D.axhline(np.pi%2C%20color%3D%22crimson%22%2C%20ls%3D%22--%22%2C%20label%3D%22%E7%9C%9F%E3%81%AE%20%CF%80%22)%0A%20%20%20%20axes_mc%5B1%5D.set_xscale(%22log%22)%0A%20%20%20%20axes_mc%5B1%5D.set_xlabel(%22%E3%82%B5%E3%83%B3%E3%83%97%E3%83%AB%E6%95%B0%20(%E5%AF%BE%E6%95%B0)%22)%0A%20%20%20%20axes_mc%5B1%5D.set_ylabel(%22%CF%80%20%E3%81%AE%E6%8E%A8%E5%AE%9A%E5%80%A4%22)%0A%20%20%20%20axes_mc%5B1%5D.set_title(%22%E3%82%B5%E3%83%B3%E3%83%97%E3%83%AB%E6%95%B0%E3%81%A8%20%CF%80%20%E6%8E%A8%E5%AE%9A%E5%80%A4%E3%81%AE%E5%8F%8E%E6%9D%9F%22)%0A%20%20%20%20axes_mc%5B1%5D.legend()%0A%20%20%20%20fig_mc.tight_layout()%0A%0A%20%20%20%20mo.md(f%22**%E6%8E%A8%E5%AE%9A%E5%80%A4**%3A%20%CF%80%20%E2%89%88%20**%7Bpi_est%3A.6f%7D**%EF%BC%88%E7%9C%9F%E3%81%AE%E5%80%A4%E3%81%A8%E3%81%AE%E8%AA%A4%E5%B7%AE%20%7Berror%3A.2e%7D%EF%BC%89%22)%0A%20%20%20%20mo.mpl.interactive(fig_mc)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20---%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%23%202.%20%E3%83%AD%E3%82%B8%E3%82%B9%E3%83%86%E3%82%A3%E3%83%83%E3%82%AF%E5%86%99%E5%83%8F%E3%81%AE%E5%88%86%E5%B2%90%E5%9B%B3%EF%BC%88%E3%82%AB%E3%82%AA%E3%82%B9%EF%BC%89%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(mo%2C%20np%2C%20plt)%3A%0A%20%20%20%20r_values%20%3D%20np.linspace(2.5%2C%204.0%2C%202_000)%0A%20%20%20%20orbit%20%3D%20np.full_like(r_values%2C%200.5)%0A%20%20%20%20for%20_%20in%20range(1_000)%3A%0A%20%20%20%20%20%20%20%20orbit%20%3D%20r_values%20*%20orbit%20*%20(1%20-%20orbit)%0A%0A%20%20%20%20fig_chaos%2C%20ax_chaos%20%3D%20plt.subplots(figsize%3D(12%2C%205))%0A%20%20%20%20for%20_%20in%20range(250)%3A%0A%20%20%20%20%20%20%20%20orbit%20%3D%20r_values%20*%20orbit%20*%20(1%20-%20orbit)%0A%20%20%20%20%20%20%20%20ax_chaos.plot(r_values%2C%20orbit%2C%20%22%2Ck%22%2C%20alpha%3D0.25)%0A%20%20%20%20ax_chaos.set_xlabel(%22%E3%83%91%E3%83%A9%E3%83%A1%E3%83%BC%E3%82%BF%20r%22)%0A%20%20%20%20ax_chaos.set_ylabel(%22%E9%95%B7%E6%9C%9F%E7%9A%84%E3%81%AA%E6%8C%AF%E3%82%8B%E8%88%9E%E3%81%84%20x%22)%0A%20%20%20%20ax_chaos.set_title(%22%E3%83%AD%E3%82%B8%E3%82%B9%E3%83%86%E3%82%A3%E3%83%83%E3%82%AF%E5%86%99%E5%83%8F%20%24x_%7Bn%2B1%7D%20%3D%20r%20%5C%5C%2C%20x_n%20(1%20-%20x_n)%24%20%E3%81%AE%E5%88%86%E5%B2%90%E5%9B%B3%22)%0A%20%20%20%20fig_chaos.tight_layout()%0A%0A%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20r%22%22%22%0A%20%20%20%20%20%20%20%20%E3%81%9F%E3%81%A3%E3%81%9F%201%20%E6%9C%AC%E3%81%AE%E7%B0%A1%E5%8D%98%E3%81%AA%E6%BC%B8%E5%8C%96%E5%BC%8F%20%24x_%7Bn%2B1%7D%20%3D%20r%20x_n%20(1%20-%20x_n)%24%20%E3%81%8C%E3%80%81%24r%24%20%E3%81%AE%E5%80%A4%E3%81%AB%E3%82%88%E3%81%A3%E3%81%A6%0A%20%20%20%20%20%20%20%20%E4%B8%80%E7%82%B9%E3%81%AB%E8%90%BD%E3%81%A1%E7%9D%80%E3%81%8D%EF%BC%88%E5%8F%8E%E6%9D%9F%EF%BC%89%E2%86%92%20%E6%8C%AF%E5%8B%95%20%E2%86%92%20%E3%82%AB%E3%82%AA%E3%82%B9%E3%81%B8%E3%81%A8%E5%8A%87%E7%9A%84%E3%81%AB%E5%A4%89%E5%8C%96%E3%81%97%E3%81%BE%E3%81%99%E3%80%82%E5%8D%98%E7%B4%94%E3%81%AA%E8%A6%8F%E5%89%87%E3%81%8C%0A%20%20%20%20%20%20%20%20%E4%BA%88%E6%B8%AC%E4%B8%8D%E8%83%BD%E3%81%AA%E8%A4%87%E9%9B%91%E3%81%95%E3%82%92%E7%94%9F%E3%82%80%E3%80%8C%E6%B1%BA%E5%AE%9A%E8%AB%96%E7%9A%84%E3%82%AB%E3%82%AA%E3%82%B9%E3%80%8D%E3%81%AE%E4%BB%A3%E8%A1%A8%E4%BE%8B%E3%81%A7%E3%81%99%E3%80%82%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20)%0A%20%20%20%20mo.mpl.interactive(fig_chaos)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20---%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%23%203.%20%E3%83%9E%E3%83%B3%E3%83%87%E3%83%AB%E3%83%96%E3%83%AD%E9%9B%86%E5%90%88%E3%82%92%E6%8E%A2%E7%B4%A2%E3%81%99%E3%82%8B%EF%BC%88%E3%83%95%E3%83%A9%E3%82%AF%E3%82%BF%E3%83%AB%EF%BC%89%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20center_x%20%3D%20mo.ui.slider(-1.0%2C%200.5%2C%20step%3D0.01%2C%20value%3D-0.5%2C%20label%3D%22%E4%B8%AD%E5%BF%83%20x%22)%0A%20%20%20%20center_y%20%3D%20mo.ui.slider(-0.5%2C%200.5%2C%20step%3D0.01%2C%20value%3D0.0%2C%20label%3D%22%E4%B8%AD%E5%BF%83%20y%22)%0A%20%20%20%20zoom%20%3D%20mo.ui.slider(0%2C%206%2C%20step%3D1%2C%20value%3D0%2C%20label%3D%22%E3%82%BA%E3%83%BC%E3%83%A0%E6%AE%B5%E9%9A%8E%22)%0A%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20r%22%22%22%0A%20%20%20%20%20%20%20%20%E8%A4%87%E7%B4%A0%E5%B9%B3%E9%9D%A2%E4%B8%8A%E3%81%AE%E5%90%84%E7%82%B9%20%24c%24%20%E3%81%AB%E5%AF%BE%E3%81%97%E3%80%81%24z_%7Bn%2B1%7D%20%3D%20z_n%5E2%20%2B%20c%24%20%E3%81%8C%E7%99%BA%E6%95%A3%E3%81%99%E3%82%8B%E3%81%8B%E3%81%A9%E3%81%86%E3%81%8B%E3%82%92%0A%20%20%20%20%20%20%20%20%E5%88%A4%E5%AE%9A%E3%81%97%E3%81%A6%E8%89%B2%E4%BB%98%E3%81%91%E3%81%97%E3%81%9F%E5%9B%B3%E3%81%A7%E3%81%99%E3%80%82%E5%A2%83%E7%95%8C%E3%81%AF%E3%81%84%E3%81%8F%E3%82%89%E6%8B%A1%E5%A4%A7%E3%81%97%E3%81%A6%E3%82%82%E5%90%8C%E3%81%98%E3%82%88%E3%81%86%E3%81%AA%E8%A4%87%E9%9B%91%E3%81%95%E3%81%8C%E7%8F%BE%E3%82%8C%E3%82%8B%0A%20%20%20%20%20%20%20%20%E3%80%8C%E8%87%AA%E5%B7%B1%E7%9B%B8%E4%BC%BC%E3%80%8D%E3%81%AE%E6%A7%8B%E9%80%A0%E3%82%92%E6%8C%81%E3%81%A1%E3%81%BE%E3%81%99%E3%80%82%E3%81%8A%E3%81%99%E3%81%99%E3%82%81%E3%81%AF%E3%83%9E%E3%83%89%E3%83%A9%E3%83%BC%E3%81%AE%E6%B5%B7%E9%A6%AC%E8%B0%B7%0A%20%20%20%20%20%20%20%20%24c%20%5C%5Capprox%20(-0.7436%2C%200.1319)%24%20%E3%81%A7%E3%81%99%E3%80%82%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20)%0A%20%20%20%20center_x%2C%20center_y%2C%20zoom%0A%20%20%20%20return%20center_x%2C%20center_y%2C%20zoom%0A%0A%0A%40app.cell%0Adef%20_(center_x%2C%20center_y%2C%20mo%2C%20np%2C%20plt%2C%20zoom)%3A%0A%20%20%20%20width%20%3D%202.5%20%2F%20(1.8**zoom.value)%0A%20%20%20%20n_pts%20%3D%20500%0A%20%20%20%20max_iter%20%3D%20120%0A%0A%20%20%20%20grid_x%20%3D%20np.linspace(center_x.value%20-%20width%2C%20center_x.value%20%2B%20width%2C%20n_pts)%0A%20%20%20%20grid_y%20%3D%20np.linspace(center_y.value%20-%20width%2C%20center_y.value%20%2B%20width%2C%20n_pts)%0A%20%20%20%20c%20%3D%20grid_x%5BNone%2C%20%3A%5D%20%2B%201j%20*%20grid_y%5B%3A%2C%20None%5D%0A%20%20%20%20z%20%3D%20np.zeros_like(c%2C%20dtype%3Dnp.complex128)%0A%20%20%20%20counts%20%3D%20np.zeros(c.shape%2C%20dtype%3Dnp.uint16)%0A%20%20%20%20alive%20%3D%20np.ones(c.shape%2C%20dtype%3Dbool)%0A%0A%20%20%20%20for%20i%20in%20range(max_iter)%3A%0A%20%20%20%20%20%20%20%20z%5Balive%5D%20%3D%20z%5Balive%5D%20**%202%20%2B%20c%5Balive%5D%0A%20%20%20%20%20%20%20%20escaped%20%3D%20np.abs(z)%20%3E%202%0A%20%20%20%20%20%20%20%20counts%5Bescaped%20%26%20alive%5D%20%3D%20i%0A%20%20%20%20%20%20%20%20alive%20%26%3D%20~escaped%0A%20%20%20%20%20%20%20%20if%20not%20alive.any()%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20break%0A%20%20%20%20counts%5Balive%5D%20%3D%20max_iter%0A%0A%20%20%20%20fig_mandel%2C%20ax_mandel%20%3D%20plt.subplots(figsize%3D(9%2C%209))%0A%20%20%20%20image%20%3D%20ax_mandel.imshow(%0A%20%20%20%20%20%20%20%20counts%2C%0A%20%20%20%20%20%20%20%20extent%3D(center_x.value%20-%20width%2C%20center_x.value%20%2B%20width%2C%20center_y.value%20-%20width%2C%20center_y.value%20%2B%20width)%2C%0A%20%20%20%20%20%20%20%20origin%3D%22lower%22%2C%0A%20%20%20%20%20%20%20%20cmap%3D%22magma%22%2C%0A%20%20%20%20%20%20%20%20interpolation%3D%22bilinear%22%2C%0A%20%20%20%20)%0A%20%20%20%20ax_mandel.set_title(f%22%E3%83%9E%E3%83%B3%E3%83%87%E3%83%AB%E3%83%96%E3%83%AD%E9%9B%86%E5%90%88%EF%BC%88%E4%B8%AD%E5%BF%83%20(%7Bcenter_x.value%7D%2C%20%7Bcenter_y.value%7D)%EF%BC%8C%E5%B9%85%20%7B2%20*%20width%3A.4f%7D%EF%BC%89%22)%0A%20%20%20%20fig_mandel.colorbar(image%2C%20ax%3Dax_mandel%2C%20label%3D%22%E7%99%BA%E6%95%A3%E3%81%BE%E3%81%A7%E3%81%AE%E5%8F%8D%E5%BE%A9%E5%9B%9E%E6%95%B0%22)%0A%20%20%20%20fig_mandel.tight_layout()%0A%20%20%20%20mo.mpl.interactive(fig_mandel)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20---%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%23%204.%20%E3%83%99%E3%83%B3%E3%83%95%E3%82%A9%E3%83%BC%E3%83%89%E3%81%AE%E6%B3%95%E5%89%87%E3%82%92%E5%AE%9F%E3%83%87%E3%83%BC%E3%82%BF%E3%81%A7%E6%A4%9C%E8%A8%BC%E3%81%99%E3%82%8B%EF%BC%88%E6%95%B0%E8%AB%96%EF%BC%89%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(load_titanic)%3A%0A%20%20%20%20df%20%3D%20load_titanic()%0A%20%20%20%20fares%20%3D%20df%5B%22fare%22%5D.dropna()%0A%20%20%20%20fares%20%3D%20fares%5Bfares%20%3E%200%5D%0A%20%20%20%20first_digits%20%3D%20fares.astype(str).str%5B0%5D.astype(int)%0A%20%20%20%20return%20(first_digits%2C)%0A%0A%0A%40app.cell%0Adef%20_(first_digits%2C%20mo%2C%20np%2C%20plt)%3A%0A%20%20%20%20observed%20%3D%20first_digits.value_counts(normalize%3DTrue).sort_index()%0A%20%20%20%20digits%20%3D%20np.arange(1%2C%2010)%0A%20%20%20%20expected%20%3D%20np.log10(1%20%2B%201%20%2F%20digits)%0A%0A%20%20%20%20fig_benford%2C%20ax_benford%20%3D%20plt.subplots(figsize%3D(9%2C%204.5))%0A%20%20%20%20ax_benford.bar(observed.index%2C%20observed.values%2C%20width%3D0.6%2C%20color%3D%22steelblue%22%2C%20alpha%3D0.75%2C%20label%3D%22%E3%82%BF%E3%82%A4%E3%82%BF%E3%83%8B%E3%83%83%E3%82%AF%E5%8F%B7%E3%81%AE%E9%81%8B%E8%B3%83%22)%0A%20%20%20%20ax_benford.plot(digits%2C%20expected%2C%20%22o-%22%2C%20color%3D%22crimson%22%2C%20label%3Dr%22%E3%83%99%E3%83%B3%E3%83%95%E3%82%A9%E3%83%BC%E3%83%89%E3%81%AE%E6%B3%95%E5%89%87%20%24P(d)%20%3D%20%5Clog_%7B10%7D(1%20%2B%201%2Fd)%24%22)%0A%20%20%20%20ax_benford.set_xticks(digits)%0A%20%20%20%20ax_benford.set_xlabel(%22%E5%85%88%E9%A0%AD%E3%81%AE%E6%95%B0%E5%AD%97%20d%22)%0A%20%20%20%20ax_benford.set_ylabel(%22%E5%87%BA%E7%8F%BE%E9%A0%BB%E5%BA%A6%22)%0A%20%20%20%20ax_benford.set_title(%22%E8%87%AA%E7%84%B6%E7%95%8C%E3%81%AE%E6%95%B0%E5%AD%97%E3%81%AF%201%20%E3%81%8B%E3%82%89%E5%A7%8B%E3%81%BE%E3%82%8B%E3%82%82%E3%81%AE%E3%81%8C%E4%B8%80%E7%95%AA%E5%A4%9A%E3%81%84%22)%0A%20%20%20%20ax_benford.legend()%0A%20%20%20%20fig_benford.tight_layout()%0A%0A%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20r%22%22%22%0A%20%20%20%20%20%20%20%20%E6%A0%AA%E4%BE%A1%E3%82%84%E4%BA%BA%E5%8F%A3%E3%80%81%E6%B2%B3%E5%B7%9D%E3%81%AE%E9%95%B7%E3%81%95%E3%81%AE%E3%82%88%E3%81%86%E3%81%AA%E3%80%8C%E8%87%AA%E7%84%B6%E3%81%AB%E7%94%9F%E3%81%BE%E3%82%8C%E3%81%9F%E6%95%B0%E5%80%A4%E3%80%8D%E3%81%A7%E3%81%AF%E3%80%81%E5%85%88%E9%A0%AD%E3%81%AE%E6%95%B0%E5%AD%97%E3%81%8C%0A%20%20%20%20%20%20%20%20%241%24%20%E3%81%A7%E3%81%82%E3%82%8B%E7%A2%BA%E7%8E%87%E3%81%8C%E7%B4%84%2030%25%20%E3%81%A8%E6%9C%80%E3%82%82%E9%AB%98%E3%81%8F%E3%80%81%249%24%20%E3%81%AB%E3%81%AA%E3%82%8B%E3%81%BB%E3%81%A9%E4%BD%8E%E3%81%8F%E3%81%AA%E3%82%8A%E3%81%BE%E3%81%99%E3%80%82%E3%82%BF%E3%82%A4%E3%82%BF%E3%83%8B%E3%83%83%E3%82%AF%E3%81%AE%0A%20%20%20%20%20%20%20%20%E9%81%8B%E8%B3%83%E3%83%87%E3%83%BC%E3%82%BF%E3%82%82%E3%81%93%E3%81%AE%E6%B3%95%E5%89%87%E3%81%AB%E3%81%8D%E3%82%8C%E3%81%84%E3%81%AB%E5%BE%93%E3%81%86%E3%81%8B%E3%82%92%E7%A2%BA%E8%AA%8D%E3%81%97%E3%81%A6%E3%81%BF%E3%81%BE%E3%81%97%E3%82%87%E3%81%86%E3%80%82%E4%BC%9A%E8%A8%88%E4%B8%8D%E6%AD%A3%E3%81%AE%E6%A4%9C%E5%87%BA%E3%81%AB%E3%82%82%0A%20%20%20%20%20%20%20%20%E4%BD%BF%E3%82%8F%E3%82%8C%E3%82%8B%E6%9C%89%E5%90%8D%E3%81%AA%E6%95%B0%E8%AB%96%E3%81%AE%E5%BF%9C%E7%94%A8%E3%81%A7%E3%81%99%E3%80%82%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20)%0A%20%20%20%20mo.mpl.interactive(fig_benford)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20---%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%23%20%E3%81%BE%E3%81%A8%E3%82%81%0A%0A%20%20%20%20%E9%81%8A%E3%81%B3%E3%82%92%E5%8D%92%E6%A5%AD%E3%81%97%E3%81%A6%E3%82%82%E3%80%81%E7%A2%BA%E7%8E%87%E3%83%BB%E3%82%AB%E3%82%AA%E3%82%B9%E3%83%BB%E3%83%95%E3%83%A9%E3%82%AF%E3%82%BF%E3%83%AB%E3%83%BB%E6%95%B0%E8%AB%96%E3%81%AB%E3%81%AF%E3%81%84%E3%81%8F%E3%82%89%E3%81%A7%E3%82%82%E3%80%8C%E8%87%AA%E7%94%B1%E5%B8%B3%E3%80%8D%E3%81%AB%0A%20%20%20%20%E6%8F%8F%E3%81%8D%E3%81%9F%E3%81%8F%E3%81%AA%E3%82%8B%E8%A9%B1%E9%A1%8C%E3%81%8C%E3%81%82%E3%82%8A%E3%81%BE%E3%81%99%E3%80%82%E6%AC%A1%E3%81%AF%E3%81%93%E3%81%AE%E3%83%8E%E3%83%BC%E3%83%88%E3%81%AE%E5%90%84%E3%82%BB%E3%82%AF%E3%82%B7%E3%83%A7%E3%83%B3%E3%82%92%E6%94%B9%E9%80%A0%E3%81%97%E3%81%A6%E3%80%81%0A%20%20%20%20%E8%87%AA%E5%88%86%E3%81%A0%E3%81%91%E3%81%AE%E8%A6%B3%E5%AF%9F%E3%82%92%E5%8A%A0%E3%81%88%E3%81%A6%E3%81%BF%E3%81%A6%E3%81%8F%E3%81%A0%E3%81%95%E3%81%84%E3%80%82%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%20%22%22%22)%0A%20%20%20%20return%0A%0A%0Aif%20__name__%20%3D%3D%20%22__main__%22%3A%0A%20%20%20%20app.run()%0A