数值解用vpasolve,格式[x1,s3,x3,x4]=vpasolve(eqns,vars,int_guess),分别是方程组，变量，初始值，具体百度一下用法就知道啦
使用的是迭代的算法，用到的一般是数值分析中的计算方法。例如可以使用牛顿迭代法 等等。