[pstricks] psplotDiffEqn: Compiles fine but postscript output is not shown
Hernando Diaz Morales
hdiazmo at unal.edu.co
Thu May 27 15:24:02 CEST 2010
Hello everyone:
I am trying to simulate the solution to a 2nd order ODE using
\psplotDiffEqn. The following code compiles just fine, but GSView
crashes with the result.
Please help!
THE CODE:
\listfiles
\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{pstricks}
\usepackage{pstricks-add}
\begin{document}
\psset{whichabs=0,whichord=1}
\psset{method=rk4,plotpoints=200,linewidth=0.5pt,linecolor=blue}
\psset{algebraic=true}
This is the differential equation whose solution is to be simulated
\begin{align}
\dot{X} &=\left(1+X-Y\right)^2\left(1-X\right)-
\theta\left(1-X+Y\right)^2X
-2\eta\left(1-X+Y\right)\left(1+X-Y\right)X\\
\dot{Y}
&=-\left(1+X-Y\right)^2Y+\theta\left(1-X+Y\right)^2\left(1-Y\right)
-2\eta\left(1-X+Y\right)\left(1+X-Y\right)Y
\end{align}
\psset{unit=8cm}
\def\rthe{1.05}
\def\reta{0.95}
\def\Func{(1-y[1])*(1+y[1]-y[2])^(2)-\rthe*y[1]*(1-y[1]+y[2])^(2)-2*\reta*y[1]*(1+y[1]-y[2])*(1-y[1]+y[2])|%
-y[2]*(1+y[1]-y[2])^(2)+\rthe*(1-y[2])*(1-y[1]+y[2])^(2)-2*\reta*y[2]*(1+y[1]-y[2])*(1-y[1]+y[2])}
\begin{center}
\begin{pspicture}(-0.05,-0.05)(1.05,1.05)
\psplotDiffEqn{0}{20}{0.8 0.2}{\Func}
\end{pspicture}
\end{center}
\end{document}
THE FILE LIST
*File List*
minimal.cls 2001/05/25 Standard LaTeX minimal class
amsmath.sty 2000/07/18 v2.13 AMS math features
amstext.sty 2000/06/29 v2.01
amsgen.sty 1999/11/30 v2.0
amsbsy.sty 1999/11/29 v1.2d
amsopn.sty 1999/12/14 v2.01 operator names
pstricks.sty 2008/11/26 v0.40 LaTeX wrapper for `PSTricks' (RN,HV)
pstricks.tex 2010/04/26 v2.05 `PSTricks' (tvz,hv)
pst-xkey.tex 2005/11/25 v1.6 PSTricks specialization of xkeyval (HA)
xkeyval.sty 2008/08/13 v2.6a package option processing (HA)
xkeyval.tex 2008/08/13 v2.6a key=value parser (HA)
pst-fp.tex 2010/04/26 v2.05 `PST-fp' (hv)
xcolor.sty 2007/01/21 v2.11 LaTeX color extensions (UK)
color.cfg 2007/01/18 v1.5 color configuration of teTeX/TeXLive
dvips.def 1999/02/16 v3.0i Driver-dependant file (DPC,SPQR)
pstricks-add.sty 2010/02/11 v. 0.14 package wrapper for
pstricks-add.tex (hv
)
pstricks-add.tex 2010/03/12 v3.39 `PSTricks-add' (hv)
***********
THE ERROR REPORT:
AFPL Ghostscript 8.14 (2004-02-20)
Copyright (C) 2004 artofcode LLC, Benicia, CA. All rights reserved.
This software comes with NO WARRANTY: see the file PUBLIC for details.
Displaying DSC file C:/Documents and Settings/Todos/My
Documents/Biologia/Malaria/Biomat2010/TransVector.ps
Displaying page 1
Error: /rangecheck in --get--
Operand stack:
0.0 --nostringval-- k0 1.0 1.0 --nostringval-- 2
Execution stack:
%interp_exit .runexec2 --nostringval-- --nostringval--
--nostringval-- 2 %stopped_push --nostringval--
--nostringval-- false 1 %stopped_push 1 3 %oparray_pop 1
3 %oparray_pop 1 3 %oparray_pop 1 3 %oparray_pop
.runexec2 --nostringval-- --nostringval-- --nostringval-- 2
%stopped_push --nostringval-- --nostringval-- --nostringval--
--nostringval-- --nostringval--
Dictionary stack:
--dict:1120/1686(ro)(G)-- --dict:0/20(G)-- --dict:116/200(L)--
--dict:178/300(L)-- --dict:37/200(L)-- --dict:125/200(L)--
--dict:13/20(L)--
Current allocation mode is local
Last OS error: No such file or directory
--- Begin offending input ---
aload pop k0 { 2 div } forall addvect ylength array astore
/y exch def x dup dx 2 div add /x exch def /k1 Func { dx mul } forall
ylength array astore def exch dup aload pop k1 { 2 div } forall addvect
y astore pop /k2 Func { dx mul } forall ylength array astore def dup
aload pop k2 aload pop addvect y astore pop exch dup dx add /x exch
def /k3 Func { dx mul } forall ylength array astore def /x exch def
dup aload pop k0 aload pop k1 aload pop k2 aload pop addvect 2 mulvect
addvect k3 aload pop addvect 6 divvect addvect y astore pop } def /xy
{ y 0 get 227.62204 mul y RK 1 get 227.62204 mul } def /ps at Exit false
def xy moveto 199 1 sub { /x x dx add def xy -1.0 0 gt { dup abs
-1.0 227.62204 mul gt { moveto }{ L } ifelse } { L } ifelse } ps at Exit
{ exit } if repeat ps at Exit not { /x x1 def xy -1.0 0 gt { dup abs
-1.0 227.62204 mul gt { moveto }{ L } ifelse } { L } ifelse } if
gsave 0.5 SLW 0 0 1 setrgbcolor 1. .setopacityalpha 0 setlinecap
stroke grestore end
@endspecial 0
TeXcolorgray 0 TeXcolorgra
--- End offending input ---
file offset = 181924
gsapi_run_string_continue returns -101
More information about the PSTricks
mailing list