function [cross]=cross_2d(a,b)
%
% cross product of vectors in a plane
%
cross =  a(1)*b(2) - a(2)*b(1);
