mercredi 11 juillet 2018

Mise à jour de pst-contourplot et pst-2fils

La nouvelle version de pst-contourplot est dans ce répertoire avec les exemples et la documentation :

Le fichier zippé contient le package, la documentation et les fichiers sources des exemples.

Les lignes de champ et les équipotentielles de 2 fils parallèles uniformément chargés tiennent compte de la nouvelle version de pst-contourplot. Les fichiers sont dans le répertoire :


12 juillet 2018 : correction de l'option [algebraic] qui est bien plus rapide maintenant, il n'est plus nécessaire de coder en postscript pour que la compilation s'effectue rapidement (superbe travail de Dominique Rodriguez pour cette option).

15 juillet 2018 : une nouvelle version qui contient des options nouvelles suivantes :
  • [showpoints]  booléen pour afficher les points 
  • [WriteData] : pour enregistrer sur le disque les points de la courbe
ainsi qu'une commande \psReadData[FileName=...] pour dessiner la courbe enregistrée.

La documentation en anglais (pst-contourplot-docEN.pdf) est incluse.

Cette version est actualisée dans le répertoire indiqué au début.

3 exemples on été rajoutés :
Le fichier "exemples-3.tex" est dans le répertoire. En voici le listing, les sources sont indiquées dans les commentaires :

\documentclass[pstricks]{standalone}
\usepackage{pst-contourplot,multido,pst-plot}
\begin{document}

\begin{pspicture}(-6.25,-6.25)(6.25,6.25)
% https://mathematica.stackexchange.com/questions/547/plotting-an-implicit-polar-equation
\psframe(-6.25,-6.25)(6.25,6.25)
\pstVerb{
/arctan {
 3 dict begin
 /x exch def
 /y exch def
 /Arc y x atan def
  Arc 180 ge {/Arc Arc 360 sub def} if
  Arc
 end
} def
}%
\psContourPlot[unit=0.25,a=0.1,function=y x arctan DegToRad dup mul 5.55 x dup mul y dup mul add sqrt RadToDeg cos mul sub,linewidth=0.2,Fill,fillcolor=orange](-25,-25)(25,25)
\end{pspicture}

\begin{pspicture}[showgrid](-0.5,-0.5)(10,10)
% https://stackoverflow.com/questions/42076864/plotting-the-implicit-function-xy-logx-logy-2-0-on-matlab
\multido{\nC=5.6+-0.5,\n=0.0+0.1}{8}{
\psContourPlot[unit=1,a=0.1,algebraic,function=-ln(x)-ln(y)+x+y-\nC,Fill,fillcolor={[hsb]{\n,1,1}}](0.01,0.01)(10,10)}
\end{pspicture}

\begin{pspicture}(-5,-5)(5,5)
% https://mathsbyagirl.wordpress.com/2015/12/04/math-isnt-cool/
\psframe*[linecolor=cyan](-5,-5)(5,5)
\psContourPlot[unit=1,a=0.02,algebraic,function=sin(x^2+y^2)-cos(x*y),Fill,fillcolor=orange](-5,-5)(5,5)
\end{pspicture}
\end{document}

18 juillet 2018 

J'emprunte ces exemples au cours de Julien Royer :

Le fichier source exemple-4.tex est dans le répertoire indiqué au début.
pst-contourplot.tex a été modifié pour assurer la compatibilité de l'option [grid] avec celle de pst-solides3d.

des lignes de niveau, en 3D, avec pst-solides3d :
En projection sur un plan horizontal avec pst-contourplot :
Le listing de ces derniers exemples :

\documentclass[pstricks]{standalone}
\usepackage{pst-solides3d,pst-contourplot}
\begin{document}
\begin{pspicture}(-5,-5)(5,5)
% Julien Royer
% https://www.math.univ-toulouse.fr/~jroyer/TD/2015-16-L2PS/L2PS-poly.pdf
\psframe(-5,-5)(5,5)
\multido{\n=-2.0+.2,\nH=0.0+0.05}{20}{
\psContourPlot[unit=1,a=0.05,algebraic,function=sin(x)-sin(y)-\n,Fill,fillcolor={[hsb]{\nH,1,1}},ReverseColors](-5,-5)(5,5)}
\end{pspicture}

\begin{pspicture}(-5,-5)(5,5)
% Julien Royer
% https://www.math.univ-toulouse.fr/~jroyer/TD/2015-16-L2PS/L2PS-poly.pdf
\psframe(-5,-5)(5,5)
\multido{\n=-2.0+.2,\nH=0.0+0.05}{20}{
\psContourPlot[unit=1,a=0.05,algebraic,function=sin(x)-sin(y)-\n,linecolor={[hsb]{\nH,1,1}}](-5,-5)(5,5)}
\end{pspicture}

\begin{pspicture}(-5,-5)(5,5)
% https://www.math.univ-toulouse.fr/~jroyer/TD/2015-16-L2PS/L2PS-poly.pdf
\psframe(-5,-5)(5,5)
\multido{\n=-1.00+.05,\nH=0.0+0.04}{25}{
\psContourPlot[unit=2,a=0.05,algebraic,function=(x^2-2*y^2)*Euler^(-2*x^2-y^2)-\n,linecolor={[hsb]{\nH,1,1}}](-2.5,-2.5)(2.5,2.5)}
\end{pspicture}

\begin{pspicture}(-5,-5)(5,5)
\psframe(-5,-5)(5,5)
\psset{viewpoint=50 40 20 rtp2xyz,Decran=70,lightsrc=50 30 35 rtp2xyz}
\psSurface[
    fillcolor=white,algebraic,
    intersectionplan={[0 0 1 -0.1] [0 0 1 0.3] [0 0 1 0.2]},
    intersectioncolor=(bleu) (red) (green),
    intersectionlinewidth=1,
    intersectiontype=0,
    ngrid=.1 .1,incolor=yellow!50,linewidth=0.01](-2.5,-2.5)(2.5,2.5){2*(x^2-2*y^2)*Euler^(-2*x^2-y^2)}
\end{pspicture}
\end{document}

20 juillet 2018 : Le package a été déposé sur le serveur du CTAN :









2 commentaires: