3_PowerRadiated/measuredPowerLienardWiechert/measuredPowerLienardWiechert.c
//magnitude of measured poynting vector for relativistic speeds:
magnitude of measured poynting vector for relativistic speeds.
magnitude of measured poynting vector for relativistic speeds .
double relMeasuredPoyntingVecMag(double R, double angle, double linearV, double effArea, double B) {
  double mu0 = 1.25663706*pow(10.0,-6.0);
  double epsilon0 = 8.8541878128*pow(10.0,-12.0);
  double c = 299792458.0;
  double e = 1.602176634*pow(10.0,-19.0);
  double eMass = 9.10938356*pow(10.0,-31.0);
  double beta = linearV/c;
  double lorentzFactor = 1.0/pow((1.0 - pow(beta,2.0)),0.5);
  double radius = (lorentzFactor*eMass*linearV)/(e*B); 
  double centripetalAccel = (e*linearV*B)/(lorentzFactor*eMass);
  double factor1 = 1.0/(mu0*c);
  double factor2 = (e/(4.0*M_PI*epsilon0*pow(c,2.0)*R))*centripetalAccel;
  double factor3 = 1.0/pow((1.0 - beta*sin(angle)),3.0);
  double factor4 = 1.0 - pow(cos(angle),2.0)/(pow(lorentzFactor,2.0)*pow((1.0 - beta*sin(angle)),2.0));
  double poyntingMag = factor1*pow(factor2,2.0)*factor3*factor4*effArea;
  return poyntingMag;
}
double relMeasuredPoyntingVecMag ( double R , double angle , double linearV , double effArea , double B ) { double mu0 = 1.25663706 * pow ( 10.0 , -6.0 ) ; double epsilon0 = 8.8541878128 * pow ( 10.0 , -12.0 ) ; double c = 299792458.0 ; double e = 1.602176634 * pow ( 10.0 , -19.0 ) ; double eMass = 9.10938356 * pow ( 10.0 , -31.0 ) ; double beta = linearV / c ; double lorentzFactor = 1.0 / pow ( ( 1.0 - pow ( beta , 2.0 ) ) , 0.5 ) ; double radius = ( lorentzFactor * eMass * linearV ) / ( e * B ) ; double centripetalAccel = ( e * linearV * B ) / ( lorentzFactor * eMass ) ; double factor1 = 1.0 / ( mu0 * c ) ; double factor2 = ( e / ( 4.0 * M_PI * epsilon0 * pow ( c , 2.0 ) * R ) ) * centripetalAccel ; double factor3 = 1.0 / pow ( ( 1.0 - beta * sin ( angle ) ) , 3.0 ) ; double factor4 = 1.0 - pow ( cos ( angle ) , 2.0 ) / ( pow ( lorentzFactor , 2.0 ) * pow ( ( 1.0 - beta * sin ( angle ) ) , 2.0 ) ) ; double poyntingMag = factor1 * pow ( factor2 , 2.0 ) * factor3 * factor4 * effArea ; return poyntingMag ; }
magnitude of measured poynting vector for relativistic speeds:
magnitude of measured poynting vector for relativistic speeds :
//Lamor radius/gyroradius/cyclotron radius for relativistic motion //double centripetalAccel = (pow(lorentzFactor,2.0)*pow(linearV, 2.0))/radius;