This paper presents a design method for handling magnetic force density distributions in 2-D linear magnetic systems. In order to control not the global force but the distribution of local force density, the discrete sensitivity analysis is combined with the magnetic charge method. For calculation of the magnetic force density, we use the magnetic charge method that is valid from the viewpoint of mechanical deformation. The design sensitivity is obtained through the direct differentiation of the finite element system matrices. Two models of a levitated rod and a cantilever electromagnet are designed to produce the desired force density distribution. Their results show the validity of the proposed method.