A finite element method for first order elliptic systems in three dimensions