#!/bin/sh

set -e

# initial run with rfcon, numerical solution in duff.dat
rfcon duff

# continuation in omega
rauto duff omf dat -IRS 1
rauto duff omb dat -IRS 1 -DS - -NMX 300 -UZR 2 1

# compute solutions branching off in BP
rauto duff po1 omb -IRS `splabs omb BP 1 1` -ISW -1 -NTST 30 -DSMAX 0.05

# compute some period-doubled solutions
rauto duff po2 po1 -IRS `splabs po1 PD 1 1` -ISW -1 \
	-NTST 60 -DS 0.01 -DSMAX 0.02 -NMX 1000
rauto duff po3 po2 -IRS `splabs po2 PD 1 1` -ISW -1 \
	-NTST 100 -DS 0.001 -DSMAX 0.005 -NMX 4000

# show graphics
echo ""
echo "plot bifurcation diagram with:"
echo 'plotbd -i 1:3 `lsruns "om*" "po*"`'
