A Novel Numerical Method of High-Order Accuracy for Flow in Unsaturated Porous Media