Computing natural frequencies and mode shapes of a non-uniform circular membrane