# torus bifurcation
par rn=-.9,be=.5,ga=-.6,r=.6,a3=.328578,b3=.933578
x'=( -(be+rn)*x + be*y - a3*x^3 + b3*(y-x)^3 )/r
y'=be*x - (be+ga)*y - z - b3*(y-x)**3
z'=y
done


3 1 0 1				NDIM,IPS,IRS,ILP
1   1   			NICP,(ICP(I),I=1,NICP)
25 4 3 2 1 0 0 0 		NTST,NCOL,IAD,ISP,ISW,IPLT,NBC,NINT
125 -0.9 -0.4 0 5 		NMX,RL0,RL1,A0,A1
500 10 2 8 5 3 1		NPR,MXBF,IID,ITMX,ITNW,NWTN,JAC
1e-07 1e-07 1e-05 		EPSL,EPSU,EPSS
0.001 0.0001 0.01 1 		DS,DSMIN,DSMAX,IADS
1 				NTHL,((I,THL(I)),I=1,NTHL)
11 0
0 				NTHU,((I,THU(I)),I=1,NTHU)
0 				NUZR,((I,UZR(I)),I=1,NUZR)