function b=sympck(B)
%function b=sympck(B)
% Packs lower half of symmetric matrix B into col vector b
% Symunpck.m reverses procedure

%Arthur Jutan 
%Dept of Chemical and Biochemical Enginering
%Univ of Western Ontario
%London Ontario Canada, N6A 5B9

B=tril(B); % get lower half top half = zeros
b=B(:); % long col
i=find(b==0); % get zero positions
b(i)=[]; % remove and shorten vector

