#ifndef vector_h
#define vector_h

#include <iostream>
#include <math.h>
using namespace std;
//矢量球坐標系數據結構
typedef struct
{
	double R;
	double thinc;
	double phinc;
}doubleVector;

//矢量直角坐標系數據結構
typedef struct
{
	double x;
	double y;
	double z;
}doublexyzVector;

//////////////////////////////////////////////////////////////////////
//                                  
//                           坐標轉化和矢量運算
//
//////////////////////////////////////////////////////////////////////
doubleVector Vectorset(double r,double thinc,double phinc);
doublexyzVector RTPtoXYZ(doubleVector r);						//將球坐標轉化為直角坐標
doubleVector XYZtoRTP(doublexyzVector x);						//將直角坐標轉化為球坐標
doublexyzVector XYZadd(doublexyzVector x,doublexyzVector y);	//兩矢量相加
doublexyzVector RTPsub(doubleVector rise,doubleVector end);	    //兩矢量相減
doublexyzVector XYZsub(doublexyzVector x,doublexyzVector y);	//兩矢量相減
doublexyzVector XYZmult(double c,doublexyzVector);				//矢量乘常數
doublexyzVector XYZVMult(doublexyzVector X,doublexyzVector Y);  //矢量杈乘
double dotmult(doubleVector X,doubleVector Y);					//矢量點乘
double XYZdotmult(doublexyzVector x,doublexyzVector y);			//直角坐標系矢量點乘
doublexyzVector Enxyz(doublexyzVector X);						//求直角坐標系中某矢量的單位矢量
doubleVector En(doubleVector X);								//求球坐標系中某矢量的單位矢量
double MO(doublexyzVector x);									//矢量的模
double area(doubleVector A,doubleVector B,doubleVector C);      //三角形的面積
void printdoubleVector(doubleVector r);                         //顯示點的球坐標系的參量
void printdoublexyzVector(doublexyzVector r);                   //顯示點的直角坐標系的參量

#endif