# a heteroclinic orbit
# unstable at u=1, stable at u=0
u'=up*per
up'=(-c*up-u*(1-u)*(u-a))*per
per'=0
c'=0
# unstable manifold
b (u-1)*lamp-up
b (u'-0)*lamn-up'
b sqrt((u-1)^2+up^2)-eps
b sqrt((u')^2+(up')^2)-eps 
lamn=(c-sqrt(c^2+4*a))/2
lamp=(c+sqrt(c^2+4*(1-a)))/2
par a=.5,eps=.1
init u=.918,up=-.0577,per=6.67,c=0
@ total=1.01,dt=.01
done