Regularized solutions

# Imports
from gfeatpy.observation import Range
from gfeatpy import plotting
from gfeatpy.gravity import GravityAnomaly
import matplotlib.pyplot as plt
import numpy as np

# Define error spectra
mwi = lambda f: 2.62 * np.sqrt(1 + (0.003/f)**2) * 1e-6
acc = lambda f: 1e-10 * np.sqrt(1 + (f/0.5)**4 + (0.005/f))
# Define observation settings
l_max = 60
Nr = 16
Nd = 1
I = np.deg2rad(60)
rho=220e3
# Set up system and solve
range = Range(l_max, Nr, Nd, I, rho)
range.set_observation_error(mwi, acc)
range.set_kaula_regularization()
range.solve()

# Plot results
plotting.synthesis(range, 360, 180, GravityAnomaly())
plotting.ground_track(Nr, Nd, I, samples_per_rev =5_000, s=0.01, color='k')
plt.show()

# Plot results (no ground-track)
plotting.synthesis(range, 360, 180, GravityAnomaly())
plt.show()