Sage code for trigonometric forcing

var('t')

omegaf = 3
limit=130

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
mainPlot.show(ymin = -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:

WordPress.com Logo

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

Google+ photo

You are commenting using your Google+ 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