Post-Newtonian SPH calculations of binary neutron star coalescence: Method and first results