func3p2.asv

来自「数值分析最常用的四十种算法」· ASV 代码 · 共 16 行

yy=linspace(0,1000,121)
xx=linspace(0,12000,121)
for i=0:120
  y=1000*sub3p2(i*100)
  xx(i+1)=i*100
  yy(i+1)=y
end
plot(xx,yy,'-')
%figure enclosed
s=dsolve('D4y=0','y(0)=0','y(12000)=1000','Dy(0)=0','Dy(12000)=0')
f=linspace
for i=0:100:12000
f=polyval(s,i)
end
plot(s,f,'-')