function Mv=Mfun1(v,A,nu)
  n=length(v)/3;
  w=A\v(n+1:2*n);
  z=(w+v(2*n+1:3*n))./nu;
  Mv=[w;A\(v(1:n)+z);z];
end