y1'=-.04*y1+1e4*y2*y3
y2'=.04*y1-1e4*y2*y3-3e7*y2^2
y3'=3e7*y2^2
init y1=1,y2=0,y3=0
aux c=y1+y2+y3
done