x=-3:14;

  i = 2;
  x1=-3:0.05:14;
  n1 = length(x1);
  for i1=1:n1,
     y1(i1)=bspline(x,i,x1(i1));
     y2(i1) = bspline(x,i+1,x1(i1));
     y3(i1) = bspline(x,i+3,x1(i1));
     y4(i1) = bspline(x,10,x1(i1));
  end

  plot(x1,y1,x1,y2,x1,y3,x1,y4) 
  title('several B-spline')
  xlabel('x')