#ifndef POINTH
#define POINTH
#include <iostream>
using namespace std;

class point {
	double x;   //Cartesian coordinates
	double y;
	double z;
public:
	point(double initial_x = 0, double initial_y = 0, double initial_z = 0)
		: x(initial_x), y(initial_y), z(initial_z) {}

	point& xrot(double theta);   //theta in radians
	point& yrot(double theta);
	point& zrot(double theta);

	friend ostream& operator<<(ostream& ost, const point& p) {
		return ost << "(" << p.x << ", " << p.y << ", " << p.z << ")";
	}
};
#endif