N=10;
i=0:N;
x=-1+2/N*i;
fx=1./（1+25*（x.^2））;
xx=-1:0.001:1;
ffx=zeros（1length（xx））;
for?i=1:N
????ffx（find（xx==x（i）））=fx（i）;
end
h=（max（x）-min（x））/N;
H=4*h*diag（ones（N1））+h*diag（ones（N-11）1）+h*diag（ones（N-11）-1）;
H=H（2:N2:N）;
for?i=1:N
????fi（i）=（fx（i+1）-fx（i））/h;
end
D=zeros（N-11）;
for?i=1:N-1
????D（i）=fi（i+1）-fi（i）;
end
ddfx=D/h;
D（1）=D（1）-h*ddfx（1）/6;
D（N-1）=D（N-1）-h*ddfx（N-1）/6;
D=6*D;
M=mathanalysis202（HD）;
M=[ddfx（1）
????M
????ddfx（N-1）];

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????989??2008-01-10?22:02??mathanalysis303.m