I’ve attached a file that has a pdf explaining the issue and the solution. I couldn't attach the .m file with the solution. Here's a brief description of the solution. In the function implement_time_step in the muscle class,
-
commented out lines 80 to 83
Slack mode in MATMyoSim02.pdf
Slack mode in MATMyoSim02.pdf
-
added a comment describing my commenting out in lines 78 and 79
-
moved the commented-out lines to line 99 to 102.
-
Made a small change to the input values to implement_step_time in the half-sarcomere class as shown here:
obj.hs(1).implement_time_step(time_step, adjustment, Ca_value, m_props)
uses “adjustment” as the harcomere length change as opposed to obj.hs(1).implement_time_step(time_step, 0, Ca_value, m_props)
which was originally happening in line 80. By calling it with “adjustment”, it also runs move_cb_distribution in implement_time_step in the half_sarcomere class (line 11).
-
commented out line 103 since this already calls move_cb_distribution
I’ve attached a file that has a pdf explaining the issue and the solution. I couldn't attach the .m file with the solution. Here's a brief description of the solution. In the function implement_time_step in the muscle class,
commented out lines 80 to 83
Slack mode in MATMyoSim02.pdf
Slack mode in MATMyoSim02.pdf
added a comment describing my commenting out in lines 78 and 79
moved the commented-out lines to line 99 to 102.
Made a small change to the input values to implement_step_time in the half-sarcomere class as shown here:
obj.hs(1).implement_time_step(time_step, adjustment, Ca_value, m_props)
uses “adjustment” as the harcomere length change as opposed to obj.hs(1).implement_time_step(time_step, 0, Ca_value, m_props)
which was originally happening in line 80. By calling it with “adjustment”, it also runs move_cb_distribution in implement_time_step in the half_sarcomere class (line 11).
commented out line 103 since this already calls move_cb_distribution