A methodology for adaptive finite element analysis: Towards an integrated computational environment. We accomplish the general implementation of the superconvergent patch recovery (SPR) and the recently proposed recovery-by-equilibrium in patches (REP). Both SPR and REP are compared and used for error estimation and for guiding the adaptive remeshing process. Moreover, the SPR is extended for calculating sensitivity quantities of first- and higher-orders. We accomplish the mesh (re-)generation process by means of modern methods, combining quadtree and Delaunay triangulation techniques. Surface mesh generation in arbitrary domains is performed automatically (i.e. with no user intervention) during the self-adaptive analysis, using quadrilateral or triangular elements. These ideas are implemented in the finite element system technology in adaptivity (FESTA) software. The effectiveness and versatility of FESTA are demonstrated by numerical examples.