Bifurcation theory studies the behavior of multiple solutions of nonlinear (differential) equations when parameters in these equations are varied, and describes how the number and type of these solutions change. It is a domain of applied mathematics which uses concepts from such diverse fields as functional analysis, group representations, ideal theory and many others. For real, e.g. physically motivated problems, the calculations necessary to determine even the simplest bifurcations become excessively complicated. Therefore, a project to build a package “bifurcation and singularity theory” in computer algebra is presented. Specifically, Gröbner bases are used to determine the codimension of a singularity, thereby extending the Buchberger Algorithm to modules. Also, a program in SMP is described, which permits determining whether a given function g is contact equivalent to a polynomial normal form h for one dimensional bifurcation problems up to codimension three.