Based on walkingrandomly article
Underdamped - The system oscillates (at reduced frequency compared to the undamped case) with the amplitude gradually decreasing to zero.
f = exp(-(0.0001) * x) f.plot(x,1,40000)
f = 10 * sin(0.01 * x) f.plot(x,1,10000)
Damped sine plot
f = 10 * exp(-(0.0001) * x) * sin(0.01 * x) f.plot(x,1,10000)