... Wade Ward <zaxfuuq@gmail.com> wrote: I think this version gives as much width on epsilon nought as double precision can do. I can't believe I got stung with a single-...+00 2. The only use of ** in your program was c**2.0_dp so why not write c**2 instead? -- John Harper, School of Mathematics, Statistics and Computer Science, Victoria University, PO ...

...of atan. So I'd leave c as is, but change the statements...4e-7_dp * PIE E0 = 1/ (U0*c**2) ! because **2 is better for...real(kind=dp):: U0,E0,PIE, c, e0_calc PIE = 4.0_dp* atan(1...* PIE E0= 8.854187E-12_dp c = 2.99792458E8_dp e0_calc = 1.0_dp/ ...c**2.0_dp) write (*,'(f35.25)') c, U0, PIE,E0, e0_calc write (*,*)... gives as much width on epsilon nought as double precision can do. I can't believe I ...

... Resnick.) It doesn't help that it's pronuounced mu nought and looks like u zero. I think Jan ... (the A being Ampere) is magnetic based, such that mu nought is based on geometry and to make the units ... on the force between wires with a specific geometry which sets the value of mu nought in the appropriate equations. Epsilon nought is then determined through Maxwell's equations to satisfy ...

... MM posting a value for c at tens of thousands of ...x 10**8, this forces either C =(E0*U0)**(-.5) to be an... but in SI units both C and U0 above are both ... at the above link that c is defined as a rational ... problem? They're still constants: c = 299792458 m.s^-1 (exact)...^-2 (exact) e0 = 1/(u0.c^2) F.m^-1 (exact)... The irrationality of mu nought and epsilon nought isn't a ...

... pi's cancel leaving the inverse square root of (10**-16 * ...! another statement function asf(arg1) = SQRT(arg1-OI)**2 E0=(1E-9...always lobbied for the existence of C-style comments in fortran, but found...x 10**8, this forces either C =(E0*U0)**(-.5) to be an..., we never calculated epsilon nought or mu nought. These were constants in... integration. The formula for the inverse of the square root of ...

...I have another problem with Wade's terminology. c is a velocity, which happens to be a...metres per second, which Paul gave. c is NOT a rational number. Rational ...) real(kind=dp):: x,U0,E0,PIE, c, e0_calc PIE = 4* atan(1.0) U0=(...4.0 * PIE E0= 8.854187E-12_dp c = 2.99792458E8_dp x = (E0*U0)**(-.5) write..."going the other way" e0_calc = 1/ (U0 * c**2.0) write (*,'(f35.25)') c, ...

John Harper wrote: (snip) I have another problem with Wade's terminology. c is a velocity, which happens to be a rational (for that matter an integer) number of metres per second, which Paul gave. c is NOT a rational number. Rational numbers are dimensionless. Remember the spacecraft that crashed on Mars because someone was ...

...Wade Ward <zaxf...@gmail.com> wrote: C =(E0*U0)**(-.5) to be an..., but in SI units both C and U0 above are both exact...see at the above link that c is defined as a rational slightly...a problem? They're still constants: c = 299792458 m.s^-1 (exact) ... problem with Wade's terminology. c is a velocity, which happens ... per second, which Paul gave. c is NOT a rational number. ...

...With MM posting a value for c at tens of thousands of meters...x 10**8, this forces either C =(E0*U0)**(-.5) to be an..., but in SI units both C and U0 above are both exact...kg m s**(-2) A**(-2), C as m/s, E0 as F...see at the above link that c is defined as a rational slightly...a problem? They're still constants: c = 299792458 m.s^-1 (exact)...^-2 (exact) e0 = 1/(u0.c^2) F.m^-1 (exact)...

...> With MM posting a value for c at tens of thousands of meters per second less than 3 x 10**8, this forces either C =(E0*U0)**(-.5) to be an approximation or ... you are using, but in SI units both C and U0 above are both exact, E0 is an.../A**2 = kg m s**(-2) A**(-2), C as m/s, E0 as F/m = kg... the net and see at the above link that c is defined as a rational slightly less than 3...