ex4.m

来自「重点内容是计算中的误差、函数方程求根、插值呆逼近、数值积分和微分、线性代数方程组」· M 代码 · 共 10 行

x=[0,70,130,210,337,578,776,1012,1142,1462,1841];
y1=[0,57,78,103,135,182,214,244,256,272,275];
y2=[0,0,57,78,103,135,182,214,244,256,272,275,0];
v=0:0.1:1841;
ywhat1=spline(x,y1);ywhat2=spline(x,y2);
y1=ppval(ywhat1,v);y2=ppval(ywhat2,v);
y1_50=ppval(ywhat1,50)
y2_50=ppval(ywhat2,50)
plot(v,y1,'r');hold on;
plot(v,y2,'k')