#!/usr/bin/octave
################################################################################
#zad 2 jak 1 ale dla splajny s_k t. ze s(x(j))=1 k=j albo 0 j<>k
################################################################################
#funkcje spline i ppval
a=-3;b=3;
x=a:b;
y=zeros(1,length(x));y(3)=1;
pp=spline(x,y);
z=linspace(a,b);
plot(z,ppval(pp,z),";splajn rowny 1 w 3cim wezle,0 w pozostalych;",x,y,"r+")
pause(1);

yh=[0,y,0];
pph=spline(x,yh);
plot(z,ppval(pp,z),";not-a-knot;",z,ppval(pph,z),";herm;",x,y,"r+")
pause(1);
z=z(1:10);
plot(z,ppval(pp,z),";not-a-knot;",z,ppval(pph,z),";herm;")
pause(3);

