将下列代码以ex511bf.m文件名保存到当前目录下，就可以调用ex511bf(tvalues,w)。tvalues和w是一个具体的数值。
ex511bf.m
functionxvalues=ex511bf(tvalues,w)
%例5-1-1用的函数文件,它应该能用元素群运算。
xvalues=sqrt(3)/2*exp(-4*tvalues).*sin(w*tvalues+pi/3);
dt=0.01;t=[0:dt:1.5];w=4*sqrt(3);
y=ex511bf(t,w);Dy=diff(y)/dt;
subplot(2,1,1),plot(t,y),grid%绘制曲线并加上坐标网格
title('绘图示例'),xlabel('时间t'),ylabel('y(t)')%加标注
%求导数并绘制导数曲线,注意数组求导数后其长度减少一
subplot(2,1,2),plot(t(1:length(t)-1),Dy),grid
ylabel('Dy(t)')%加标注
