三次样条差值的matlab仿真

三次样条差值-三弯矩法的matlab仿真，用于数值分析

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