1; % some physical constants hbar=1; im_one=0+1i; el_charge=1;