%《MATLAB及其在理工课程中的应用指南》(十一五规划版)第五章例5-1-8程序
% 斜率计算和切线绘制
% 西安电子科技大学出版社出版 陈怀琛编著 2007年5月
%
clear,clf,
x=linspace(-1/2,3,100);
subplot(1,2,1),plot(x,x.^3+2*x),grid on,hold on
% 在一个点上,取不同步长求斜率,画出的多根切线的程序。
h=[1,0.1,0.001];
for k=1:3
x1=[1,1+h(k)];f= x1.^3+2*x1;
q=diff(f)/diff(x1),
plot(x,f(1)+q*(x-x1(1)),':'),
end
% 在多个点上,取同样步长h=0.1求斜率,画出的多根切线的程序。
for k=0:2 x1=[k,k+0.1];f= x1.^3+2*x1;
q=diff(f)/diff(x1),
subplot(1,2,2),
plot(x,x.^3+2*x),grid on,hold on
plot(x,f(1)+q*(x-x1(1)),':r'),
end
pause
syms x
y=x^3+2*x,
Dy=diff(y),
Dy1=subs(Dy,x,[0,1,2])