function [ mdot ] = LL( m )
global gamma H lambda nu s
H1=H+nu*s;
mdot =  -gamma.*cross(m,H1)+lambda.*cross(m,cross(m,H));
end

