using thetan*thetaeq/tau_n not working