Springs are widely utilized in machine element. To find out stiffness of coil spring, the space beam theory using the finite element method is adopted in this paper. In three dimensional space, a space frame element is a straight bar of uniform cross section which is capable of resisting axial forces, bending moments about two principal axes in the plane of its cross section and twisting moment about its centroidal axis. The corresponding displacement degrees of freedom are twelve. The displacements of nodal points due to small increment of force are calculated by the finite element method and the calculated nodal displacements are added to coordinates of nodal points. The new stiffness matrix of the system using the new coordinates of nodal points is adopted to calculated the another increments of nodal displacements, that is, the step by step method is used in this paper. The results of the finite element method are fairly well agreed with those of various experiments. Using MATLAB program developed in this paper, spring constants can be predicted by input of few factors.