Nonautomonous ordinary differential equations, depending on two parameters (\mu_{1}) and (\mu_{2}), are considered in (\mathbb{R}^{n}). It is assumed that when both parameters are zero the differential equation is autonomous with a hyperbolic equilibrium and a homoclinic solution. No restriction is placed on the dimension of the phase space, (\mathbb{P}^{\prime \prime}), or on the dimension of intersection of the stable and unstable manifolds. By means of the method of Lyapunov-Schmidt a bifurcation function, (H), is constructed between two finite dimensional spaces where the zeros of (H) correspond to homoclinic solutions at nonzero parameter values. The independent variables of (H) consist of scalars (\mu_{1}, \mu_{2}, \xi) and a vector (\beta) where (\xi) is a phase angle and (\beta) corresponds to directions, other than along the original homoclinic solution, tangent 10 both the stable and unstable manifolds. When (\xi) is fixed the equation (H=0) yields, in general, several bifurcation curves through the origin in the (\mu_{1}-\mu_{2}) plane along which there exists a homoclinic solution. When (\xi) is varied these become a number of wedge-shaped regions. The theory is applied to two examples, one in (\mathbb{R}^{6}) where the invariant manifolds meet in dimension three and a second in (\mathbb{P}^{4}) where these manifolds agree. β 1995 Academic Press. Inc