def METHOD_NAME(xrs, map_data, n_steps=300, step_size=0.01):
  atom_center = xrs.scatterers()[0].site
  a_cell = xrs.unit_cell().parameters()[0]
  dist = flex.double()
  rho = flex.double()
  x=-10
  while x<=10:
    x += step_size
    point = x/a_cell, 0, 0
    density_value_at_point_e = map_data.eight_point_interpolation(point)
    dist.append(x)
    rho.append(density_value_at_point_e)
  return dist, rho