)+x.*exp(x);%四阶导函数plot(x,f)%画出原函数x=2;f=4*exp(x)+x.*exp(x);%计算导函数最大值h4=0.5*10^(-7)*180*16/f;h=sqrt(sqrt(abs(h4)))n=1/h;%求分段区间个数n=2*ceil(1/h)+1%选取的点数估计结果步长h及结点数n分别为h=0.13411n=47h=0.76542n=35h=0.18433n=29h=0.18546n=49积分计算的Matlab程序:formatlonggpromps={'请选择积分公式,若用复化梯形,请输入T,用复化simpson,输入S,用复化Gauss_Legendre,输入GL:'};result=inputdlg(promps,'charpt4',1,{'T'});Nb=char(result);if(Nb~='T'&Nb~='S'&Nb~='GL')errordlg('积分公式选择错误');return;endresult=inputdlg({'请输入积分式题号1-4:'},'实验4.1',1,{'1'});Nb_f=str2num(char(result));if(Nb_f<1|Nb_f>4)errordlg('没有该积分式');return;endswitchNb_fcase1fun=inline('-2./(x.^2-1)');a=2;b=3;case2fun=inline('4./(x.^2+1)');a=0;b=1;case3fun=inline('3.^x');a=0;b=1;case4fun=inline('x.*exp(x)');a=1;b=2;endif(Nb=='T')%用复化梯形公式promps={'请输入用复化梯形公式应取的步长:'};result=inputdlg(promps,'实验4.2',1,{'0.01'});
全文预览
