#***********************************************************************
#
#       use with test problem Square 
#
#       nonlinear eigenvalue problem
#
#***********************************************************************
#
#
#       initialize coefficient functions and create mesh
#
u c0=0, icont=15
t iadapt=4, irefn=8
#
#       follow trivial branch, locate first eigenvalue
#
s rltrgt=0, ispd=1, mxcg=30
s rltrgt=10
s itask=1, rltrgt=22
w
#
#       follow bifurcating branch for first eigenvalue
#
s itask=2
s rtrgt=.01
s rltrgt=25
s rltrgt=50
s rltrgt=100
s rltrgt=150
s rltrgt=300
s rltrgt=500
#
#       refine the mesh
#
t iadapt=1, nvtrgt=1000
s itask=7
i inplsw=5, lines=0
f nx=-1, ny=-1, ifun=0, lines=2
#
#       follow trivial branch, locate second eigenvalue
#
r
s rltrgt=30
s rltrgt=50
s itask=1, rltrgt=60
w
#
#       follow bifurcating branch for second eigenvalue
#
s itask=2
s rtrgt=.01
s rltrgt=60
s rltrgt=110
s rltrgt=150
s rltrgt=300
s rltrgt=500
#
#       refine the mesh
#
t iadapt=1, nvtrgt=1000
s itask=7
i inplsw=5, lines=0
f nx=-1, ny=-1, ifun=0, lines=2
#
#       follow trivial branch, locate third eigenvalue
#
r
s rltrgt=70
s rltrgt=80
s itask=1, rltrgt=85
w
#
#       follow bifurcating branch for third eigenvalue
#
s itask=2
s rtrgt=.01
s rltrgt=100
s rltrgt=150
s rltrgt=200
s rltrgt=300
s rltrgt=400
s rltrgt=500
#
#       refine the mesh
#
t iadapt=1, nvtrgt=1000
s itask=7
i inplsw=5, lines=0
f nx=-1, ny=-1, ifun=0, lines=2
#
#       follow trivial branch, locate fourth eigenvalue
#
r
s rltrgt=100
s rltrgt=101
s itask=1, rltrgt=103
w
#
#       follow bifurcating branch for fourth eigenvalue
#
s itask=2
s rtrgt=.005
s rltrgt=105
s rltrgt=125
s rltrgt=150
s rltrgt=200
s rltrgt=300
s rltrgt=400
s rltrgt=500
#
#       refine the mesh
#
t iadapt=1, nvtrgt=1000
s itask=7
i inplsw=5, lines=0
f nx=-1, ny=-1, ifun=0, lines=2
#
g igrsw=3
g igrsw=0
g igrsw=1
g igrsw=2
#q