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

int main() {
    int i, j;
    array double theta[37], phi[19], r[703];
    class CPlot plot;

    linspace(theta, 0, 2*M_PI);
    linspace(phi, -M_PI/2, M_PI/2);
    r = (array double [703])1;
    plot.dimension(3);
    plot.data3D(theta, phi, r);
    plot.coordSystem(PLOT_COORD_SPHERICAL);
    plot.axisRange(PLOT_AXIS_XY, -2.5, 2.5);
    plot.plotting();
}