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

int main() {
    array double theta1[37], phi1[37], r1[37];
    array double theta2[37], phi2[37], r2[37];
    array double theta3[20], phi3[20], r3[20];
    class CPlot plot;

    linspace(theta1, 0, 2*M_PI);
    phi1 = (array double [37])0;
    r1 = (array double [37])1;
    theta2 = (array double [37])M_PI/2;
    linspace(phi2, 0, 2*M_PI);
    r2 = (array double [37])1;
    theta3 =  (array double [20])0;
    linspace(phi3, -M_PI/2, M_PI/2);
    r3 = (array double [20])4;
    plot.dimension(3);
    plot.data3D(theta1, phi1, r1);
    plot.data3D(theta2, phi2, r2);
    plot.data3D(theta3, phi3, r3);
    plot.point(0, 0, 0);
    plot.coordSystem(PLOT_COORD_SPHERICAL);
    plot.axisRange(PLOT_AXIS_XY, -2, 7);
    plot.size3D(1, 2);
    chplot.hiddenLineRemoval(PLOT_OFF);
    plot.plotting();
}