The synthesis of spatial mechanisms is formulated as a mathematical programming problem and solved using a penalty function approach. The objective function as well as the equality and inequality constraints are determined explicitly from conditions such as those required for linkage closure, mobility, and transmissibility.