f  = sin(x) + x^2*cos(x) + sqrt(x)
fp = f.plot(0, 15, thickness=2)
tp = sum((f.taylor(x,0,i).plot(0, 15, color='green')) for i in range(1, 4))
show(tp + fp, gridlines=True, figsize=[12,4])

%var x y z
g = golden_ratio; r = 4.77
p = 2 - (cos(x + g*y) + cos(x - g*y) + cos(y + g*z) +
         cos(y - g*z) + cos(z - g*x) + cos(z + g*x))
show(implicit_plot3d(p, (x, -r, r), (y, -r, r), (z, -r, r),
     plot_points=30, color='green', mesh=1, opacity=.9))