diff --git a/src/jat/examples/ThreeBodyExample/ThreeBodyAPL.java b/src/jat/examples/ThreeBodyExample/ThreeBodyAPL.java deleted file mode 100644 index 41a565c..0000000 --- a/src/jat/examples/ThreeBodyExample/ThreeBodyAPL.java +++ /dev/null @@ -1,47 +0,0 @@ -package jat.examples.ThreeBodyExample; - -import org.apache.commons.math3.ode.FirstOrderDifferentialEquations; - -import jat.coreNOSA.cm.ThreeBody; -import jat.coreNOSA.math.MatrixVector.data.VectorN; - -public class ThreeBodyAPL extends ThreeBody implements FirstOrderDifferentialEquations{ - - public ThreeBodyAPL(double G, double m1, double m2, double m3) { - super(G, m1, m2, m3); - } - - @Override - public void computeDerivatives(double arg0, double[] arg1, double[] arg2) { - // TODO Auto-generated method stub - - } - - @Override - public int getDimension() { - // returns a dimension of 6 - return 6; - } - public double[] randv() { - double[] randv = new double[18]; - randv[0] = 2.0; - randv[1] = 6.0; - randv[2] = 7.0; - randv[3] = 4.0; - randv[4] = 8.0; - randv[5] = 7.0; - randv[6] = 5.0; - randv[7] = 8.0; - randv[8] = 3.0; - randv[9] = 6.0; - randv[10] = 8.0; - randv[11] = 9.0; - randv[12] = 0.0; - randv[13] = 7.0; - randv[14] = 4.0; - randv[15] = 5.0; - randv[16] = 7.0; - randv[17] = 8.0; - return randv; - } -} diff --git a/src/jat/examples/ThreeBodyExample/ThreeBodyExample.java b/src/jat/examples/ThreeBodyExample/ThreeBodyExample.java deleted file mode 100644 index a77a8db..0000000 --- a/src/jat/examples/ThreeBodyExample/ThreeBodyExample.java +++ /dev/null @@ -1,94 +0,0 @@ -package jat.examples.ThreeBodyExample; - -import jat.core.cm.TwoBodyAPL; -import jat.core.plot.plot.FrameView; -import jat.core.plot.plot.Plot2DPanel; -import jat.core.plot.plot.PlotPanel; -import jat.core.plot.plot.plots.ScatterPlot; - -import java.awt.Color; -import java.text.DecimalFormat; - -import javax.swing.JFrame; - -import org.apache.commons.lang3.ArrayUtils; -import org.apache.commons.math3.linear.ArrayRealVector; -import org.apache.commons.math3.linear.RealVectorFormat; -import org.apache.commons.math3.ode.FirstOrderIntegrator; -import org.apache.commons.math3.ode.nonstiff.DormandPrince853Integrator; - -public class ThreeBodyExample { - - public ThreeBodyExample() { - } - - public static void main(String[] args) { - - ThreeBodyExample x = new ThreeBodyExample(); - - // create a TwoBody orbit using orbit elements - TwoBodyAPL sat = new TwoBodyAPL(7000.0, 0.3, 0.0, 0.0, 0.0, 0.0); - - double[] y = sat.randv(); - - ArrayRealVector v = new ArrayRealVector(y); - - DecimalFormat df2 = new DecimalFormat("#,###,###,##0.00"); - RealVectorFormat format = new RealVectorFormat(df2); - System.out.println(format.format(v)); - - // find out the period of the orbit - double period = sat.period(); - - // set the final time = one orbit period - double tf = period; - - // set the initial time to zero - double t0 = 0.0; - - // propagate the orbit - FirstOrderIntegrator dp853 = new DormandPrince853Integrator(1.0e-8, 100.0, 1.0e-10, 1.0e-10); - dp853.addStepHandler(sat.stepHandler); - // double[] y = new double[] { 7000.0, 0, 0, .0, 8, 0 }; // initial - // state - - dp853.integrate(sat, 0.0, y, 8000, y); // now y contains final state at - // tf - - Double[] objArray = sat.time.toArray(new Double[sat.time.size()]); - double[] timeArray = ArrayUtils.toPrimitive(objArray); - double[] xsolArray = ArrayUtils.toPrimitive(sat.xsol.toArray(new Double[sat.time.size()])); - double[] ysolArray = ArrayUtils.toPrimitive(sat.ysol.toArray(new Double[sat.time.size()])); - - double[][] XY = new double[timeArray.length][2]; - - // int a=0; - // System.arraycopy(timeArray,0,XY[a],0,timeArray.length); - // System.arraycopy(ysolArray,0,XY[1],0,ysolArray.length); - - for (int i = 0; i < timeArray.length; i++) { - XY[i][0] = xsolArray[i]; - XY[i][1] = ysolArray[i]; - } - - Plot2DPanel p = new Plot2DPanel(); - - // Plot2DPanel p = new Plot2DPanel(min, max, axesScales, axesLabels); - - ScatterPlot s = new ScatterPlot("orbit", Color.RED, XY); - // LinePlot l = new LinePlot("sin", Color.RED, XY); - // l.closed_curve = false; - // l.draw_dot = true; - p.addPlot(s); - p.setLegendOrientation(PlotPanel.SOUTH); - double plotSize = 10000.; - double[] min = { -plotSize, -plotSize }; - double[] max = { plotSize, plotSize }; - p.setFixedBounds(min, max); - - new FrameView(p).setDefaultCloseOperation(JFrame.EXIT_ON_CLOSE); - - System.out.println("end"); - - } -}