octave:2> S.a=1 S = { a = 1 } octave:3> S.b=[1,2] S = { a = 1 b = 1 2 } octave:4> S.a ans = 1 octave:5> S.b ans = 1 2 octave:6> x=linspace(0,10,11) x = 0 1 2 3 4 5 6 7 8 9 10 octave:7> x=linspace(0,10,11)' x = 0 1 2 3 4 5 6 7 8 9 10 octave:8> y=sin(x); octave:9> pp=spline(x,y) pp = { x = 0 1 2 3 4 5 6 7 8 9 10 P = -0.04186 -0.26124 1.14457 0.00000 -0.04186 -0.38682 0.49651 0.84147 0.14694 -0.51240 -0.40272 0.90930 0.16035 -0.07157 -0.98669 0.14112 0.03722 0.40946 -0.64880 -0.75680 -0.12340 0.52112 0.28178 -0.95892 -0.16837 0.15093 0.95384 -0.27942 -0.06404 -0.35418 0.75059 0.65699 0.11895 -0.54630 -0.14989 0.98936 0.11895 -0.18945 -0.88564 0.41212 n = 10 k = 4 d = 1 } octave:10> xx=linspace(0,10,100); octave:11> yy=ppval(pp,xx); octave:12> plot(xx,yy) octave:13> plot(xx,yy,x,y,'*') octave:14> x x = 0 1 2 3 4 5 6 7 8 9 10 octave:15> x(2:11)-x(1:10) ans = 1 1 1 1 1 1 1 1 1 1 octave:16> y(2:11)-y(1:10) ans = 0.841471 0.067826 -0.768177 -0.897923 -0.202122 0.679509 0.936402 0.332372 -0.577240 -0.956140 octave:17> pp=spline(x,y); octave:18> yy=ppval(pp,xx); octave:19> plot(xx,yy,x,y,'*') octave:20> ppn=splinenaturale(x,y); warning: time stamp for `/home/accounts/altri/caliari/aa0809/calcolo_numerico/interpolazione/splinenaturale.m' is in the future octave:21> yyn=ppval(ppn,xx); octave:22> plot(xx,yy,xx,yyn,x,y,'*') octave:23> axis([0,2]) octave:24> quit