微分方程（组）边值问题的变分原理及MATLAB求解