script to compute neutrino pz solution