function [u,resid] = EE(M,a,mu,dx,ap,am,h,u)

R = rhs(M,a,mu,dx,ap,am,u);	
resid = norm(R);
u = update(M,h,u,R);