[Maple]Problème Tracé de champ magnétique

Question

Bonjour à tous,

à vrai dire je ne sais pas très bien où placer ce post ... 

J'essaye de tracer les lignes de champ magnétique crées par une spire de courant avec Maple V Release 4 mais je suis bloqué.

Voici ce que j'ai écris :

"restart; 
> with(plots):

> Bx:=(x,y)->9*10^(-1)*int(-y*cos(t)/((x^2+y^2+(15*10^(-2))^2-2*15*10^(-2)*x*cos(t))^(3/2)),t=0..Pi);

                         Pi
                        /
                       |                    y cos(t)
 Bx := (x, y) -> 9/10  |    - ------------------------------------ dt
                       |        2    2                         3/2
                      /       (x  + y  + 9/400 - 3/10 x cos(t))
                        0

> By:=(x,y)->9*10^(-1)*int((x*cos(t)-15*10^(-2))/(x^2+y^2+(15*10^(-2))^2-2*15*10^(-2)*x*cos(t))^(3/2),t=0..Pi);

                          Pi
                         /
                        |              x cos(t) - 3/20
  By := (x, y) -> 9/10  |    ------------------------------------ dt
                        |      2    2                         3/2
                       /     (x  + y  + 9/400 - 3/10 x cos(t))
                         0

> fieldplot([Bx,By],20*10^(-3)..20*10^(-2),20*10^(-3)..15*10^(-2));

Error, (in plots/fieldplot/arrowsf2d) division by zero
"

Ici j'essaye directement de tracer le champ de vecteurs mais je ne vois pas pourquoi maple me répond cette erreur.

J'ai essayé sinon une autre piste en essayant de tracer les courbes intégrales d'une équation différentielle que vérifie les lignes de champ :

"eq:=diff(x(y),y)=By(x(y),y)/(Bx(x(y),y));
                    Pi
                   /
                  |                x(y) cos(t) - 3/20
                  |    ------------------------------------------ dt
                  |         2    2                            3/2
                 /     (x(y)  + y  + 9/400 - 3/10 x(y) cos(t))
      d            0
eq :-- x(y) -----------------------------------------------------
      dy           Pi
                  /
                 |                       y cos(t)
                 |    - ------------------------------------------ dt
                 |           2    2                            3/2
                /       (x(y)  + y  + 9/400 - 3/10 x(y) cos(t))
                  0

> dsolve({eq,x(0)=10^(-3)},x(y));
Error, (in int) wrong number (or type) of arguments
> 
> with(DEtools):
> DEplot({eq},\{y(x)},x=-30*10^(-2)..30*10^(-2),[[y(0)=10^(-3)]],y=-30*10^(-2)..30*10^(-2),stepsize=.05);
Error, (in DEtools/DEplot/CheckDE) Symbolic name, t, could not evaluate to float.
" 


J'arrive encore à des erreurs ...


Si quelqu'un voit ce qui ne va pas.

Merci d'avance.



PS: Dsl pour les grosses intégrales qui ne rentre pas bien dans le petit cadre

Nethernat · Answer

Erf Désolé çà s'est tout packé. 

çà doit être illisible. En copiant collant sur maple peut être que çà devrait remettre le format initial...

[Maple]Problème Tracé de champ magnétique

1 réponse

Votre réponse

Discussions similaires