Thanks for reaching out.
In the MATLAB Live Editor file, we have taken some assumptions into consideration. So, if you are looking for perfect fit you can try the following things:
- Include coefficient C in the magic formula instead of assuming it to be equal to 1/B
- You can also try with exact value of D
- Try few more simulations with different initial conditions
Let us know how does it work.