The free transverse vibrations of thin nonhomogeneous rectangular plates of variable thickness have been studied using boundary characteristic orthogonal polynomials in the Rayleigh-Ritz method. Gram-Schmidt process has been used to generate these orthogonal polynomials in two variables. The thickness variation is bidirectional and is the cartesian product of linear variations along two concurrent edges of the plate. The nonhomogeneity of the plate is assumed to arise due to linear variations in Young's modulus and density of the plate material with the in-plane coordinates. Numerical results have been computed for four different combinations of clamped, simply supported and free edges. Effect of the nonhomogeneity and thickness variation with varying values of aspect ratio on the natural frequencies of vibration is illustrated for the first three modes of vibration. Three dimensional mode shapes for all the four boundary conditions have been presented. A comparison of results with those available in the literature has been made.