1112.sagews

x = PolynomialRing(QQ, 'x').gen()
f(x) = hermite(5,x)
p = plot(f(x), (x,-2,2))
p.show()

var('x,y')
f(x,y) = x^2*cos(x-y) - y^2*sin(x+y)
p = plot3d(f(x,y), (x,0,2*pi), (y,0,2*pi))
p.show()

f1(t) = 6*cos(2*pi*100*t-5*pi*1)
f2(t) = 6*cos(2*pi*110*t-5*pi*1)
f(t) = f1(t) + f2(t)
q = plot(f(x), (t,0,0.5))
q.show()
plot

<function plot at 0x2ea0668>

f(x) = (3 * x^4 + 4 * x^3 + 16 * x^2 + 20 * x + 9) / ((x + 2) * (x^2 + 3)^2)
g(x) = f.partial_fraction(x)
show(g)

\newcommand{\Bold}[1]{\mathbf{#1}}x \ {\mapsto}\ \frac{2 \, x}{x^{2} + 3} + \frac{1}{x + 2} + \frac{4 \, x}{{\left(x^{2} + 3\right)}^{2}}