x = ([0,0.785,1.570,2.356,3.141,3.927,4.712,5.])
y=([0,0.707,1.000,0.707,0.000,-0.707,-1.000,-0.707,-0.000])
xi=linspace(0.2*pi,50)
yi=interp1(x,y,xi,'cubic')

Strang-5.2.5-1 (last edited 2010-12-18 10:15:02 by 10)