I have problem with quodratic tern in one equation of my mod file. The equation is as follows:
mud=Omega*((1-taustar)R( 1-tau_adj*( (Dstar-steady_state(Dstar))^2)/ (1-tauB) )-ebar(+1)/ebarRstar);
tau_adj is a parameter, and Dstar is
I found that when I set the solution order equals 1, there is no difference when I introduce the quadratic term:
But when I set the solution order equals 2, the introduced quadratic term has an effect, i.e. , the result with tau_adj>0 is different with that of tau_adj=0.
Where is the problem?