#decay.ode #the initial condition is 1 X(0)=1 #the parameters param tau=2 #the equation dX/dt=-X/tau #an auxiliary function aux rate=-X/tau param init=1 aux exact=init*exp(-t/tau) end #end of file