In this paper, a treatment for determining the eigenfunctions and eigenvalues for the vibration of a circular membrane with fixed boundary and density varying as exp[− εγ2] is made by the method of perturbation, and the expressions for perturbed eigenvalues and eigenfunctions are obtained up to second‐order terms. The first‐order perturbed term of the eigenvalue is effected by numerical integration.