Model nuclear matter calculations with a new fermion lowest order constrained variational method