Watch this Author's files
How do you deal with non-close contour? e.g., >> [X,Y] = meshgrid(-10:10);Z=X.^2;
Excellent work!
Contact us at files@mathworks.com