# monod_fit.py — verstrekte synthetische snelheidsmetingen
import numpy as np
from scipy.optimize import least_squares
S=np.array([.05,.1,.2,.5,1,2,5]) # g/L
gemeten=.4*S/(.5+S)+np.array([-.001,.002,-.002,.001,.002,-.001,0])
fit=least_squares(lambda z:z[0]*S/(z[1]+S)-gemeten,
    [.35,.4],bounds=([.001,.001],[2,10]))
mu,K=fit.x
print(f"mu_max={mu:.6f} h^-1; Ks={K:.6f} g/L")
print("residuen:",fit.fun)
print("max residu h^-1:",np.max(abs(fit.fun)))
print("Geen biologisch betrouwbaarheidsinterval uit alleen deze zeven lespunten.")