Diatomic eigenvalue problem and wavepacket propagation for collinear H+H2 using non-orthognal coordinates: the Bond Order formulation