function c = comm(a,b)



c = a*b-b*a;