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