#!/bin/sh

set -e

# initial run forward, a HB is detected
rauto roess alf

# switch to periodic solutions
rrun=alf
rauto roess po1 $rrun -ICP 1 11 -IRS `splabs $rrun HB` \
	-DS 0.1 -DSMIN 0.05 -DSMAX 0.5 -ISW -1 -IPS 2 -NPR 1

# compute cascade of period-doublings
for (( i=2 ; $i <= ${1:-3} ; ++i )) do
	run=po$i
	rrun=po$(( $i - 1 ))
	h0=$(printac $rrun DS )
	hmin=$(printac $rrun DSMIN )
	hmax=$(printac $rrun DSMAX )
	h0=$( dc -e "8 k $h0 4.67 / p" )
	hmin=$( dc -e "8 k $hmin 4.67 / p" )
	hmax=$( dc -e "8 k $hmax 4.67 / p" )
	NTST=$(( 2*$(printac $rrun NTST) ))
	echo "NTST=$NTST, DS=$h0, DSMIN=$hmin, DSMAX=$hmax"
	rauto roess $run $rrun -ICP 1 11 -IRS `splabs $rrun PD` \
		-DS $h0 -DSMIN $hmin -DSMAX $hmax -NTST $NTST \
		-ISW -1 -IPS 2 -NPR 5
done

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

echo ""
echo "plot family of periodic solutions with:"
echo 'plotsol -i 0:2:3 po1'
echo "or"
echo 'plotsol -i -bl 2 0:2:3 po1'
