A hagyományos Dopplert elrontottam, helyesen:

r:7.044850718840330e-16 m E:1.397949934516888e+02 MeV 1.001610755619064e+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.324918460055154e-10 m
r(qm):5.291772082952704e-11 m




#include <stdio.h>
#include <stdlib.h>
#include <math.h>


int main(){

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*1.0)/(e*1e6);E/=2.0;
printf("r:%.15e m ",r);
printf("E:%.15e MeV %.15e n",E,E/mpion );


v=c/alpha;


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*(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);





l=h/(v*me);
printf("l(qm):%.15e mn",l);
r=l/(M_PI*2.0);
printf("r(qm):%.15e mn",r);



return 0;
}