#include<geometry.h>

int main()
{
cout << "number of dimensions: ";
int nd; cin >> nd; cout << nd << endl;
if (nd == 2) {
   cout << "point; ends of line: ";
   Vector2d p, a, b; cin >> p >> a >> b;
   cout << p << '\t' << a << '\t' << b << endl;
   double angsub = angsub2d(p,a,b);
   cout << "angle subtended: " << angsub << " radians\n";
} else if (nd == 3) {
   cout << "point; vertices of triangle: ";
   Vector3d p, a, b, c; cin >> p >> a >> b >> c;
   cout << p << '\t' << a << '\t' << b << '\t' << c << endl;
   double angsub = angsub3d(p,a,b,c);
   cout << "angle subtended: " << angsub << " radians\n";
} else { cerr << "Error: 2D/3D only!\n"; return(1); }
return 0;
}