b=[0,0,1,1,1/2,1/2;0,1,1,0,0,1/2;
b=[0,0,1,1,1/2,1/2;0,1,1,0,0,1/2];
plot(b(1,:),b(2,:),'o')
plot(b(1,:),b(2,:),'o',BScurve(b,1,t))
??? Error using ==> plot
Invalid line property: 'o'.

t=linspace(0,1);
S=BScurve(b,1,t);
plot(b(1,:),b(2,:),'o',S(1,:),S(2,:),'x')
S=BScurve(b,2,t);
plot(b(1,:),b(2,:),'o',S(1,:),S(2,:),'x')
S=BScurve(b,3,t);
plot(b(1,:),b(2,:),'o',S(1,:),S(2,:),'x')
S=BScurve(b,4,t);
plot(b(1,:),b(2,:),'o',S(1,:),S(2,:),'x')
S=BScurve(b,5,t);
plot(b(1,:),b(2,:),'o',S(1,:),S(2,:),'x')
S=BScurve(b,6,t);
plot(b(1,:),b(2,:),'o',S(1,:),S(2,:),'x')
y=decasteljau(b,t);
plot(b(1,:),b(2,:),'o',S(1,:),S(2,:),'x',y{1}(1,:),y{1}(2,:))
b1=[0,0,0.8,1,1/2,1/2;0,1,0.8,0,0,1/2];
S1=BScurve(b1,6,t);
plot(S(1,:),S(2,:),S1(1,:),S1(2,:)),legend('S','S1')
S1=BScurve(b1,3,t);
S=BScurve(b,3,t);
plot(S(1,:),S(2,:),S1(1,:),S1(2,:)),legend('S','S1')
S1=BScurve(b1,2,t);
S=BScurve(b,2,t);
plot(S(1,:),S(2,:),S1(1,:),S1(2,:)),legend('S','S1')
% Modify expression to add input arguments.
% Example:
%   a = [1 2 3; 4 5 6]; 
%   foo(a);

font
% Modify expression to add input arguments.
% Example:
%   a = [1 2 3; 4 5 6]; 
%   foo(a);

font
% Modify expression to add input arguments.
% Example:
%   a = [1 2 3; 4 5 6]; 
%   foo(a);

periodicfont
% Modify expression to add input arguments.
% Example:
%   a = [1 2 3; 4 5 6]; 
%   foo(a);

periodicfont3
% Modify expression to add input arguments.
% Example:
%   a = [1 2 3; 4 5 6]; 
%   foo(a);

periodicfont3
% Modify expression to add input arguments.
% Example:
%   a = [1 2 3; 4 5 6]; 
%   foo(a);

periodicfont3
% Modify expression to add input arguments.
% Example:
%   a = [1 2 3; 4 5 6]; 
%   foo(a);

periodicfont3
% Modify expression to add input arguments.
% Example:
%   a = [1 2 3; 4 5 6]; 
%   foo(a);

periodicfont3
% Modify expression to add input arguments.
% Example:
%   a = [1 2 3; 4 5 6]; 
%   foo(a);

periodicfont3
% Modify expression to add input arguments.
% Example:
%   a = [1 2 3; 4 5 6]; 
%   foo(a);

periodicfont3
% Modify expression to add input arguments.
% Example:
%   a = [1 2 3; 4 5 6]; 
%   foo(a);

periodicfont3
% Modify expression to add input arguments.
% Example:
%   a = [1 2 3; 4 5 6]; 
%   foo(a);

periodicfont3
% Modify expression to add input arguments.
% Example:
%   a = [1 2 3; 4 5 6]; 
%   foo(a);

font
% Modify expression to add input arguments.
% Example:
%   a = [1 2 3; 4 5 6]; 
%   foo(a);

font
% Modify expression to add input arguments.
% Example:
%   a = [1 2 3; 4 5 6]; 
%   foo(a);

font
% Modify expression to add input arguments.
% Example:
%   a = [1 2 3; 4 5 6]; 
%   foo(a);

font
% Modify expression to add input arguments.
% Example:
%   a = [1 2 3; 4 5 6]; 
%   foo(a);

font
% Modify expression to add input arguments.
% Example:
%   a = [1 2 3; 4 5 6]; 
%   foo(a);

font
% Modify expression to add input arguments.
% Example:
%   a = [1 2 3; 4 5 6]; 
%   foo(a);

font
% Modify expression to add input arguments.
% Example:
%   a = [1 2 3; 4 5 6]; 
%   foo(a);

font
% Modify expression to add input arguments.
% Example:
%   a = [1 2 3; 4 5 6]; 
%   foo(a);

font
% Modify expression to add input arguments.
% Example:
%   a = [1 2 3; 4 5 6]; 
%   foo(a);

font
size(b)

ans =

     2    29

% Modify expression to add input arguments.
% Example:
%   a = [1 2 3; 4 5 6]; 
%   foo(a);

font
??? Operation terminated by user during ==>
BSplineEval

In ==> <a href="matlab: opentoline('/home/accounts/personale/clrmrc90/aa1213/advanced_numerical_analysis/BSplineEval.m',10,0)">BSplineEval at 10</a>
    yy =
    yy+(t(j+k)-xx)/(t(j+k)-t(j+1)).*BSplineEval(j+1,k-1,t,xx);
    
In ==> <a href="matlab: opentoline('/home/accounts/personale/clrmrc90/aa1213/advanced_numerical_analysis/BSplineEval.m',10,0)">BSplineEval at 10</a>
    yy =
    yy+(t(j+k)-xx)/(t(j+k)-t(j+1)).*BSplineEval(j+1,k-1,t,xx);
    
In ==> <a href="matlab: opentoline('/home/accounts/personale/clrmrc90/aa1213/advanced_numerical_analysis/BSplineEval.m',7,0)">BSplineEval at 7</a>
    yy =
    (xx-t(j))/(t(j+k-1)-t(j)).*BSplineEval(j,k-1,t,xx);
    
In ==> <a href="matlab: opentoline('/home/accounts/personale/clrmrc90/aa1213/advanced_numerical_analysis/BSplineEval.m',10,0)">BSplineEval at 10</a>
    yy =
    yy+(t(j+k)-xx)/(t(j+k)-t(j+1)).*BSplineEval(j+1,k-1,t,xx);
    
In ==> <a href="matlab: opentoline('/home/accounts/personale/clrmrc90/aa1213/advanced_numerical_analysis/BSplineEval.m',10,0)">BSplineEval at 10</a>
    yy =
    yy+(t(j+k)-xx)/(t(j+k)-t(j+1)).*BSplineEval(j+1,k-1,t,xx);
    
In ==> <a href="matlab: opentoline('/home/accounts/personale/clrmrc90/aa1213/advanced_numerical_analysis/BSplineEval.m',7,0)">BSplineEval at 7</a>
    yy =
    (xx-t(j))/(t(j+k-1)-t(j)).*BSplineEval(j,k-1,t,xx);
    
In ==> <a href="matlab: opentoline('/home/accounts/personale/clrmrc90/aa1213/advanced_numerical_analysis/BSplineEval.m',10,0)">BSplineEval at 10</a>
    yy =
    yy+(t(j+k)-xx)/(t(j+k)-t(j+1)).*BSplineEval(j+1,k-1,t,xx);
    
In ==> <a href="matlab: opentoline('/home/accounts/personale/clrmrc90/aa1213/advanced_numerical_analysis/BSplineEval.m',10,0)">BSplineEval at 10</a>
    yy =
    yy+(t(j+k)-xx)/(t(j+k)-t(j+1)).*BSplineEval(j+1,k-1,t,xx);
    
In ==> <a href="matlab: opentoline('/home/accounts/personale/clrmrc90/aa1213/advanced_numerical_analysis/BSplineEval.m',7,0)">BSplineEval at 7</a>
    yy =
    (xx-t(j))/(t(j+k-1)-t(j)).*BSplineEval(j,k-1,t,xx);
    
In ==> <a href="matlab: opentoline('/home/accounts/personale/clrmrc90/aa1213/advanced_numerical_analysis/BSplineEval.m',10,0)">BSplineEval at 10</a>
    yy =
    yy+(t(j+k)-xx)/(t(j+k)-t(j+1)).*BSplineEval(j+1,k-1,t,xx);
    
In ==> <a href="matlab: opentoline('/home/accounts/personale/clrmrc90/aa1213/advanced_numerical_analysis/BSplineEval.m',10,0)">BSplineEval at 10</a>
    yy =
    yy+(t(j+k)-xx)/(t(j+k)-t(j+1)).*BSplineEval(j+1,k-1,t,xx);
    
In ==> <a href="matlab: opentoline('/home/accounts/personale/clrmrc90/aa1213/advanced_numerical_analysis/BSplineEval.m',10,0)">BSplineEval at 10</a>
    yy =
    yy+(t(j+k)-xx)/(t(j+k)-t(j+1)).*BSplineEval(j+1,k-1,t,xx);
    
In ==> <a href="matlab: opentoline('/home/accounts/personale/clrmrc90/aa1213/advanced_numerical_analysis/BSplineEval.m',7,0)">BSplineEval at 7</a>
    yy =
    (xx-t(j))/(t(j+k-1)-t(j)).*BSplineEval(j,k-1,t,xx);
    
In ==> <a href="matlab: opentoline('/home/accounts/personale/clrmrc90/aa1213/advanced_numerical_analysis/BSplineEval.m',10,0)">BSplineEval at 10</a>
    yy =
    yy+(t(j+k)-xx)/(t(j+k)-t(j+1)).*BSplineEval(j+1,k-1,t,xx);
    
In ==> <a href="matlab: opentoline('/home/accounts/personale/clrmrc90/aa1213/advanced_numerical_analysis/BSplineEval.m',7,0)">BSplineEval at 7</a>
    yy =
    (xx-t(j))/(t(j+k-1)-t(j)).*BSplineEval(j,k-1,t,xx);
    
In ==> <a href="matlab: opentoline('/home/accounts/personale/clrmrc90/aa1213/advanced_numerical_analysis/BSplineEval.m',7,0)">BSplineEval at 7</a>
    yy =
    (xx-t(j))/(t(j+k-1)-t(j)).*BSplineEval(j,k-1,t,xx);
    
In ==> <a href="matlab: opentoline('/home/accounts/personale/clrmrc90/aa1213/advanced_numerical_analysis/BSplineEval.m',10,0)">BSplineEval at 10</a>
    yy =
    yy+(t(j+k)-xx)/(t(j+k)-t(j+1)).*BSplineEval(j+1,k-1,t,xx);
    
In ==> <a href="matlab: opentoline('/home/accounts/personale/clrmrc90/aa1213/advanced_numerical_analysis/BSplineEval.m',7,0)">BSplineEval at 7</a>
    yy =
    (xx-t(j))/(t(j+k-1)-t(j)).*BSplineEval(j,k-1,t,xx);
    
In ==> <a href="matlab: opentoline('/home/accounts/personale/clrmrc90/aa1213/advanced_numerical_analysis/B.m',10,0)">B at 10</a>
ret = BSplineEval(i,k,knots,t);

In ==> <a href="matlab: opentoline('/home/accounts/personale/clrmrc90/aa1213/advanced_numerical_analysis/BScurve.m',9,0)">BScurve at 9</a>
  s = s+diag(b(:,j))*repmat(B(j,n,k,t),size(b,1),1);

In ==> <a href="matlab: opentoline('/home/accounts/personale/clrmrc90/aa1213/advanced_numerical_analysis/font.m',7,0)">font at 7</a>
s = BScurve(b,k,t);

b=rand(2,10);
S=BSspline(b,3,t);
??? Undefined function or method 'BSspline' for
input arguments of type 'double'.

S=BScurve(b,3,t);
size(t)

ans =

           1        1001

plot(S(1,:),S(2,:))
y=decasteljau(b,t);
plot(S(1,:),S(2,:),y{1}(1,:),y{1}(2,:))
b1=b;
b1(:,5)=b1(:,5)/2;
b

b =

  Columns 1 through 5

    0.8147    0.1270    0.6324    0.2785    0.9575
    0.9058    0.9134    0.0975    0.5469    0.9649

  Columns 6 through 10

    0.1576    0.9572    0.8003    0.4218    0.7922
    0.9706    0.4854    0.1419    0.9157    0.9595

b1

b1 =

  Columns 1 through 5

    0.8147    0.1270    0.6324    0.2785    0.4788
    0.9058    0.9134    0.0975    0.5469    0.4824

  Columns 6 through 10

    0.1576    0.9572    0.8003    0.4218    0.7922
    0.9706    0.4854    0.1419    0.9157    0.9595

S1=BScurve(b1,3,t);
y1=decasteljau(b1,t);
plot(S(1,:),S(2,:),y{1}(1,:),y{1}(2,:),S1(1,:),S1(2,:),y1{1}(1,:),y1{1}(2,:))
legend('BS b','Bezier b','BS b1','Bezier b1')
plot(S(1,:),S(2,:),y{1}(1,:),y{1}(2,:),S1(1,:),S1(2,:),'x',y1{1}(1,:),y1{1}(2,:))
legend('BS b','Bezier b','BS b1','Bezier b1')
plot(S(1,:),S(2,:),y{1}(1,:),y{1}(2,:),S1(1,:),S1(2,:),'x',y1{1}(1,:),y1{1}(2,:),'x')
legend('BS b','Bezier b','BS b1','Bezier b1')
diary off