Problem with pst-ode

[hvoss49]

pst-ode. uses a global dict and I suppose that this leads to some problems. The first two images are correct, but the last three (on two pages!) are not. Here is only the first page dvips versus luapstricks. The other wrong image is on page two.

dvips:

luapstricks:

\DocumentMetadata{}
\documentclass{article}
\usepackage{pst-3dplot,pst-ode}
\begin{document}

\psset{xunit=0.45cm,yunit=0.45cm}
\begin{pspicture}(-17,-2.25)(5.75,7.85)
%\psframe(-17,-2.25)(5.75,7.85)
\psaxes[labelFontSize=\color{darkgray}\scriptstyle,linecolor=darkgray,
  linewidth=0.75pt,tickcolor=darkgray,Dx=3,Dy=2]{<->}(0,0)(-17,-2)(5,7)%
    [\textcolor{darkgray}{$t$},0][\textcolor{darkgray}{$x$},90]
\pstODEsolve[algebraic]{TY}{(t) 0}{-15}{3}{200}{0 1}{y[1]|t*y[0]}
\listplot[plotstyle=line,linecolor=blue]{TY}
\end{pspicture}

\psset{unit=1cm,algebraic=false}
\begin{pspicture}(-5,-2.35)(5,3.35)
\psset{xunit=0.85cm,yunit=0.85cm,Alpha=25,Beta=8}
\pstThreeDCoor[xMin=-1,xMax=5,yMin=-1,yMax=5,zMin=-1,zMax=3,linewidth=0.2pt]
\psplotThreeD[linecolor=gray,plotstyle=line,yPlotpoints=50,xPlotpoints=150,linewidth=0.25pt](-4,4)(-4,4)
{5 x 2.729 x 2 exp y 2 exp add neg exp mul mul}
\pstThreeDSquare[linecolor=blue,linewidth=0.05pt,fillstyle=solid,fillcolor=blue,opacity=0.15](-4,-4,1.5)(8,0,0)(0,8,0)
\parametricplotThreeD[xPlotpoints=100,linewidth=0.5pt,linecolor=blue](0.33601,1.16101)
{1.5 5 div 2.729 t 2 exp exp mul
t 2 exp 1.5 5 div 2.729 t 2 exp exp mul 2 exp neg add sqrt
1.5}
\parametricplotThreeD[xPlotpoints=100,linewidth=0.5pt,linecolor=blue](0.33601,1.16101)
{1.5 5 div 2.729 t 2 exp exp mul
t 2 exp 1.5 5 div 2.729 t 2 exp exp mul 2 exp neg add sqrt neg
1.5}
\pstThreeDSquare[linecolor=blue,linewidth=0.05pt,fillstyle=solid,fillcolor=blue,opacity=0.15](4,4,0.5)(-8,0,0)(0,-8,0)
\parametricplotThreeD[xPlotpoints=200,linewidth=0.5pt,linecolor=blue](0.102,1.6756)
{0.5 5 div 2.729 t 2 exp exp mul
t 2 exp 0.5 5 div 2.729 t 2 exp exp mul 2 exp neg add sqrt
0.5}
\parametricplotThreeD[xPlotpoints=200,linewidth=0.5pt,linecolor=blue](0.102,1.6756)
{0.5 5 div 2.729 t 2 exp exp mul
t 2 exp 0.5 5 div 2.729 t 2 exp exp mul 2 exp neg add sqrt neg
0.5}
\pstThreeDSquare[linecolor=blue,linewidth=0.05pt,fillstyle=solid,fillcolor=blue,opacity=0.15](4,4,-1)(-8,0,0)(0,-8,0)
\parametricplotThreeD[xPlotpoints=200,linewidth=0.5pt,linecolor=blue](0.20897,1.38953)
{1 neg 5 div 2.729 t 2 exp exp mul
t 2 exp 1 neg 5 div 2.729 t 2 exp exp mul 2 exp neg add sqrt
1 neg}
\parametricplotThreeD[xPlotpoints=200,linewidth=0.5pt,linecolor=blue](0.20897,1.38953)
{1 neg 5 div 2.729 t 2 exp exp mul
t 2 exp 1 neg 5 div 2.729 t 2 exp exp mul 2 exp neg add sqrt neg
1 neg}
\end{pspicture}


