exm05122_1.m

这是本人收集的一些程序源代码

%exm05122_1.m
x=(0:0.1:1)*2*pi;y=sin(x);
pp=spline(x,y);
int_pp=fnint(pp);
der_pp=fnder(pp);
xx=(0:0.01:1)*2*pi;
err_yy=max(abs(ppval(pp,xx)-sin(xx)))
err_int=max(abs(ppval(int_pp,xx)-(1-cos(xx))))
err_der=max(abs(ppval(der_pp,xx)-cos(xx)))  
%（2）
DefiniteIntegral.bySpline=ppval(int_pp,[1,2])*[-1;1];		%	<2>
DefiniteIntegral.byTheory=(1-cos(2))-(1-cos(1));
Derivative.bySpline=fnval(der_pp,3);
Derivative.byTheory=cos(3);
Derivative.byDiference=(sin(3.01)-sin(3))/0.01;
DefiniteIntegral,Derivative  
%（3）
clf,shg
fnplt(pp,'b-');hold on
fnplt(int_pp,'m:'),fnplt(der_pp,'r--');hold off
legend('y(x)','S(x)','dy/dx')