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