16-2 Line Integrals.sagews

var('t', 'x', 'y', 'z')
#f(x, y) = 1 + x + y/2
#f(x, y) = 1 + sin(x) + cos(y)
f(x, y) = 1.4^x - 1.8^y
A = plot3d(f(x, y), (x, 0, 3), (y, 0, 3))
C(t) = [2*cos(t), 2*sin(t)]
B = parametric_plot3d([C(t)[0], C(t)[1], 0], (t, 0, pi/2), color = 'red')
C = parametric_plot3d([C(t)[0], C(t)[1], z*f(C(t)[0], C(t)[1])], (t, 0, pi/2), (z, 0, 1), color = 'green')
D = plot3d(0, (x, 0, 3), (y, 0, 3), color = 'blue', opacity = 0.1)
show(D + A + B + C)

(t, x, y, z)

3D rendering not yet implemented