import numpy as np

gamma = 1/26.24e-9
h = 6.626e-34
c = 3e8
lam = 780.24e-9

p0= h*c*gamma/lam/8
p_fit = 25.6e-12
p_fit_err = 0.6e-12

spatial_lambda = p0/p_fit
spatial_lambda_err = p_fit_err*p0/p_fit**2	

print('lambda:'+ str(spatial_lambda) + ' +- '+ str(spatial_lambda_err))