import numpy as np import math from datetime import datetime, timedelta from pymoo.core.problem import elementwiseproblem k_legal = 0.0005624701263657791 b=6371 b=7371 phi=30.35 pi=3.1415926 standard_meridian=120 mean1=9.38167 dev1=0.22401 mean2=4816.70238 dev2=164.0126 def calculate_length(phi, delta, omega, b, b): phi_rad = np.radians(phi) delta_rad = np.radians(delta) omega_rad = np.radians(omega) r_prime = b * np.sin(phi_rad - delta_rad) r_double_prime = b * np.cos(phi_