function x = isodd (number)
% ISODD(x)
% returns 1 if the x is odd, 0 if it is even.

% Divide input by 2:

a = number./2;
whole = floor(a);
part = a-whole;
sz_rs=size(number);
part_rs=reshape(part,[],1);
x_rs=ones(size(part_rs))*NaN;
for ii=1:length(part_rs)
    if part_rs(ii) > 0;
        x_rs(ii) = 1;
    else
        x_rs(ii) = 0;
    end
end
x=logical(reshape(x_rs,sz_rs));