contour2di 2次元プロット上の曲面の等高線を計算

## 呼び出し手順

    [xc, yc] = contour2di(x, y, z, nz)

## 引数

x,y

  大きさ`n1`および`n2`の 実数行ベクトル: グリッド.

z

  大きさ `(n1,n2)`の実数行列, 関数の値.

nz

  等高線の値または等高線の数.

  `nz`が整数の場合,

    この値は等高線の数を 指定し, `zmin`から`zmax`の間で 以下のように等間隔に配置されます:

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

    `zmin` および `zmax`の 等高線は描画されません ( 通常,これらは点に縮小されます)が,以下のコマンドで追加できることに 注意してください.

        [im,jm] = find(z == zmin);     // または zmax

        plot2d(x(im)',y(jm)',-9,"000")

  `nz` がベクトルの場合,

    `nz(i)` は, `i`番目の等高線の値を指定します.

xc,yc

  等高線の定義を含む同じ大きさのベクトル. 詳細は以下を参照.

## 説明

`contour2di` は,2次元プロット上の曲面 `z=f(x,y)` の等高線を計算します. `f(x,y)` の値は,行列`z`により `x` および `y`で定義される グリット点において指定します.

`xc(1)` には最初の等高線のパスに関する値が含まれ, `yc(1)`は, この等高線のパスを定義する`N1`個の点を有します. (`xc(1+(1:N1))`, `yc(1+(1:N1))` )nは パスの点の座標を有します. 2番目のパスは,`xc(2+N1)` および `yc(2+N1)` で始まるといったようになります.

## 例

    [xc,yc]=contour2di(1:10,1:10,rand(10,10),5);

    k=1;n=yc(k);c=1;

    while k+yc(k)<size(xc,'*')

      n=yc(k);

      plot2d(xc(k+(1:n)),yc(k+(1:n)),c)

      c=c+1;

      k=k+n+1;

    end

## 参照

- contour
- contour2d
- plot2d