A method to calculate subtraction constants application to the η → π° γγ problem