Consider a bifurcation problem, namely, its bifurcation equation. There is a diffeomorphism $\Phi$ linking the actual solution set with an unfolded normal form of the bifurcation equation. The differential $D\Phi \left(0\right)$ of this diffeomorphism is a valuable information for a numerical analysis of the imperfect bifurcation. The aim of this paper is to construct algorithms for a computation of $D\Phi \left(0\right)$. Singularity classes containing bifurcation points with $\mathrm{c}odim\le 3$, $\mathrm{c}orank=1$ are considered.

The Recursive Projection Method is a technique for continuation of both the steady states and the dominant invariant subspaces. In this paper a modified version of the RPM called projected RPM is proposed. The modification underlines the stabilization effect. In order to improve the poor update of the unstable invariant subspace we have applied subspace iterations preconditioned by Cayley transform. A statement concerning the local convergence of the resulting method is proved. Results of numerical...