print("Simulating model for figure 4B...")
t = linspace(0, 6*3600, 6*60+1) // 6 hours
x = odesolve(model, t)
x_norm = c_[x["Bid_unbound"], x["PARP_unbound"], x["mSmac_unbound"]]
x_norm = 1 - x_norm / x_norm[0, :] // gets away without max() since first values are largest
After Change
print("Simulating model for figure 4B...")
t = linspace(0, 6*3600, 6*60+1) // 6 hours
x = ScipyOdeSimulator(model).run(tspan=t).all
x_norm = c_[x["Bid_unbound"], x["PARP_unbound"], x["mSmac_unbound"]]
x_norm = 1 - x_norm / x_norm[0, :] // gets away without max() since first values are largest