Dear all

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)/ebar*Rstar);
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:
tau_adj*( (Dstar-steady_state(Dstar))^2)

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?

Thank you.