function obj = Obbl(T,C,N,a,b,sigma,r0)

[n m] = size(C);

for i = 1:n
  B(i) = exp( - T(i) .* Vasicek(T(i),a,b,sigma,r0));
end 

C(n) = C(n) + N;

obj = C' * B';

endfunction