数值分析第四章选做题
从习题2,7任选一题做,求解过程如下:
习题2.7第3小题
(1)用作图法确定初始值x01和x02。
在命令窗口输入以下程序:
>> x=-2:0.00001:2;
>>e=exp(1);
>>y=x.*exp(x)-1;
>> plot(x,y),grid,
>> gtext('y=x.*exp(x)-1')
图1
运行后画出函数的图形(见图1),可见在区间(0,1)内有一实根,取初值x01 0.5,x02 0.6。
(2)建立名为fnq.m的M文件
function y=fnq(x)
e=exp(1);
y=x.*exp(x)-1;
(3)保存名为gexian.m的M文件;
(4)在MATLAB命令窗口输入以下程序:
[k,piancha,xdpiancha,xk,yk]=gexian(0.5,0.6,1e-5,1e-5,100)