A novel numerical method of high-order accuracy for flow in unsaturated porous media