用于求解常微分方程边值问题的求解,具体请参看程序