1## [xx, yy, z] = peaks() 2## 3## return the "well" known MatLab (R) peaks function 4## evaluated in the [-3,3] x,y range. 5 6function [xx, yy, z] = peaks() 7 8 xx = yy = -3:0.1:3; 9 [x,y] = meshgrid (xx,yy); 10 11 z = 3 * (1-x).^2 .* exp(-(x.^2) - (y+1).^2) - ... 12 10 * (x/5 - x.^3 - y.^5) .* exp(-x.^2 - y.^2)- ... 13 1/3 * exp(-(x+1).^2 - y.^2); 14 15 ## z(z<-1) = -1; ## Jungfraujoch/Interlaken 16 17endfunction 18