We present the status of our study on the equation of state in 2+1 flavor QCD with non-perturbatively improved Wilson quarks coupled with the RG improved glue. We apply the T-integration method to non-perturbatively calculate the equation of state by the fixed-scale approach.