Fichiers Sage
%md # Paraboloïde
var('x, y, z') F(x,y,z) = x^2+y^2-z show(F)
r = 10 G = implicit_plot3d(F(x,y,z)==0, (x, -r, r), (y, -r, r), (z, -r, r), color='orange', mesh=1, opacity=.7, spin=1) show(G)
%md # Droite normale
gradF = F.gradient() show(gradF)
gradF0 = gradF(1, 1, 2) show(gradF0)
var('t') t_max=5 G += parametric_plot3d([1+2*t, 1+2*t, 2-t], (t, -t_max, t_max), color='red') show(G)
%md # Intersections
var('t') x(t) = 1+2*t y(t) = 1+2*t z(t) = 2-t eq = x(t)^2+y(t)^2 == z(t) show(eq)
solve(eq, t)
P(t) = ([x(t), y(t), z(t)]) show(P)
P0 = P(0) show(P0)
Q = P(-9/8) show(Q)