%TITRE{Maxima - Surfaces (1)}
%VTEX{\entete}
%ITEX{env=quote}{../note.tex}
%VTEX{\vspace{1cm}}
<[
<style>
 .out {
  background-color: white;
  color: blue;
 }
</style>
]>

======================================== Script pour adapter la couleur
FICHIER:chgcouleur.pl::n:
#!/usr/bin/perl
open(SRC,"maxout.ps");
@lignes = <SRC>;
close(SRC);
open(DEST,">$ARGV[0]");
$lignes[41] = "/xr $ARGV[1] def \n";
$lignes[42] = "/xg $ARGV[2] def \n";
$lignes[43] = "/xb $ARGV[3] def \n";
foreach $l (@lignes) {
    print DEST $l;
}
close(DEST);
§
SH:chmod u+x chgcouleur.pl

======================================== La macro de la fabrique
FICHIER:eps::n:
%!
§
FICHIER:image.txt::n:
<[ 
<pre class="out">
]> 
<[fig%%n%%.txt]> 
<[ 
</pre> 
]> 
%ITEX{env=verbatim}{fig%%n%%.txt}
FICHIER:fig%%n%%.tex:
\begin{center} %
\includegraphics[scale=0.5]{fig%%n%%.eps} %
\end{center} %
%§
SH:maxima -b fig%%n%%.txt
SH:./chgcouleur.pl fig%%n%%.ps %%couleur%%
SH:convert -crop 0x0 -geometry 350x350 fig%%n%%.ps fig%%n%%.jpg
SH:gs -sDEVICE=bbox -dNOPAUSE -dQUIET fig%%n%%.ps quit.ps 2> bbox
SH:cat eps bbox fig%%n%%.ps > fig%%n%%.eps
SH:rm -f maxout.ps fig%%n%%.ps bbox
<[
<div align="center">
 <img src="fig%%n%%.jpg" alt="fig%%n%%">
</div>
]>
§

%P{Voici une première série de surfaces générées par §gMaxima§. Les images
ci-dessous sont obtenues à partir du fichier §vmaxout.ps§, leurs couleurs
ont été adaptées, pour introduire un peu de variété...}
 
%S{Surfaces en Z}

FICHIER:fig1.txt::n:
plot3d(cos(x)*cos(y),[x,-%pi,%pi],[y,-%pi,%pi],
       ['colour_z,true],['plot_format,ps]);
§

M:image: n="1" couleur=".6 .9 .2"

FICHIER:fig2.txt::n:
plot3d(sin(x^2+y^2)/(x^2+y^2), [x,-%pi,%pi],[y,-%pi,%pi],
       ['plot_format,ps],['colour_z,true]);
§

M:image: n="2" couleur=".3 1 .4"

%S{Le ruban de Möbius}

FICHIER:fig3.txt::n:
plot3d([cos(x)*(3+y*cos(x/2)),sin(x)*(3+y*cos(x/2)),y*sin(x/2)],[x,-%pi,%pi],[y,-1,1],
       ['plot_format,ps],['colour_z,true],['grid,50,15]);
§

M:image: n="3" couleur=".4 .1 1"

%S{La bouteille de Klein}

FICHIER:fig4.txt::n:
plot3d([5*cos(x)*(cos(x/2)*cos(y)+sin(x/2)*sin(2*y)+3.0)-10.0,
        -5*sin(x)*(cos(x/2)*cos(y)+sin(x/2)*sin(2*y)+3.0),
        5*(-sin(x/2)*cos(y)+cos(x/2)*sin(2*y))],[x,-%pi,%pi],[y,-%pi,%pi],
        ['plot_format,ps],['colour_z,true],['grid,30,30]);
§

M:image: n="4" couleur="1 .2 0"

%S{Le tore}

FICHIER:fig5.txt::n:
plot3d([cos(y)*(10.0+6*cos(x)),sin(y)*(10.0+6*cos(x)),-6*sin(x)],[x,0,2*%pi],[y,0,2*%pi],
       ['plot_format,ps],['colour_z,true],['grid,31,31]);
§

M:image: n="5" couleur="1 .1 .8"       

%S{Surface définie en polaires}

FICHIER:fig6.txt::n:
plot3d(r^.3333*cos(t/3),[r,0,1],[t,0,6*%pi],['transform_xy,polar_to_xy],
       ['colour_z,true],['view_direction,1,1,1.4],['plot_format,ps],['grid,12,75]);
§

M:image: n="6" couleur="0 0 .5"

%S{Une selle}

FICHIER:fig7.txt::n:
plot3d(2^(x^2-y^2),[x,-1.5,1.5],[y,-2,2],['plot_format,ps],['colour_z,true]);
§

M:image: n="7" couleur="0 .5 .2"