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

int main() {
    int numpoints = 36;
    array double theta[36], z[20], r[720];
    int i, j;
    class CPlot plot;

    linspace(theta, 0, 360);
    linspace(z, 0, 2*M_PI);
    for(i=0; i<36; i++) {
        for(j=0; j<20; j++) {
            r[i*20+j] = 2+cos(z[j]);
        }
    }
    plot.dimension(3);
    plot.coordSystem(PLOT_COORD_CYLINDRICAL, PLOT_ANGLE_DEG);
    plot.data3D(theta, z, r);
    plot.axisRange(PLOT_AXIS_XY, -4, 4);
    plot.plotting();
}