Hill-type muscle model.ipynb

Forces chapter. "$F_{SEE}$: parallel element force. calculated"
Isn't it serial element force?

Yup, that’s correct. Just a typo.

Thank you for writing this code, which greatly helped me to learn Hill model.
It seems that there is a problem with the cell 9:

def f_PEE(L_CE_norm, V_CE_norm): 
    return F_max * ((L_CE_norm - L_CE_opt) / (L_CE_opt * w)) * f_V(V_CE_norm)

should be

def f_PEE(L_CE_norm, V_CE_norm): 
    return F_max * np.power((L_CE_norm - L_CE_opt) / (L_CE_opt * w), 2) * f_V(V_CE_norm)

Refer to equation 7

Hmm. That's possibly a typo or I may have removed it on purpose. It's been so long I don't remember.

Thank you for writing this code, which greatly helped me to learn Hill model. It seems that there is a problem with the cell 9:

def f_PEE(L_CE_norm, V_CE_norm): 
    return F_max * ((L_CE_norm - L_CE_opt) / (L_CE_opt * w)) * f_V(V_CE_norm)

should be

def f_PEE(L_CE_norm, V_CE_norm): 
    return F_max * np.power((L_CE_norm - L_CE_opt) / (L_CE_opt * w), 2) * f_V(V_CE_norm)

Refer to equation 7

The power is calculated in the following if