SCS-SA

import%20marimo%0A%0A__generated_with%20%3D%20%220.11.23%22%0Aapp%20%3D%20marimo.App(width%3D%22medium%22)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(%22%22%22%23**SCS-SA**%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20p01%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%22Parameter%22%3A%20mo.md(%22Area%22)%2C%0A%20%20%20%20%20%20%20%20%22Value%22%3A%20mo.ui.text()%2C%0A%20%20%20%20%20%20%20%20%22Units%22%3A%20mo.md(%22km%3Csup%3E2%3C%2Fsup%3E%22)%2C%0A%20%20%20%20%7D%0A%0A%20%20%20%20p02%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%22Parameter%22%3A%20mo.md(%22Catchment%20Lag%20Time%22)%2C%0A%20%20%20%20%20%20%20%20%22Value%22%3A%20mo.ui.text()%2C%0A%20%20%20%20%20%20%20%20%22Units%22%3A%20mo.md(%22hours%22)%2C%0A%20%20%20%20%7D%0A%0A%20%20%20%20p03%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%22Parameter%22%3A%20mo.md(%22Rainfall%20Depth%22)%2C%0A%20%20%20%20%20%20%20%20%22Value%22%3A%20mo.ui.number()%2C%0A%20%20%20%20%20%20%20%20%22Units%22%3A%20mo.md(%22mm%22)%2C%0A%20%20%20%20%7D%0A%0A%20%20%20%20p04%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%22Parameter%22%3A%20mo.md(%22Curve%20Number%22)%2C%0A%20%20%20%20%20%20%20%20%22Value%22%3A%20mo.ui.text()%2C%0A%20%20%20%20%20%20%20%20%22Units%22%3A%20mo.md(%22%22)%2C%0A%20%20%20%20%7D%0A%0A%20%20%20%20p05%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%22Parameter%22%3A%20mo.md(%22Rainfall%20intensity%20distribution%20type%22)%2C%0A%20%20%20%20%20%20%20%20%22Value%22%3A%20mo.ui.dropdown(options%3D%5B%221%22%2C%20%222%22%2C%20%223%22%2C%20%224%22%5D)%2C%0A%20%20%20%20%20%20%20%20%22Units%22%3A%20mo.md(%22%22)%2C%0A%20%20%20%20%7D%0A%20%20%20%20return%20p01%2C%20p02%2C%20p03%2C%20p04%2C%20p05%0A%0A%0A%40app.cell%0Adef%20_(mo%2C%20p01%2C%20p02%2C%20p03%2C%20p04%2C%20p05)%3A%0A%20%20%20%20table%20%3D%20mo.ui.table(%0A%20%20%20%20%20%20%20%20data%3D%5Bp01%2C%20p02%2C%20p03%2C%20p04%2C%20p05%5D%2C%0A%20%20%20%20%20%20%20%20label%3D%22SCS%20Parameters%22%2C%0A%20%20%20%20%20%20%20%20selection%3DNone%2C%0A%20%20%20%20%20%20%20%20show_column_summaries%3DFalse%2C%0A%20%20%20%20%20%20%20%20show_download%3DFalse%2C%0A%20%20%20%20)%0A%0A%20%20%20%20_pcallout%20%3D%20mo.callout(table)%0A%20%20%20%20button%20%3D%20mo.ui.run_button(label%3D%22Calculate%22)%0A%0A%20%20%20%20_pstack%20%3D%20mo.vstack(%5B_pcallout%2C%20mo.right(button)%5D)%0A%20%20%20%20mo.accordion(%7B%22Input%20Parameters%22%3A%20_pstack%7D)%0A%20%20%20%20return%20button%2C%20table%0A%0A%0A%40app.cell%0Adef%20_(%0A%20%20%20%20button%2C%0A%20%20%20%20calculate%2C%0A%20%20%20%20clear_button%2C%0A%20%20%20%20empty_file%2C%0A%20%20%20%20mo%2C%0A%20%20%20%20os%2C%0A%20%20%20%20pickle%2C%0A%20%20%20%20plot%2C%0A%20%20%20%20table%2C%0A)%3A%0A%20%20%20%20def%20execute(table%2C%20all_data%3DNone%2C%20events%3DNone)%3A%0A%0A%20%20%20%20%20%20%20%20_flood%20%3D%20calculate(table)%0A%20%20%20%20%20%20%20%20_peak_flow%20%3D%20max(_flood.QNonZero)%0A%0A%20%20%20%20%20%20%20%20_rec%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20'No'%3A%20None%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20'Area'%3A%20float(table.data%5B0%5D%5B%22Value%22%5D.value)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20'Lag'%3A%20float(table.data%5B1%5D%5B%22Value%22%5D.value)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20'Depth'%3A%20table.data%5B2%5D%5B%22Value%22%5D.value%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20'CN'%3A%20float(table.data%5B3%5D%5B%22Value%22%5D.value)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20'Type'%3A%20int(table.data%5B4%5D%5B%22Value%22%5D.value%5B0%5D)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20'Peak%20Flow'%3A%20round(_peak_flow%2C%201)%0A%20%20%20%20%20%20%20%20%7D%0A%0A%20%20%20%20%20%20%20%20events.append(_flood)%0A%20%20%20%20%20%20%20%20all_data.append(_rec)%0A%20%20%20%20%20%20%20%20for%20i%2C%20_item%20in%20enumerate(all_data)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_item%5B'No'%5D%20%3D%20i%2B1%0A%0A%20%20%20%20%20%20%20%20_plt%20%3D%20plot(events)%0A%0A%20%20%20%20%20%20%20%20summary_table%20%3D%20mo.ui.table(%0A%20%20%20%20%20%20%20%20%20%20%20%20data%3Dall_data%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20label%3D'SCS%20Flood%20Events'%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%23_callout%20%3D%20mo.callout(summary_table)%0A%20%20%20%20%20%20%20%20clear_button%20%3D%20mo.ui.run_button(label%3D%22Clear%20Table%22)%0A%20%20%20%20%20%20%20%20clear_button%20%3D%20mo.right(clear_button)%0A%20%20%20%20%20%20%20%20_results%20%3D%20mo.vstack(%5B_plt%2C%20summary_table%5D)%0A%0A%20%20%20%20%20%20%20%20return%20(_results%2C%20all_data%2C%20events)%0A%0A%20%20%20%20results%20%3D%20None%0A%0A%20%20%20%20if%20button.value%3A%0A%20%20%20%20%20%20%20%20if%20not%20os.path.isfile(%22.%2Fevents.pkl%22)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20empty_file()%0A%20%20%20%20%20%20%20%20with%20open('.%2Fevents.pkl'%2C%20'rb')%20as%20dbfile%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20db%20%3D%20pickle.load(dbfile)%0A%20%20%20%20%20%20%20%20if%20db%5B'Data'%5D%20is%20not%20None%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20all_data%20%3D%20db%5B'Data'%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20floods%20%3D%20db%5B'Events'%5D%0A%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20all_data%20%3D%20list()%0A%20%20%20%20%20%20%20%20%20%20%20%20floods%20%3D%20list()%0A%20%20%20%20%20%20%20%20results%2C%20all_data%2C%20floods%20%3D%20execute(table%2C%20all_data%2C%20floods)%0A%20%20%20%20%20%20%20%20db%20%3D%20%7B'Data'%3A%20all_data%2C%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20'Events'%3A%20floods%7D%0A%20%20%20%20%20%20%20%20with%20open('.%2Fevents.pkl'%2C%20'wb')%20as%20dbfile%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20pickle.dump(db%2C%20dbfile)%0A%0A%20%20%20%20if%20clear_button.value%3A%0A%20%20%20%20%20%20%20%20results%20%3D%20None%0A%0A%20%20%20%20results%0A%20%20%20%20return%20all_data%2C%20db%2C%20dbfile%2C%20execute%2C%20floods%2C%20results%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20clear_button%20%3D%20mo.ui.run_button(label%3D%22Clear%20Table%22)%0A%20%20%20%20mo.right(clear_button)%0A%20%20%20%20return%20(clear_button%2C)%0A%0A%0A%40app.cell%0Adef%20_(pickle)%3A%0A%20%20%20%20def%20empty_file()%3A%0A%20%20%20%20%20%20%20%20_db%20%3D%20%7B%22Data%22%3A%20None%2C%20%22Events%22%3A%20None%7D%0A%20%20%20%20%20%20%20%20with%20open(%22.%2Fevents.pkl%22%2C%20%22wb%22)%20as%20_dbfile%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20pickle.dump(_db%2C%20_dbfile)%0A%20%20%20%20return%20(empty_file%2C)%0A%0A%0A%40app.cell%0Adef%20_(clear_button%2C%20empty_file)%3A%0A%20%20%20%20if%20clear_button.value%3A%0A%20%20%20%20%20%20%20%20empty_file()%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(SCS)%3A%0A%20%20%20%20def%20calculate(table)%3A%0A%0A%20%20%20%20%20%20%20%20_area%20%3D%20float(table.data%5B0%5D%5B%22Value%22%5D.value)%0A%20%20%20%20%20%20%20%20_lag%20%3D%20float(table.data%5B1%5D%5B%22Value%22%5D.value)%0A%20%20%20%20%20%20%20%20_depth%20%3D%20table.data%5B2%5D%5B%22Value%22%5D.value%0A%20%20%20%20%20%20%20%20_cn%20%3D%20float(table.data%5B3%5D%5B%22Value%22%5D.value)%0A%20%20%20%20%20%20%20%20_type%20%3D%20int(table.data%5B4%5D%5B%22Value%22%5D.value%5B0%5D)%0A%0A%20%20%20%20%20%20%20%20_flood%20%3D%20SCS(%0A%20%20%20%20%20%20%20%20%20%20%20%20Area%3D_area%2C%20Lag%3D_lag%2C%20RainDepth%3D_depth%2C%20%0A%20%20%20%20%20%20%20%20%20%20%20%20DistType%3D_type%2C%20CN%3D_cn%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_flood.DepthInc()%0A%20%20%20%20%20%20%20%20_flood.Discharge()%0A%0A%20%20%20%20%20%20%20%20return%20_flood%0A%20%20%20%20return%20(calculate%2C)%0A%0A%0A%40app.cell%0Adef%20_(plt)%3A%0A%20%20%20%20def%20plot(floods%3A%20list)%3A%0A%0A%20%20%20%20%20%20%20%20plt.figure(figsize%3D(12%2C%205))%0A%20%20%20%20%20%20%20%20for%20i%2C%20_flood%20in%20enumerate(floods)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20plt.plot(_flood.TNonZero%2C%20_flood.QNonZero%2C%20label%20%3D%20str(i%2B1))%0A%0A%20%20%20%20%20%20%20%20plt.xlabel('T%20(min)')%0A%20%20%20%20%20%20%20%20plt.ylabel('Q%20(m3%2Fs)')%0A%20%20%20%20%20%20%20%20plt.title('Runoff%20Hydrograph')%0A%20%20%20%20%20%20%20%20plt.legend(loc%3D'best')%0A%20%20%20%20%20%20%20%20plt.grid()%0A%20%20%20%20%20%20%20%20%23%20plt.savefig('%2Fhome%2Fandre%2FDocuments%2FJobs%2FWork%20Search%2FSoweto%2FOutput%2FUnit%20Hydrograph%204%20Extra.svg')%0A%20%20%20%20%20%20%20%20%23%20plt.show()%0A%0A%20%20%20%20%20%20%20%20return%20plt.gca()%0A%20%20%20%20return%20(plot%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_()%3A%0A%20%20%20%20class%20SCS%3A%0A%0A%20%20%20%20%20%20%20%20def%20__init__(self%2C%20Area%2C%20Lag%2C%20RainDepth%2C%20DistType%2C%20CN)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Assign%20the%20Area%2C%20Rainfall%20Depth%20and%20CN%0A%20%20%20%20%20%20%20%20%20%20%20%20self.CN%20%3D%20CN%0A%20%20%20%20%20%20%20%20%20%20%20%20self.Area%20%3D%20Area%0A%20%20%20%20%20%20%20%20%20%20%20%20self.RainDepth%20%3D%20RainDepth%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Determine%20time%20increments%0A%20%20%20%20%20%20%20%20%20%20%20%20self.TInc%20%3D%20Lag%20*%2060.0%20%2F%205.5%0A%20%20%20%20%20%20%20%20%20%20%20%20Xstep%20%3D%20(720.0%20%2B%200.5%20*%20self.TInc)%20%2F%20self.TInc%0A%20%20%20%20%20%20%20%20%20%20%20%20Nstep%20%3D%20int(Xstep)%0A%20%20%20%20%20%20%20%20%20%20%20%20self.Time%20%3D%20(Xstep%20-%20Nstep)%20*%20self.TInc%0A%20%20%20%20%20%20%20%20%20%20%20%20self.NoInc%20%3D%202%20*%20Nstep%20%2B%201%0A%20%20%20%20%20%20%20%20%20%20%20%20self.Lag%20%3D%20Lag%0A%20%20%20%20%20%20%20%20%20%20%20%20self.S%20%3D%20(25400.0%20%2F%20self.CN)%20-%20254.0%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Constants%20for%20the%20Rainfall%20Distribution%20type%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20DistType%20%3D%3D%201%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20self.A%2C%20self.B%2C%20self.C%20%3D%200.29935%2C%200.059%2C%200.62%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20DistType%20%3D%3D%202%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20self.A%2C%20self.B%2C%20self.C%20%3D%200.45321%2C%200.1%2C%200.75%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20DistType%20%3D%3D%203%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20self.A%2C%20self.B%2C%20self.C%20%3D%200.73402%2C%200.23%2C%200.9%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20DistType%20%3D%3D%204%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20self.A%2C%20self.B%2C%20self.C%20%3D%201.0133%2C%200.32%2C%201.0%0A%0A%20%20%20%20%20%20%20%20def%20Q(self%2C%20P)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20c%20%3D%200.1%20%20%23%20Suggested%20value%20for%20initial%20losses%0A%20%20%20%20%20%20%20%20%20%20%20%20Ia%20%3D%20c%20*%20self.S%0A%20%20%20%20%20%20%20%20%20%20%20%20Q%20%3D%20((P%20-%20Ia)%20**%202)%20%2F%20(P%20-%20Ia%20%2B%20self.S)%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20Q%0A%0A%20%20%20%20%20%20%20%20def%20DepthInc(self)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Compute%20the%20Rainfall%20Depth%20for%20each%20time%20increment%0A%20%20%20%20%20%20%20%20%20%20%20%20Tmid%20%3D%2012.0%0A%20%20%20%20%20%20%20%20%20%20%20%20self.T%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20self.Qp%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20S%20%3D%20self.S%0A%20%20%20%20%20%20%20%20%20%20%20%20QPrev%20%3D%200%20%20%23%20Not%20initialized%20in%20example%20program%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20a%20in%20range(1%2C%20self.NoInc%20%2B%2014)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20self.T.append(self.Time)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20Diff%20%3D%20Tmid%20-%20self.Time%20%2F%2060.0%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20Diff%20%3E%200%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20Deld%20%3D%202.0%20*%20Diff%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20Ratio%20%3D%20(self.A%20*%20Deld)%20%2F%20((self.B%20%2B%20Deld)%20**%20self.C)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20Prat%20%3D%20(1.0%20-%20Ratio)%20%2F%202.0%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20elif%20Diff%20%3C%200%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20Deld%20%3D%20-2.0%20*%20Diff%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20Ratio%20%3D%20(self.A%20*%20Deld)%20%2F%20((self.B%20%2B%20Deld)%20**%20self.C)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20Ratio%20%3E%201%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%20Ratio%20%3D%201.0%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20Prat%20%3D%201.0%20-%20((1.0%20-%20Ratio)%20%2F%202.0)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20self.Time%20%2B%3D%20self.TInc%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20PAcc%20%3D%20Prat%20*%20self.RainDepth%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%23%20Compute%20incremental%20Runoff%20depth%20for%20each%20time%20step%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20QPres%20%3D%20self.Q(PAcc)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20(PAcc%20-%200.1%20*%20self.S)%20%3C%200%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20QPres%20%3D%200.0%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20Q%20%3D%20QPres%20-%20QPrev%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20self.Qp.append(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%200.2083%20*%20self.Area%20*%20Q%20%2F%20(self.Lag%20%2B%20self.TInc%20%2F%20120.0)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20QPrev%20%3D%20QPres%0A%0A%20%20%20%20%20%20%20%20def%20Discharge(self)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20AK%20%3D%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%200.167%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%200.333%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%200.5%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%200.667%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%200.833%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%201.0%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%200.9%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%200.8%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%200.7%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%200.6%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%200.5%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%200.4%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%200.3%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%200.2%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%200.1%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20self.Qsyn%20%3D%20%5B0.0%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20self.Qpmax%20%3D%200%0A%20%20%20%20%20%20%20%20%20%20%20%20self.QNonZero%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20self.TNonZero%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20b%20in%20range(1%2C%20self.NoInc%20%2B%2014)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20b%20%3C%3D%2015%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20L%20%3D%201%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20b%20%3E%2015%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20L%20%3D%20b%20-%2014%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20for%20c%20in%20range(L%20-%201%2C%20b)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20k%20%3D%20b%20%2B%201%20-%20c%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20self.Qsyn%5Bb%20-%201%5D%20%2B%3D%20AK%5Bk%20-%202%5D%20*%20self.Qp%5Bc%5D%20%20%23%20m**3%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20self.Qsyn%5Bb%20-%201%5D%20%3E%20self.Qpmax%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20self.Qpmax%20%3D%20self.Qsyn%5Bb%20-%201%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20self.Qsyn%5Bb%20-%201%5D%20%3E%200%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20self.QNonZero.append(self.Qsyn%5Bb%20-%201%5D)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20self.TNonZero.append(self.T%5Bb%20-%201%5D)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20self.Qsyn.append(0.0)%0A%20%20%20%20return%20(SCS%2C)%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20import%20marimo%20as%20mo%0A%20%20%20%20return%20(mo%2C)%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20import%20os%0A%20%20%20%20import%20math%0A%20%20%20%20import%20pickle%0A%0A%20%20%20%20import%20micropip%0A%20%20%20%20import%20matplotlib.pyplot%20as%20plt%0A%20%20%20%20return%20math%2C%20micropip%2C%20os%2C%20pickle%2C%20plt%0A%0A%0A%40app.cell%0Aasync%20def%20_(micropip)%3A%0A%20%20%20%20import%20sys%0A%0A%20%20%20%20if%20sys.platform%20%3D%3D%20'emscripten'%3A%0A%20%20%20%20%20%20%20%20await%20micropip.install('tabulate')%0A%20%20%20%20%20%20%20%20from%20tabulate%20import%20tabulate%0A%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20from%20tabulate%20import%20tabulate%0A%20%20%20%20return%20sys%2C%20tabulate%0A%0A%0A%40app.cell%0Adef%20_(micropip)%3A%0A%20%20%20%20from%20pprint%20import%20pprint%0A%20%20%20%20pprint(dir(micropip))%0A%20%20%20%20return%20(pprint%2C)%0A%0A%0A%40app.cell%0Adef%20_()%3A%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