You can plot Fig. 12.8.5 with this:
I2 = 1./sqrt((r.^2 - 1).^2 + 0.01);
ylabel('amplitude of i_{2}(t)')
title('The resonance curve')
But actually, the values for the parameters M, R,
,
,
,
, are not given. For some unknown reasons, your teacher expects (but didn't say it out loud) you guys to cleverly assume the values like the author did:
and
.Perhaps the values can be found in one of the references [24] – [27] cited.
Anyhow, the rest you can follow the steps as explained in detailed by the author.
Guess your teacher wants you guys to follow the steps (in ebook) and demostrate them using syms, laplace, solve, and ilaplace to find the solution in Eq. 12.8.60.