\psset{unit=0.3cm}
\begin{pspicture}(-0.5,-13.5)(12.5,14)
\psaxes[labelFontSize=\color{gray}\scriptstyle,linecolor=gray,linewidth=0.75pt,
  tickcolor=gray,dx=3.1415,Dy=3]{<->}(0,0)(-0.5,-13.5)(13.5,13.5)[\textcolor{gray}{$t$},0][\textcolor{gray}{$x$},90]
\pstODEsolve[algebraic]{TY}{(t) 0}{0}{12.6063}{150}{0 0}{y[1]|-y[0]+2*SIN(0.9*t)}
\listplot[plotstyle=line,linecolor=blue]{TY}
\pstODEsolve[algebraic]{TY}{(t) 0}{0}{12.6063}{150}{0 0}{y[1]|-y[0]+2*SIN(t)}
\listplot[plotstyle=line,linecolor=red]{TY}
\pstODEsolve[algebraic]{TY}{(t) 0}{0}{12.6063}{150}{0 0}{y[1]|-y[0]+2*SIN(1.1*t)}
\listplot[plotstyle=line,linecolor=cyan]{TY}
\psline[linecolor=blue,linewidth=0.75pt](1.5,11.5)(2.5,11.5)
\rput[l](2.75,11.5){\scriptsize{$\omega=0,9$}}
\psline[linecolor=red,linewidth=0.75pt](1.5,9.5)(2.5,9.5)
\rput[l](2.75,9.5){\scriptsize{$\omega=1,0$}}
\psline[linecolor=cyan,linewidth=0.75pt](1.5,7.5)(2.5,7.5)
\rput[l](2.75,7.5){\scriptsize{$\omega=1,1$}}
\end{pspicture}
%
\hfill
%
\psset{trigLabels,xunit=0.5cm,yunit=0.15cm}
\begin{pspicture}(-0.5,-20.5)(12.5,14)
\psaxes[labelFontSize=\color{lightgray}\scriptstyle,linecolor=gray,linewidth=0.75pt,tickcolor=lightgray,dx=3.1415,Dy=3]{<->}(0,0)(-0.5,-20)(13.5,13.5)[\textcolor{gray}{$t$},0][\textcolor{gray}{$x$},90]
\pstODEsolve[algebraic]{TY}{(t) 0}{0}{12.6063}{150}{0 0}{y[1]|-y[0]+2*SIN(t)+0.1*y[1]}
\listplot[plotstyle=line,linecolor=blue]{TY}
\pstODEsolve[algebraic]{TY}{(t) 0}{0}{12.6063}{150}{0 0}{y[1]|-y[0]+2*SIN(t)}
\listplot[plotstyle=line,linecolor=cyan]{TY}
\pstODEsolve[algebraic]{TY}{(t) 0}{0}{12.6063}{150}{0 0}{y[1]|-y[0]+2*SIN(t)-1.1*y[1]}
\listplot[plotstyle=line,linecolor=red]{TY}
\psline[linecolor=blue,linewidth=0.75pt](1.5,-7.5)(2.5,-7.5)
\rput[l](2.75,-7.5){\scriptsize{$a=-0,1$}}
\psline[linecolor=cyan,linewidth=0.75pt](1.5,-9.5)(2.5,-9.5)
\rput[l](2.75,-9.5){\scriptsize{$a=0$}}
\psline[linecolor=red,linewidth=0.75pt](1.5,-11.5)(2.5,-11.5)
\rput[l](2.75,-11.5){\scriptsize{$a=1,1$}}
\end{pspicture}



\psset{unit=1cm}
\begin{pspicture}(-2,0)(2,4.5)
\pstVerb{
 /alphaP 10 def
 /betaP 28 def
 /gammaP 8 3 div def
}%
\pstODEsolve[algebraic]{lorenzXYZ}{0 1 2}{0}{50}{6000}{10 10 30}{
 alphaP*(x[1]-x[0]) |
 x[0]*(betaP-x[2]) - x[1] |
 x[0]*x[1] - gammaP*x[2]
}%
\pstODEsolve[algebraic]{lorenzRED}{0 1 2}{0}{50}{6000}{12 12 28}{
 alphaP*(x[1]-x[0]) |
 x[0]*(betaP-x[2]) - x[1] |
 x[0]*x[1] - gammaP*x[2]
}
\psset{unit=0.1cm,Alpha=-2,Beta=15}
\listplotThreeD[linecolor=red!70,linewidth=0.25pt]{lorenzRED}
\listplotThreeD[linecolor=blue!70,linewidth=0.25pt]{lorenzXYZ}
\end{pspicture}

\end{document}