#include <chplot.h>
#include <math.h>

int main() {
    array double x[20], y[20], z[2][400];
    int i, j;
    class CPlot plot;

    linspace(x, -2, 2);
    linspace(y, -2, 2);
    for (i=0; i<20; i++) {
        for(j=0; j<20; j++) {
            z[0][i*20+j] = x[i]*exp(-x[i]*x[i]-y[j]*y[j]);
            z[1][i*20+j] = z[0][i*20+j] +2;
        }
    }
    plot.dimension(3);
    plot.data3D(x, y, z);
    plot.plotting();
}