ex4.m

重点内容是计算中的误差、函数方程求根、插值呆逼近、数值积分和微分、线性代数方程组解法、常微分方程初值问题的数值解法、数学软件。

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')