@mmcdara Thank you for your swift response. When running your code I noticed several error messages. As you did not have them is have attached them as pictures so you can see which errors I got:

The data was fit to the model by integrating equation 20 of the paper.

The constants are difficult to find in the article as they are not in one location, but I believe they are correct. In the mean time I calculated Gamma_i for IL6 and that one seems to be very close to the value of the paper. I agree with your comment on s versus minutes which is the reason why I would like to incorporate units in the model.

Thank you once again for your precious help!

@mmcdara First, thank you for the additional time to optimize my original attempt. Your background is certainly much better than mine! I started the project when reading the attached "simple model". I digitized the data for TNF (those that I provided in the Maple worksheet) in order to achieve the gamma_i for TNF given in Table 1. My result is far away from the value given by the authors and if I plot their value the plot result is worse than the value I obtain. So I try to understand whether I made a mistake in interpreting the technique used by the authors or that the value given in the article is incorrect. 

If you were willing and of course have some time, could you verify whether I made the wrong interpretation or that the published value is wrong?

DiLeo.pdf

Additionally, I was wondering whether I could use units instead of correcting them as done now.

PS I tried to attach the pdf and the zipped pdf but I am not sure they are attached.

@mmcdara Many thanks

@dharr Thank you for your time and help with my question! 

@mmcdara your assumptions in the beginning of your reply were right. I am very thankful for your guidance as I am still struggling quite a bit with learning Maple. I do have some minor additional questions:

How can you write normalisatie in one line

Why do you have to use 'copy'?

Thanks again for your time and valuable help!

@Thomas Richard  Thanks to all who took the time to respond. I wish all a nice weekend.

@sursumCorda thank you for your swift reply! What I really want to do is to convert the Python code into Maple code. Is that possible? Thanks in advance.

@Carl Love I do apologize that I initially did reduce the data points but once solutions were proposed it made sense to use all data points for a correct comparison. Thank you for your time and insights.

@Carl Love thank you for your response. I am still amazed of what can be done with Maple. The response given by @C_R 1454 is for me closer to what I was used to doing. I appreciate your support to this forum as you help so many people.

@C_R First of all thank you for your patience! You were right in your assumptions. For clearness I have added following explantion:

and

In my original example I didn't provide all data points. I added them in your worksheet and the results gets pretty close to that of Minerr. SolveBlockQuestion-AddData.mw

I have some additional questions: Is this the most efficient way in doing this type of calculation or would it better to make H a function avoiding all ~ and potential typing errors? Thank you again for all help!

@C_R when I do this elementwise I obtain:

@C_R thank you for trying to help me out. The bold = is equal to a normal = in maple. The arrow above stands for vectorization and does an elementwise approach. So in maple it would translate as (1-~(Q/Q0). The Q0, kf, kl and Qs are above the "given" are guess values. Hopefully this makes it a bit more clear.

@acer thank you I was not aware of that.

@Carl Love thank you for the additonal explanation as it was very helpful. I would love to vote you up but do not know how I can do that. Could you explain?

@Carl Love thank you so much, that did the trick. Maple is a great program, but the syntax possibilities are huge. Thanks again for taking time to help me out.