ASY-Forum dédié à Asymptote

[cjorssen]

Bonjour,

Je cherche à représenter la fonction suivante :



pour 0 < d < 5 et surtout (c'est là que je ne sais pas faire) pour -d < z < d (les bornes sur la première variable dépendent de la seconde).

Un idée ?

Merci d'avance

--
Christophe

[OG]

Bonsoir

et si tu composais : avec ?
évidemment il faudra remplacer les intervalles ouverts par des intervalles
fermé ou de ce genre.

O.G.

[GM]

Le sujet a été déplacé dans le forum d'aide.
Comme le suggère la remarque en rouge sur l'index, on va conserver les sous-forums thématiques pour des propositions finalisées qui peuvent être intéressantes.
L'exemple suivant peut en faire partie mais on le mettra quand il sera abouti, tel que tu le veux.

*** Pour masquer/découvrir le code Asymptote qui a permis de créer la figure, il faut cliquer dessus. ;-) ***

CODE ASYMPTOTE de la figure ci-dessus : Tout sélectionner

 
import graph3;
size(7.5cm,0);
size3(7.5cm,IgnoreAspect);
real f(pair z) {return 1/(1+(z.x+z.y/2)^2)^(3/2)+1/(1+(z.x-z.y/2)^2)^(3/2);}
bool cond(pair z) {return abs(z.x)<z.y;}
currentprojection=orthographic(5,5,5);
draw(surface(f,(-5,-5),(5,5),nx=100,Spline,cond),orange+white,meshpen=gray);
 

[GM]

*** Pour masquer/découvrir le code Asymptote qui a permis de créer la figure, il faut cliquer dessus. ;-) ***

CODE ASYMPTOTE de la figure ci-dessus : Tout sélectionner

 
import graph3;
import contour;
currentlight=Viewport;
size(7.5cm,0);
size3(7.5cm,IgnoreAspect);
real f(pair z) {return 1/(1+(z.x+z.y/2)^2)^(3/2)+1/(1+(z.x-z.y/2)^2)^(3/2);}
bool cond(pair z) {return abs(z.x)<z.y;}
currentprojection=orthographic(5,3,8);
draw(surface(f,(-5,-5),(5,5),nx=200,Spline,cond),yellow+white);
 

[OG]

Oui, je n'avais pas pensé à cela : surface(f, blabla) mais aussi graph pour la 2D
accepte une fonction à valeurs booléennes, selon vrai ou faux l'élément est tracé ou non.
Pour une surface cela fait un découpage selon les carreaux de Bézier.

O.G.

[pivaldi]

Oui, je n'avais pas pensé à cela : surface(f, blabla) mais aussi graph pour la 2D
accepte une fonction à valeurs booléennes, selon vrai ou faux l'élément est tracé ou non.
Pour une surface cela fait un découpage selon les carreaux de Bézier.

D'ailleurs, cela ne fait pas très joli. Il faudra penser un jour à faire une routine de lissage.

[cjorssen]

Super. Merci. Pour info, il s'agit de la représentation du champ magnétique sur l'axe de deux spires. La position d=1 correspond à la position dite des bobines de Helmoltz (le champ pour -d/2 < z < d/2, c'est-à-dire entre les deux spires, est quasi constant pour d=1).

PS : promis, la prochaine fois, je poste au bon endroit

--
Christophe

[GM]

Super. Merci. Pour info, il s'agit de la représentation du champ magnétique sur l'axe de deux spires. La position d=1 correspond à la position dite des bobines de Helmoltz (le champ pour -d/2 < z < d/2, c'est-à-dire entre les deux spires, est quasi constant pour d=1).

PS : promis, la prochaine fois, je poste au bon endroit

C'est une figure intéressante, ne serait-ce que pour le coup du booléen donc... reviens mettre le code de ta figure telle que tu la souhaites dans ton cours et on la mettra dans le sous-forum des exemples 3D. Et si tu donnes une petite explication... on l'ajoutera pour les moteurs de recherche.

Tu noteras la différence de nx dans les deux exemples, liée au fait que l'on demande ou pas les lignes : on ne peut pas faire ce que l'on veut car sinon... out of memory.

[OG]

Le coup du booléen est intéressant mais pour des exemples plus compliqués à mon avis.
Là il suffit de faire une surface paramétrée. Pas besoin d'un nx démesuré, et la coupe est lisse.

*** Pour masquer/découvrir le code Asymptote qui a permis de créer la figure, il faut cliquer dessus. ;-) ***

CODE ASYMPTOTE de la figure ci-dessus : Tout sélectionner

 
 
import graph3;
import contour;
currentlight=Viewport;
size(7.5cm,0);
size3(7.5cm,IgnoreAspect);
real f(pair z) {return 1/(1+(z.x+z.y/2)^2)^(3/2)+1/(1+(z.x-z.y/2)^2)^(3/2);}
triple g(pair z)
{
  return (z.x*z.y,z.y,f((z.x*z.y,z.y)));
}
 
bool cond(pair z) {return abs(z.x)<z.y;}
currentprojection=orthographic(5,3,8);
draw(surface(g,(-1,0),(1,5),nu=20,Spline),yellow+white);
//draw(surface(f,(-5,-5),(5,5),nx=100,Spline,cond),red+white);
 

O.G.

[GM]

la coupe est lisse.

C'est vrai que c'est plus joli.