QRCodeDemo

import%20marimo%0A%0A__generated_with%20%3D%20%220.16.5%22%0Aapp%20%3D%20marimo.App(width%3D%22medium%22)%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20import%20numpy%20as%20np%0A%20%20%20%20import%20marimo%20as%20mo%0A%20%20%20%20import%20seaborn%20as%20sns%0A%20%20%20%20import%20matplotlib.pyplot%20as%20plt%0A%20%20%20%20from%20scipy.constants%20import%20pi%0A%20%20%20%20from%20scipy.stats%20import%20norm%2C%20cauchy%0A%0A%20%20%20%20EPSILON%20%3D%201e-3%0A%20%20%20%20sns.set_theme(style%3D%22white%22%2Ccontext%3D%22talk%22%2Cfont_scale%3D1.25%2Crc%3D%7B%22axes.edgecolor%22%3A%20%22black%22%2C%20%22grid.color%22%3A%20%22silver%22%7D)%0A%20%20%20%20plt.rcParams%5B%22figure.figsize%22%5D%20%3D%20(10%2C%2010)%0A%20%20%20%20plt.rcParams%5B%22figure.dpi%22%5D%20%3D%20140%0A%20%20%20%20plt.rcParams%5B%22lines.linewidth%22%5D%20%3D%201.5%0A%20%20%20%20plt.rcParams%5B%22axes.titlesize%22%5D%20%3D%2026%0A%20%20%20%20plt.rcParams%5B%22axes.titlepad%22%5D%20%3D%2015%0A%20%20%20%20plt.rcParams%5B%22xtick.labelsize%22%5D%20%3D%2018%0A%20%20%20%20plt.rcParams%5B%22ytick.labelsize%22%5D%20%3D%2018%0A%20%20%20%20plt.rcParams%5B%22xtick.bottom%22%5D%20%3D%20True%0A%20%20%20%20plt.rcParams%5B%22ytick.left%22%5D%20%3D%20True%0A%20%20%20%20plt.rcParams%5B%22xtick.minor.visible%22%5D%20%3D%20True%0A%20%20%20%20plt.rcParams%5B%22ytick.minor.visible%22%5D%20%3D%20True%0A%20%20%20%20plt.rcParams%5B%22xtick.minor.size%22%5D%20%3D%206%0A%20%20%20%20plt.rcParams%5B%22ytick.minor.size%22%5D%20%3D%206%0A%20%20%20%20plt.rcParams%5B%22xtick.color%22%5D%20%3D%20%22black%22%0A%20%20%20%20plt.rcParams%5B%22ytick.color%22%5D%20%3D%20%22black%22%0A%20%20%20%20return%20EPSILON%2C%20cauchy%2C%20mo%2C%20norm%2C%20np%2C%20pi%2C%20plt%0A%0A%0A%40app.cell%0Adef%20_(EPSILON%2C%20cauchy%2C%20norm%2C%20np%2C%20pi)%3A%0A%20%20%20%20def%20n_kc(x%2C%20a%2C%20b%2C%20gamma%2C%20mu)%3A%0A%20%20%20%20%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20norm.pdf(cauchy.ppf((1%20-%20np.power(1%20-%20np.power(x%2C%20a)%2C%20b))%2C%20scale%3Dgamma)%2C%20loc%3Dmu)%0A%20%20%20%20%20%20%20%20%20%20%20%20*%20(gamma%20*%20(1%20%2F%20np.power(np.cos(pi%20*%20((1%20-%20np.power(1%20-%20np.power(x%2C%20a)%2C%20b))%20-%200.5))%2C%202))%20*%20pi)%0A%20%20%20%20%20%20%20%20%20%20%20%20*%20(a%20*%20b%20*%20np.power(x%2C%20a%20-%201)%20*%20np.power(np.maximum((1.0%20-%20np.power(x%2C%20a))%2C%20EPSILON)%2C%20b%20-%201))%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20return%20(n_kc%2C)%0A%0A%0A%40app.cell%0Adef%20_(EPSILON%2C%20mo%2C%20n_kc%2C%20np%2C%20plt)%3A%0A%20%20%20%20def%20n_kc_fig(a%2C%20b%2C%20gamma%2C%20mu)%3A%0A%20%20%20%20%20%20%20%20x%20%3D%20np.linspace(EPSILON%2C%201%20-%20EPSILON%2C%20500)%0A%20%20%20%20%20%20%20%20fig%2C%20ax%20%3D%20plt.subplots()%0A%20%20%20%20%20%20%20%20ax.plot(x%2C%20n_kc(x%2C%20a%2C%20b%2C%20gamma%2C%20mu))%0A%20%20%20%20%20%20%20%20ax.set_xlim(0%2C%201)%0A%20%20%20%20%20%20%20%20ax.set_ylim(0%2C%205.6)%0A%20%20%20%20%20%20%20%20ax.set_xlabel(xlabel%3Dr%22%24X%24%20Value%22)%0A%20%20%20%20%20%20%20%20ax.set_ylabel(%22Density%22)%0A%20%20%20%20%20%20%20%20ax.set_title(r'%24%5C%2CN%5C!-%5C!K%5C%7BC%5C%7D%5C%2C%24%20PDF')%0A%20%20%20%20%20%20%20%20ax.grid()%0A%20%20%20%20%20%20%20%20return%20fig%0A%0A%20%20%20%20a%20%3D%20mo.ui.slider(start%3D0.5%2C%20stop%3D2.5%2C%20step%3D0.10%2C%20value%3D1%2C%20label%3D'%24a%24'%2C%20show_value%3DTrue)%0A%20%20%20%20b%20%3D%20mo.ui.slider(start%3D0.5%2C%20stop%3D2.5%2C%20step%3D0.10%2C%20value%3D1%2C%20label%3D'%24b%24'%2C%20show_value%3DTrue)%0A%20%20%20%20gamma%20%3D%20mo.ui.slider(start%3D0.5%2C%20stop%3D4.5%2C%20step%3D0.10%2C%20value%3D2%2C%20label%3Dr'%24%5Cgamma%24'%2C%20show_value%3DTrue)%0A%20%20%20%20mu%20%3D%20mo.ui.slider(start%3D-1.5%2C%20stop%3D1.5%2C%20step%3D0.10%2C%20value%3D0%2C%20label%3Dr'%24%5Cmu%24'%2C%20show_value%3DTrue)%0A%20%20%20%20return%20a%2C%20b%2C%20gamma%2C%20mu%2C%20n_kc_fig%0A%0A%0A%40app.cell%0Adef%20_(a%2C%20b%2C%20gamma%2C%20mo%2C%20mu%2C%20n_kc_fig)%3A%0A%20%20%20%20mo.vstack(%5B%0A%20%20%20%20%20%20%20%20mo.md(%22%23%20Parameters%20%3A%22)%2C%0A%20%20%20%20%20%20%20%20mo.md(%22%23%22)%2C%0A%20%20%20%20%20%20%20%20mo.vstack(%5Ba%2C%20b%2C%20gamma%2C%20mu%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20justify%3D'center'%2C%20align%3D'center'%2C%20gap%3D0.5)%2C%0A%20%20%20%20%20%20%20%20mo.as_html(n_kc_fig(a.value%2C%20b.value%2C%20gamma.value%2C%20mu.value))%0A%20%20%20%20%20%20%20%20%5D)%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