finite element technology

Hello everybody,



I'm a PhD student and as part of my research I need to implement a triangular finite element using the displacement/pressure mixed formulation. I'm implementing two field incompressible elasticity from the book of Zienkiewicz and Taylor "The finite element method, Vol. I". The mixed approximation form is given as:



|A*u + C*p = f1

|

|C^T*u - V*p = 0



which leads to (A + C*V^-1*C^T)*u=f1.