Mivel a régi Dopplerrel számolva pontosabbak a kapott eredmények, szinte biztos, hogy az egész tényleg csak játék a számokkal, hiszen a relativitás a pontosabb modellje a valóságnak.
r:7.044850718840330e-16 m E:5.581579759870600e+02 MeV 1.999560278517446e+00
l:6.065775539802652e-13 m E:5.109989099021028e-01 MeV 1.000000000106580e+00 x m(e)
r:5.291772082388662e-11 m E:27.21138 eV
l(qm):3.324829930671497e-10 m
r(qm):5.291631184062524e-11 m
double
h=6.62606896e-34,
me=9.10938215e-31,
c=2.99792458e8,
k=8.9875517874e9,
e=1.602176487e-19,
alpha=1.0/7.2973525376e-3,
mpion=139.57018,
E,l,m,v,b,f,f1,f2,f0,r,F;
int i;
v=c;
r=k*e*e/(me*v*v*4.0);// F=m*v*v/r; F=k*e*e/(4.0*r*r);
l=r*2.0*M_PI;
f0=v/l;
m=h*f0/(v*v);//E=hf E=mv2
E=(m*c*c - me*c*c*2.0)/(e*1e6);
printf("r:%.15e m ",r);
printf("E:%.15e MeV %.15e n",E*2.0,E/mpion );
v=c/alpha;
//f1=f0*sqrt((1.0-v/c)/(1.0+v/c));
//f2=f0*sqrt((1.0+v/c)/(1.0-v/c));
f1=f0*(1.0-v/c);
f2=f0*(1.0+v/c);
l=2.0*c/(f2-f1);
printf("l:%.15e m ",l);
f=c/l;
f/=4.0;
m=h*f/(c*c);
E=(m*c*c)/(e*1e6);
printf("E:%.15e MeV %.15e x m(e) n",E,m/me);
//f1=f*sqrt((1.0-v/c)/(1.0+v/c));
//f2=f*sqrt((1.0+v/c)/(1.0-v/c));
f1=f*(1.0-v/c);
f2=f*(1.0+v/c);
l=2.0*c/(f2-f1);
r=l/(M_PI*2.0);
printf("r:%.15e m ",r);
f=v/l;
m=h*f/(c*c);
E=(m*c*c)/e;
printf("E:%.5f eV nn",E);
b=1.0/sqrt(1.0-v*v/(c*c));
l=h/(v*me*b);
printf("l(qm):%.15e mn",l);
r=l/(M_PI*2.0);
printf("r(qm):%.15e mn",r);