Upload New File

ODE_template.py

+#!/usr/bin/env python3
+# -*- coding: utf-8 -*-
+"""
+Created on Mon Feb 19 12:12:47 2024
+
+@author: delvigne
+"""
+import numpy as np
+from scipy.integrate import solve_ivp
+from matplotlib.pyplot import figure
+from matplotlib.pyplot import plot
+from matplotlib.pyplot import xlabel
+from matplotlib.pyplot import ylabel
+
+'''
+This code doesn't work, it is only a template'
+'''
+
+#Specify your parameters
+
+V = ...
+Q = ...
+P = ...
+X = ...
+...
+
+#Enter your ODE(s) as a new function in Python using def
+
+def PFR(t, x):
+    yourODEs
+
+r1 = solve_ivp(fun=yourODEs, t_span=[...,...], y0=[...], t_eval=np.arange(...,...,...))
+
+
+#Plot the simulation
+figure(1)
+plot(r1.t,r1.y.T)
+xlabel('Time (min)')
+ylabel('Concentration')
+