f=x^5-5*x-1; solve(f) 例8分别建立命令式M文件和函数式M文件求向量的平均值和标准差 解:命令式函数文件 在文件编辑窗口输入如下命令: x=rand (20)i stdev=sqrt(sum(x. 2)/n-mean. 2) 随便以一个文件名命名即可,例如保存为 programme. m,在命令窗口输入 programme 即可得到结果。 函数式函数文件 在文件编辑窗口输入如下命令: function [mean, stdev]=s tat(x) mean=sum(x)/n stdevesgrt(sum(x. 2)/n-mean. 2) 保存为文件 stat m,在命令窗口输入 x=rand(20)i [me an, stdev]=stat(x) 即可得到结 例9画出f(x)=e2在x=0(P(0,1))处的切线及若干条割线,观察割线的变化趋势。 解:在曲线y=e2上另取一点M(h,e"),则PM的方程是 x+1,取h=3,2,10.1,0.01,分别作出几条割线。 h 然后作出y=e在x=0处的切线y=x+1。 h=[3,2,1,0.1,0.01]; a=(exp(h)-1)./h 1:0.1:3; plot(x, exp(x),'r) hold on for i=l: 5 plot(h(i), exp(h(i)),'r') plot(x, a(i)*x+l axis square plot(x, x+l,'r') 从图上看出割线越来越接近切线f=x^5-5*x-1; solve(f) 例8 分别建立命令式M文件和函数式M文件求向量的平均值和标准差。 解：命令式函数文件： 在文件编辑窗口输入如下命令： x=rand(20); n=length(x); mean=sum(x)/n stdev=sqrt(sum(x.^2)/n-mean.^2) 随便以一个文件名命名即可，例如保存为programme.m，在命令窗口输入programme 即可得到结果。 函数式函数文件： 在文件编辑窗口输入如下命令： function [mean,stdev]=stat(x) n=length(x); mean=sum(x)/n stdev=sqrt(sum(x.^2)/n-mean.^2) 保存为文件stat.m，在命令窗口输入： x=rand(20); [mean,stdev]=stat(x) 即可得到结果。 例 9 画出 x f (x) = e 在 x = 0 （ P(0,1) ）处的切线及若干条割线，观察割线的变化趋势。 解：在曲线 x y = e 上另取一点 ( , ) h M h e ，则 PM 的方程是： 1 1 + − = x h e y h ，取 h = 3,2,1,0.1,0.01 ，分别作出几条割线。 然后作出 x y = e 在 x = 0 处的切线 y = x +1。 h=[3,2,1,0.1,0.01]; a=(exp(h)-1)./h; x=-1:0.1:3; plot(x,exp(x),'r'); hold on for i=1:5 plot(h(i),exp(h(i)),'r.') plot(x,a(i)*x+1) end axis square plot(x,x+1,'r') 从图上看出割线越来越接近切线