x0 = get_x(petab_prob) res = calibrate(petab_prob, x0, IPNewton()) @show res.xmin using Plots plot(res, petab_prob; linewidth = 2.0) Results in the following ...