3_PowerRadiated/measuredPowerLienardWiechert/measuredPowerLienardWiechert.c
//magnitude of measured poynting vector for non-relativistic speeds:
magnitude of measured poynting vector for non-relativistic speeds.
magnitude of measured poynting vector for non - relativistic speeds .
double nonRelMeasuredPoyntingVecMag(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 radius = (eMass*linearV)/(e*B); 
  double centripetalAccel = pow(linearV, 2.0)/radius;
  double factor1 = 1.0/(mu0*c);
  double factor2 = (e/(4.0*M_PI*epsilon0*pow(c,2.0)*R))*centripetalAccel;
  double poyntingMag = factor1*pow(factor2,2.0)*pow(sin(angle),2.0)*effArea;
  return poyntingMag;
}
double nonRelMeasuredPoyntingVecMag ( 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 radius = ( eMass * linearV ) / ( e * B ) ; double centripetalAccel = pow ( linearV , 2.0 ) / radius ; double factor1 = 1.0 / ( mu0 * c ) ; double factor2 = ( e / ( 4.0 * M_PI * epsilon0 * pow ( c , 2.0 ) * R ) ) * centripetalAccel ; double poyntingMag = factor1 * pow ( factor2 , 2.0 ) * pow ( sin ( angle ) , 2.0 ) * effArea ; return poyntingMag ; }
magnitude of measured poynting vector for non-relativistic speeds:
magnitude of measured poynting vector for non - relativistic speeds :
//Lamor radius/gyroradius/cyclotron radius