what are some problems and bugs with the following code? ""import numpy as np import matplotlib.pyplot as plt from astropy import constants as const from astropy import units as u from astroquery.simbad import Simbad from astroML.plotting import scatter_contour

Define classes for Vec2 and Particle similar to the previous program Define class for binary system (combining Hellion and Sarnan) class BinarySystem: def init(self, hellion_data, sarnan_data, orbital_data): # Initialize the binary system using the provided data # Set up the stars (Hellion and Sarnan) with their properties # Set up the orbital parameters (semi-major axis, eccentricity, etc.) pass

Check the units of the input data for key, value in orbital_data.items(): if not isinstance(value, u.Quantity): raise ValueError(f'The value {key} must be a Quantity.')

... (rest of the class methods)
Main code to create the binary system and perform simulations def main(): # Data for Hellion hellion_data = { 'Age': 5.2e9 * u.year, 'Radius': 699226.7 * u.km, 'Hydrogen Abundance': 71.2 * 0.01, # Convert from percentage to fraction 'Helium Abundance': 8.3 * 0.01, # Convert from percentage to fraction 'Metallicity': 0.38 * 1e-2, # Convert from percentage to fraction 'Mass': 2.148228e30 * u.kg, 'Density': 1.6736e9 * u.kg / u.km3, 'Volume': 1.31316e30 * u.km3, 'Average Density': 1.634e9 * u.kg / u.km3, 'Total Binding Energy (U1)': 2.3096e36 * u.joule, 'Central Density (Hellion)': 1500.1605932993 * u.kg / u.m3, 'Density at r = 500,000 km': 4.633 * u.kg / u.km**3, 'Composition': { 'Hydrogen (H)': 71.2 * 0.01, 'Helium (He)': 8.3 * 0.01, 'Carbon (C)': 5.2 * 0.01, 'Nitrogen (N)': 2.4 * 0.01, 'Oxygen (O)': 0.77 * 0.01, 'Silicon (Si)': 6 * 0.01, 'Iron (Fe)': 0.16 * 0.01, 'Sulfur (S)': 0.04 * 0.01, 'NAME_1 (B)': 2.9 * 0.01, 'Sodium (Na)': 1.5 * 0.01, 'Magnesium (Mg)': 4.8 * 0.01 * 0.98, # Reduced by 2% from its previous value } }

Data for Sarnan
sarnan_data = { 'Age': 5.2e9 * u.year, 'Radius': 699496.94 * u.km, 'Hydrogen Abundance': 73.1 * 0.01, # Convert from percentage to fraction 'Helium Abundance': 8.0 * 0.01, # Convert from percentage to fraction 'Metallicity': 0.4 * 1e-2, # Convert from percentage to fraction 'Mass': 2.14751e30 * u.kg, 'Density': 1.6730e9 * u.kg / u.km3, 'Volume': 1.31377e30 * u.km3, 'Average Density': 1.634e9 * u.kg / u.km3, 'Total Binding Energy (U2)': 2.3101e36 * u.joule, 'Central Density (Sarnan)': 1497.921756058127 * u.kg / u.m3, 'Composition': { 'Hydrogen (H)': 73.1 * 0.01, 'Helium (He)': 8.0 * 0.01, 'Carbon (C)': 5.2 * 0.01, 'Nitrogen (N)': 2.4 *