The tensor optimized Fermi sphere (TOFS) method is applied first for the study of the property of nuclear matter using the Argonne V4'
NN potential. In the TOFS method, the correlated nuclear matter wave function is taken to be a power-series-type of the correlation function
F, where
F can induce central, spin-isospin, tensor, etc.~correlations. This expression has been ensured by a cluster expansion theory. In the TOFS calculation, we take into account the contributions from all the many body terms arising from the product of the nuclear-matter Hamiltonian
H and
F. It is found that the density dependence of the energy per particle in nuclear matter is reasonably reproduced, in comparison with other methods such as the Brueckner-Hartree-Fock (BHF) approach.