N charges ponctuelles dans l'espace

Posté par brocoli brocoli

Bonsoir,

J'ai fait un programme maple, qui me permet de tracer dans le plan les lignes de champs et les lignes d'équipotentielles.
Cela fonctionne très bien, mais à présent je voudrait le faire dans l'espace. L'ennui c'est que j'utilise la fonction fieldplot, qui ne fonctionne que dans le plan...
DOnc si quelqu'un pouvait me dire quelle modification faire pour arranger ca!
voici le code :

   restart:with(plots):with(linalg):
    roll := rand(-2..2):
    > rolll:=rand(-5..5):
    > rollll:=rand(-5..5):
    > for i from 0 to 15 do
    > rq[i]:=[rolll(),rollll(),roll()]:
    > od:
    > V:=(x0,y0,q,x,y)->q/sqrt((x-x0)^2+(y-y0)^2):
    > V1:=sum(V(rq[k][1],rq[k][2],rq[k][3],x,y),k=0..15):
    > leqpot:=contourplot(V1,x=-10..10,y=-10..10,grid=[100,100],scaling=constrained):
    > E:=grad(-V1,[x,y]):
    > Eunitaire:=evalm(E/norm(E,2)):
    > ldc:=fieldplot(Eunitaire,x=-10..10,y=-10..10,grid=[100,100],scaling=constrained):
    > display(leqpot,ldc);