function fu=abu_f(u,u0,t,f)



[fu,Dfu] = per_Df(u,u0,t,f,@phase_condition);