The following MATLAB functions that compute GS-filters are available at the MATLAB File Exchange:
File IDs [4038; 5661; 4270; 6121; 7003; 5556; 6047; 6151]
This zip-folder contains two functions: ztrans1d.m and ztrans2d.m that compute the z-transform of any GS filter and plot its amplitude spectrum.
Following is a demo of these functions in action for 1D and 2D GS-filters.