u1'=u2
u2'=-r*exp(u1)
b u1
b u1'
par r=0
@ total=1
done