It works well but it's slow because of the optim procedure : lam_Nt > calibrate vr > Leslie matrix as f(Nt)