A formulation, using the double time Green’s function method, is presented for the anisotropic spin-one ferromagnetic system in the presence of Ising part of the biquadratic exchange and an expression is given for the critical temperatureT_c in terms of two parameters. One is associated with the uniaxial anisotropy (treated exactly) and the other one corresponds to the biquadratic interaction (treated in RPA). The bilinear Heisenberg exchange term is treated in Callen’s approximation emphasizing the role of intersite correlation for small parameters. The prediction thatT_c increases with the biquadratic exchange is found to be consistent with the recent calculation based on high temperature series expansion.