Scilab Function

contour - level curves on a 3D surface

### Calling Sequence

  `contour(x,y,z,nz,[theta,alpha,leg,flag,ebox,zlev])`
  `contour(x,y,z,nz,<opt_args>)`

### Parameters

- `x,y` : two real row vectors of size n1 and n2.
- `z` : real matrix of size (n1,n2), the values of the function.
- `nz` : the level values or the number of levels.

- `-` If `nz` is an integer, its value gives the number of level curves equally spaced from zmin to zmax as follows:

     z= zmin + (1:nz)*(zmax-zmin)/(nz+1)

Note that the `zmin` and `zmax` levels are not drawn (generically they are reduced to points) but they can be added with

     [im,jm] = find(z == zmin);     // or zmax
     plot2d(x(im)',y(jm)',-9,"000")

- `-` If `nz` is a vector, `nz(i)` gives the value of the ith level curve. Note that it can be useful in order to see `zmin` and `zmax` level curves to add an epsilon tolerance: `nz=[ zmin+%eps,..,zmax-%eps]`.

- `<opt_args>` : a sequence of statements `key1=value1, key2=value2`, ... where keys may be `theta`,`alpha`,`leg`,`flag`, `ebox`,`zlev` (see below). In this case, the order has no special meaning.
- `theta, alpha` : real values giving in degree the spherical coordinates of the observation point.
- `leg` : string defining the captions for each axis with @ as a field separator, for example "X@Y@Z".
- `flag` : a real vector of size three `flag=[mode,type,box]`.

- `mode` : string (treatment of hidden parts).

- `mode>0` the hidden parts of the surface are removed and the surface is painted with color `mode`.
- `mode=0` the hidden parts of the surface are drawn.
- `mode<0` only the shadow of the surface is painted with color or pattern id `-mode`. Use `xset()` to see the meaning of the ids.

- `type` : an integer (scaling).

- `type=0` the plot is made using the current 3D scaling (set by a previous call to `param3d`, `plot3d`, `contour` or `plot3d1`).
- `type=1` rescales automatically 3d boxes with extreme aspect ratios, the boundaries are specified by the value of the optional argument `ebox`.
- `type=2` rescales automatically 3d boxes with extreme aspect ratios, the boundaries are computed using the given data.
- `type=3` 3d isometric with box bounds given by optional `ebox`, similarily to `type=1`
- `type=4` 3d isometric bounds derived from the data, to similarily`type=2`
- `type=5` 3d expanded isometric bounds with box bounds given by optional `ebox`, similarily to `type=1`
- `type=6` 3d expanded isometric bounds derived from the data, similarily to `type=2`

- `box` : an integer (frame around the plot).

- `box=0` nothing is drawn around the plot.
- `box=1` unimplemented (like box=0).
- `box=2` only the axes behind the surface are drawn.
- `box=3` a box surrounding the surface is drawn and captions are added.
- `box=4` a box surrounding the surface is drawn, captions and axes are added.

- `ebox` : used when `type` in `flag` is 1. It specifies the boundaries of the plot as the vector `[xmin,xmax,ymin,ymax,zmin,zmax]`.
- `zlev` : real number.

### Description

`contour` draws level curves of a surface z=f(x,y). The level curves are drawn on a 3D surface. The optional arguments are the same as for the function `plot3d` (except `zlev`) and their meanings are the same. They control the drawing of level curves on a 3D plot. Only `flag(1)=mode` has a special meaning.

- mode=0 : the level curves are drawn on the surface defined by (x,y,z).
- mode=1 : the level curves are drawn on a 3D plot and on the plan defined by the equation z=zlev.
- mode=2 : the level curves are drawn on a 2D plot.

You can change the format of the floating point number printed on the levels by using `xset("fpf",string)` where `string` gives the format in C format syntax (for example `string="%.3f"`). Use `string=""` to switch back to default format and Use `string=" "` to suppress printing.

Usually we use `contour2d` to draw levels curves on a 2D plot.

Enter the command `contour()` to see a demo.

### Examples

    t=%pi*[-10:10]/10;
    deff("[z]=surf(x,y)","z=sin(x)*cos(y)"); z=feval(t,t,surf);
    rect=[-%pi,%pi,-%pi,%pi,-1,1];
    contour(t,t,z,10,35,45," ",[0,1,0],rect)
    // changing the format of the printing of the levels
    xset("fpf","%.2f")
    xbasc()
    contour(t,t,z,10,flag=[0,1,0],ebox=rect)

### See Also

`contour2d` `fcontour` `fcontour2d` `plot3d` `xset`

### Author

J.Ph.C.