FEM working code (only steady)