# lampul.ode -- boundary value problem
r'=L*r*u
u'=(phi*phi-u*u-lambda(r)/d)*L
phi'=((omega-q*r*r)/d-2*u*phi)*L
omega'=0

# initial data
r(0)=0.9031584151492208
u(0)=0
phi(0)=0
omega(0)=.575192015732142
# parameters
par q=1,d=.1,a=.25,L=2.5
lambda(r)=(1-r*r)*(r*r-a*a)
# boundary conditions
bndry u
bndry phi
bndry u'
bndry phi'
@ dt=.01,total=1.001
done