code: def fa_trace_to_lambdas(fa=0.08, trace=0.0021):
      lambda1 = (trace / 3.0) * (
          1 + 2 * fa / (3 - 2 * fa**2) ** (1 / 2.0)
      )
      lambda2 = (trace / 3.0) * (
          1 - fa / (3 - 2 * fa**2) ** (1 / 2.0)
      )
      evals = np.array([lambda1, lambda2, lambda2])
      return evals
  
imports: import warnings
  import numpy
  import numbers
  from scipy.integrate import quad
  from scipy.special import lpn, gamma
  import scipy.linalg
  import scipy.linalg.lapack
  from dipy.data import small_sphere, get_sphere, default_sphere
  from dipy.core.geometry import cart2sphere
  from dipy.core.ndindex import ndindex
  from dipy.sims.voxel import single_tensor
  from dipy.reconst.multi_voxel import multi_voxel_fit
  from dipy.reconst.dti import TensorModel, fractional_anisotropy
  from dipy.reconst.shm import sph_harm_ind_list, real_sh_descoteaux_from_index, sph_harm_lookup, lazy_index, SphHarmFit, real_sh_descoteaux, sh_to_rh, forward_sdeconv_mat, SphHarmModel
  from dipy.reconst.utils import _roi_in_volume, _mask_from_roi
  from dipy.direction.peaks import peaks_from_model
  from dipy.core.geometry import vec2vec_rotmat
  from dipy.utils.deprecator import deprecate_with_version, deprecated_params
label_desc: Use image node_id 14 for calling a global function with example usage: fa_trace_to_lambdas(fa, trace) and returns: evals