000000000000011839I
fNaNNaNNaNNaN
第2页共2页
f由上面的输出结果可以看出，所建立的迭代函数并不收敛，即方程无解。为此我们建立加速迭代函数：（此处为y1的加速迭代函数）建立M文件：第二题y1的加速迭代函数：x10x24x1510x1315x1415x122whileabsx1x2105x1x2x24x1510x1315x1415x122e
dx1输出结果：x107685其迭代过程为：第二题y1加速迭代函数：x107fork130xk14xk510xk315xk415xk22e
dt0130plottxtitley1gridx输出结果为：x07000077890768707685076850768507685076850768507685076850768507685076850768507685076850768507685076850768507685076850768507685076850768507685076850768507685
由上图可知该迭代函数是收敛的！并可得x076853．求解下列方程组
第3页共3页
f1
2x1x2ex1xx12x2e2
22x125x27x3123x1x2x1x311x102xx40x0123
2
直接使用MATLAB命令：solve和fsolve对方程组求解。（1）编写M文件如下：第三题1solve解答：x1x2solve2x1x2expx1x12x2expx2fsolve解答：fu
ctio
pQ0301xp12x1x2expx1p2x12x2expx2fsolveQ0301111输出结果：Solve：x156714329040978387299996866221036x256714329040978387299996866221036Fsolve：a
s0567105671（2）编写M文件如下：第三题2solve解答：x1x2x3solvex125x227x33123x1x2x1x311x102x2x340x10doublex1doublex2doublex3fsolve解答：fu
ctio
lQ0302xl1x125x227x3312l23x1x2x1x311x1l32x2x340x1fsolveQ03021111输出结果：Solve：a
s10e00700000
第4页共4页
f10570000000000000000i0000000000ia
s10e003000015i000015i00083964000300004000007i0004000007ia
s10e00400000010000100001i0000100001i25178000020000100002i0000100002iFsolve：a
s000000000011968【注解：以上两种解法不能同时运行，需先注释一个再运行另外一个即可】4．编写用二分法求方程根的函数M文件。第四题：fu
ctio
xEfsm
Efm为做端点，
为右端点，Ef为待解方程，x为方程根y1fevalEfmy2fevalEf
pm
2y3fevalEfpwhileabsm
108
第5页共5页
fify10xmbreakelseify20x
breakelseify30xpbreakelsey1fevalEfmy2fevalEf
pm
2y3fevalEfpify1y30
pelsempe
de
de
de
de
ddispTheequatio
rootisxm
2在Matlabr