function   s=sig(a)

%This function calculates a ``Paulion" sigma.a

global_declarations;

s= a(1)*sigx + a(2)*sigy + a(3)*sigz;