# dae1.ode # simple DAE x'=y 0= y+x-1 solve y=1 init x=0 aux yy=y done