double magBol = 12.1; // (Compute from bolometric corrections in Allen)
int pop = 6; // (Besancon model. 6: 5-7 Gyr)
String spType = "M4 V"; // String defining the spectral type
// Generate simulation
Simulation sim = new Simulation();
// Generate solar system, set parameters, generate planet and star; finalize
SimulatedSystem system = sim.addSystem("GJ 876");
system.setAstrometry(ra, dec, parallax, muRa, muDec, vRad);
system.createStar(magMv, distance, vMinusI, absV, feH, alphaE, mass, spType, magBol, pop);
system.addPlanet(massPlanet, period, timePeriastron, eccentricity, inclination, omega2, nodeAngle);
system.finalizeSystem();//*/