write a script in matlab

The mechanical power output P in a contracting muscle is given by Umax 0 uma where T' is the muscle tension, v is the shortening velocity (max of vmaz), To is the isometric tension (i.e., tension at zero velocity), and k is a non-dimensional constant that ranges between 0.15 and 0.25 for most muscles. The equation can be written in the non-dimensional form ku(1-u) where p = (Tu)/(T,Umar) and u = u/umaz. A figure with k format long for parts (c) and (d) 0.25 is shown below. Use 0.1 0.08 0.06 0.04 0.02 0 0.2 04 0. 08I a) Create a vector u ranging from 0 to 1 with increments of 0.05 (b) Using k 0.25, calculate the value of p for each value of u. (c) Find the maximum value of p. (d) Repeat the first three steps with increments of 0.01 and calculate the percent relative error, defined by pma 0.01-pma0.051 100 Pmar0.05 Solution; (c) pmaz0.05-0.095454545454545 and Pmaz0.01-0.095491071428571, and (d) E 0.038265306122427