MATLAB: I want these kinetic models to fit in MATLAB but I don’t know how to proceed. I am in real worry. If anyone can solve/code, it would be of great help. Thanks a million!!

kinetic model

You need to define Haldane as follows (your text-book expression is misleading)

Haldane = @(S) 0.478*S / ((1.137 + S) .* (1 + (S/1.139)));

Then you get:

which matches the line on your text-book graph.

doc ode45 % please read the documentation

Delete these two lines in your function and it works:

syms x1  x2  x3  x4  x5  v1  v2  v3  v4  v5;
x = [ x1 ,  x2 ,  x3 ,  x4 ,  x5 ,  v1 ,  v2 ,   v3 ,   v4 , v5]';