Sage code for trigonometric forcing


omegaf = 3

A(t) = 1-(1/(4-omegaf^2))*cos((omegaf-2)*t)
h(t) = cos(2*t)

AmpPlot = plot(A(t), (t,0,limit),linestyle='dashed', color='red',thickness=2)
hPlot = plot(h(t)*A(t), (t,0,limit))

mainPlot = AmpPlot+ hPlot = -5,ymax = 5, axes_labels=['$t$',' '],ticks = [[],[]])<span id="mce_SELREST_start" style="overflow:hidden;line-height:0;"></span>
This entry was posted in Differential equations, Sage, Uncategorized. Bookmark the permalink.

Leave a comment here

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s