package astro.galaxy;

import java.awt.*;

import myutil.*;

public class Data 
{	
	static public double m[];
	static public double x[];
	static public double y[];
	static public double vx[];
	static public double vy[];
	static public double rad[];
	static public Color  col[];
	static public double[][] as;
	static public double[][] ac;
	static public double[][] bs;
	static public double[][] bc;
	static public int M=0, N=0, K=0;

        static public void init(String name, int n, double G)
        {
	    int i,j,k;
	    double r=0, th=0, v=0;
	    double v_escape=0;
	    double r0=0.01;
	    int n7;

	    n7 = (n<7)?7:n;

	    m = new double[n7];
	    x = new double[n7];
	    y = new double[n7];
	    vx = new double[n7];
	    vy = new double[n7];
	    rad = new double[n7];
	    col = new Color[n7];
	    as = null;
	    ac = null;
	    bs = null;
	    bc = null;

		if (name.equals("LagrangeN")) {
		  for (j=0; j<n; j++) {
		    r = 0.1;
		    th = (24.0*j)/n;
		    v_escape = Math.sqrt(3*G*(n-1)/r0); 
		    v = 0.05*v_escape;
		    x[j] = r*Math.cos(th*2*Math.PI/24);
		    y[j] = r*Math.sin(th*2*Math.PI/24);
		    vx[j] = -v*Math.sin(th*2*Math.PI/24);
		    vy[j] = v*Math.cos(th*2*Math.PI/24);
		  }
		} else
		if (name.equals("LagrangeNa")) {
		  for (j=0; j<n; j++) {
		    r = 0.0;
		    th = (24.0*j)/n;
		    v_escape = Math.sqrt(3*G*(n-1)/r0); 
		    v = 0.55*v_escape;
		    x[j] = r*Math.cos(th*2*Math.PI/24);
		    y[j] = r*Math.sin(th*2*Math.PI/24);
		    vx[j] = -v*Math.sin(th*2*Math.PI/24);
		    vy[j] = v*Math.cos(th*2*Math.PI/24);
		  }
		} 

	        for (j=0; j<n; j++) {
		    // m[j] = 300000;
		    m[j] = 1;
		    rad[j] = r0;
		    col[j] = new Color((float)(0.5+Math.random()/2), 
				       (float)(0.5+Math.random()/2), 
			               (float)(0.5+Math.random()/2)) ;
		}


		if (name.equals("Rosette4")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    x[0]=   -1.19261140209982; y[0]=   -0.67343666483354; 
		    x[1]=   -0.13745545400244; y[1]=   -0.31273531249149; 
		    x[2]=    0.41488315615986; y[2]=    1.06708402421544; 
		    x[3]=    0.91518369994239; y[3]=   -0.08091204689040; 
		    vx[0]=   -0.15823018219881; vy[0]=    0.03389131793419; 
		    vx[1]=   -0.76384745689789; vy[1]=    1.32594980724075; 
		    vx[2]=    0.03336874686611; vy[2]=   -0.79939485612001; 
		    vx[3]=    0.88870889223059; vy[3]=   -0.56044626905494; 
		    M=80; N=4; K=12;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
		    as[0][1]=  -0.28187; ac[0][1]=  -1.21537; 
		    as[0][3]=   0.03455; ac[0][3]=   0.02686; 
		    as[0][5]=  -0.00066; ac[0][5]=  -0.00379; 
		    as[0][7]=  -0.00277; ac[0][7]=  -0.00016; 
		    as[0][9]=   0.00128; ac[0][9]=  -0.00025; 
		    as[1][1]=  -0.70679; ac[1][1]=  -0.11466; 
		    as[1][3]=  -0.03455; ac[1][3]=  -0.02686; 
		    as[1][5]=   0.00066; ac[1][5]=   0.00379; 
		    as[1][7]=   0.00295; ac[1][7]=   0.00053; 
		    as[1][9]=   0.00136; ac[1][9]=  -0.00015; 
		    as[2][1]=  -0.05603; ac[2][1]=   0.45255; 
		    as[2][3]=   0.02686; ac[2][3]=  -0.03455; 
		    as[2][5]=   0.00379; ac[2][5]=  -0.00066; 
		    as[2][7]=   0.00025; ac[2][7]=  -0.00304; 
		    as[2][9]=  -0.00137; ac[2][9]=   0.00024; 
		    as[3][1]=   1.04468; ac[3][1]=   0.87747; 
		    as[3][3]=  -0.02686; ac[3][3]=   0.03455; 
		    as[3][5]=  -0.00379; ac[3][5]=   0.00066; 
		    as[3][7]=  -0.00043; ac[3][7]=   0.00267; 
		    as[3][9]=  -0.00127; ac[3][9]=   0.00016; 
		    bs[0][1]=   0.11466; bc[0][1]=  -0.70679; 
		    bs[0][3]=  -0.02686; bc[0][3]=   0.03455; 
		    bs[0][5]=  -0.00379; bc[0][5]=   0.00066; 
		    bs[0][7]=   0.00053; bc[0][7]=  -0.00295; 
		    bs[0][9]=   0.00015; bc[0][9]=   0.00136; 
		    bs[1][1]=   1.21537; bc[1][1]=  -0.28187; 
		    bs[1][3]=   0.02686; bc[1][3]=  -0.03455; 
		    bs[1][5]=   0.00379; bc[1][5]=  -0.00066; 
		    bs[1][7]=  -0.00016; bc[1][7]=   0.00277; 
		    bs[1][9]=   0.00025; bc[1][9]=   0.00128; 
		    bs[2][1]=  -0.87747; bc[2][1]=   1.04468; 
		    bs[2][3]=   0.03455; bc[2][3]=   0.02686; 
		    bs[2][5]=  -0.00066; bc[2][5]=  -0.00379; 
		    bs[2][7]=   0.00267; bc[2][7]=   0.00043; 
		    bs[2][9]=  -0.00016; bc[2][9]=  -0.00127; 
		    bs[3][1]=  -0.45255; bc[3][1]=  -0.05603; 
		    bs[3][3]=  -0.03455; bc[3][3]=  -0.02686; 
		    bs[3][5]=   0.00066; bc[3][5]=   0.00379; 
		    bs[3][7]=  -0.00304; bc[3][7]=  -0.00025; 
		    bs[3][9]=  -0.00024; bc[3][9]=  -0.00137; 
		} else
		    /*
		if (name.equals("Test1")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    M=80; N=4; K=12;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
		    x[0]=   -1.01446419647856; y[0]=   -0.21907960526831; 
		    x[1]=    1.06293519241817; y[1]=   -0.00659066238804; 
		    x[2]=    0.03879701740017; y[2]=   -0.82631353811970; 
		    x[3]=   -0.08726786991278; y[3]=    1.05198377927308; 
		    vx[0]=    0.78449213652644; vy[0]=   -0.01932608782880; 
		    vx[1]=   -0.06332385039542; vy[1]=    1.03931139418777; 
		    vx[2]=    0.31062039501190; vy[2]=   -0.92276247111238; 
		    vx[3]=   -1.03178867006676; vy[3]=   -0.09722282870034; 
		    as[0][1]=   0.80840; ac[0][1]=  -1.00424; 
		    as[0][3]=  -0.01754; ac[0][3]=  -0.01246; 
		    as[0][5]=  -0.00018; ac[0][5]=   0.00206; 
		    as[1][1]=   0.00961; ac[1][1]=   1.06451; 
		    as[1][3]=  -0.00354; ac[1][3]=   0.00164; 
		    as[1][5]=  -0.00676; ac[1][5]=  -0.00286; 
		    as[2][1]=   0.24309; ac[2][1]=   0.02551; 
		    as[2][3]=   0.02337; ac[2][3]=   0.01397; 
		    as[2][5]=  -0.00016; ac[2][5]=  -0.00104; 
		    as[3][1]=  -1.06109; ac[3][1]=  -0.08578; 
		    as[3][3]=  -0.00229; ac[3][3]=  -0.00316; 
		    as[3][5]=   0.00711; ac[3][5]=   0.00184; 
		    bs[0][1]=   0.02551; bc[0][1]=  -0.24309; 
		    bs[0][3]=  -0.01397; bc[0][3]=   0.02337; 
		    bs[0][5]=  -0.00104; bc[0][5]=   0.00016; 
		    bs[1][1]=   1.06451; bc[1][1]=  -0.00961; 
		    bs[1][3]=  -0.00164; bc[1][3]=  -0.00354; 
		    bs[1][5]=  -0.00286; bc[1][5]=   0.00676; 
		    bs[2][1]=  -1.00424; bc[2][1]=  -0.80840; 
		    bs[2][3]=   0.01246; bc[2][3]=  -0.01754; 
		    bs[2][5]=   0.00206; bc[2][5]=   0.00018; 
		    bs[3][1]=  -0.08578; bc[3][1]=   1.06109; 
		    bs[3][3]=   0.00316; bc[3][3]=  -0.00229; 
		    bs[3][5]=   0.00184; bc[3][5]=  -0.00711; 
		} else
		    */
		if (name.equals("Hill3_0.5")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    M=80; N=3; K=32;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
		    /*
x[0]=   -0.66221014523924; y[0]=    0.00000000000958; 
x[1]=   -0.18104205513439; y[1]=   -0.00000000000573; 
x[2]=    0.84325196486745; y[2]=   -0.00000000000385; 
vx[0]=    0.07436813439526; vy[0]=   -1.42067395234400; 
vx[1]=   -0.05255540739883; vy[1]=    1.51224086533532; 
vx[2]=   -0.02181272329741; vy[2]=   -0.09156667752388; 
as[0][2]=   0.00000; ac[0][2]=  -0.62857; 
as[0][4]=   0.00000; ac[0][4]=  -0.00964; 
as[0][8]=   0.00000; ac[0][8]=  -0.02054; 
as[0][10]=   0.00000; ac[0][10]=  -0.00036; 
as[0][14]=   0.00000; ac[0][14]=  -0.00253; 
as[0][20]=  -0.00000; ac[0][20]=  -0.00043; 
as[1][1]=   0.00000; ac[1][1]=  -0.55418; 
as[1][2]=   0.00000; ac[1][2]=   0.31429; 
as[1][4]=   0.00000; ac[1][4]=   0.00482; 
as[1][5]=   0.00000; ac[1][5]=   0.03898; 
as[1][7]=   0.00000; ac[1][7]=  -0.00090; 
as[1][8]=  -0.00000; ac[1][8]=   0.01027; 
as[1][10]=   0.00000; ac[1][10]=   0.00018; 
as[1][11]=   0.00000; ac[1][11]=   0.00337; 
as[1][14]=  -0.00000; ac[1][14]=   0.00127; 
as[1][17]=   0.00000; ac[1][17]=   0.00051; 
as[1][20]=  -0.00000; ac[1][20]=   0.00022; 
as[2][1]=   0.00000; ac[2][1]=   0.55418; 
as[2][2]=  -0.00000; ac[2][2]=   0.31429; 
as[2][4]=  -0.00000; ac[2][4]=   0.00482; 
as[2][5]=   0.00000; ac[2][5]=  -0.03898; 
as[2][7]=   0.00000; ac[2][7]=   0.00090; 
as[2][8]=   0.00000; ac[2][8]=   0.01027; 
as[2][10]=  -0.00000; ac[2][10]=   0.00018; 
as[2][11]=   0.00000; ac[2][11]=  -0.00337; 
as[2][14]=   0.00000; ac[2][14]=   0.00127; 
as[2][17]=   0.00000; ac[2][17]=  -0.00051; 
as[2][20]=   0.00000; ac[2][20]=   0.00022; 
bs[0][2]=  -0.62857; bc[0][2]=  -0.00000; 
bs[0][4]=   0.00964; bc[0][4]=   0.00000; 
bs[0][8]=  -0.02054; bc[0][8]=  -0.00000; 
bs[0][10]=   0.00036; bc[0][10]=   0.00000; 
bs[0][14]=  -0.00253; bc[0][14]=  -0.00000; 
bs[0][20]=  -0.00043; bc[0][20]=  -0.00000; 
bs[1][1]=   0.55418; bc[1][1]=  -0.00000; 
bs[1][2]=   0.31429; bc[1][2]=   0.00000; 
bs[1][4]=  -0.00482; bc[1][4]=  -0.00000; 
bs[1][5]=   0.03898; bc[1][5]=   0.00000; 
bs[1][7]=   0.00090; bc[1][7]=  -0.00000; 
bs[1][8]=   0.01027; bc[1][8]=  -0.00000; 
bs[1][10]=  -0.00018; bc[1][10]=  -0.00000; 
bs[1][11]=   0.00337; bc[1][11]=   0.00000; 
bs[1][14]=   0.00127; bc[1][14]=  -0.00000; 
bs[1][17]=   0.00051; bc[1][17]=   0.00000; 
bs[1][20]=   0.00022; bc[1][20]=  -0.00000; 
bs[2][1]=  -0.55418; bc[2][1]=  -0.00000; 
bs[2][2]=   0.31429; bc[2][2]=  -0.00000; 
bs[2][4]=  -0.00482; bc[2][4]=   0.00000; 
bs[2][5]=  -0.03898; bc[2][5]=   0.00000; 
bs[2][7]=  -0.00090; bc[2][7]=  -0.00000; 
bs[2][8]=   0.01027; bc[2][8]=   0.00000; 
bs[2][10]=  -0.00018; bc[2][10]=   0.00000; 
bs[2][11]=  -0.00337; bc[2][11]=   0.00000; 
bs[2][14]=   0.00127; bc[2][14]=   0.00000; 
bs[2][17]=  -0.00051; bc[2][17]=   0.00000; 
bs[2][20]=   0.00022; bc[2][20]=   0.00000; 
*/
x[0]=   -0.66223207933681; y[0]=    0.00000000000305; 
x[1]=   -0.18102759775194; y[1]=   -0.00000000000183; 
x[2]=    0.84325862470966; y[2]=   -0.00000000000122; 
vx[0]=    0.07475415806299; vy[0]=   -1.42135918612408; 
vx[1]=   -0.05284701840503; vy[1]=    1.51275550361619; 
vx[2]=   -0.02190712312858; vy[2]=   -0.09139526528612; 
as[0][2]=   0.00000; ac[0][2]=  -0.62857; 
as[0][4]=  -0.00000; ac[0][4]=  -0.00964; 
as[0][8]=  -0.00000; ac[0][8]=  -0.02054; 
as[0][10]=   0.00000; ac[0][10]=  -0.00036; 
as[0][14]=  -0.00000; ac[0][14]=  -0.00253; 
as[0][20]=  -0.00000; ac[0][20]=  -0.00043; 
as[1][1]=   0.00000; ac[1][1]=  -0.55419; 
as[1][2]=   0.00000; ac[1][2]=   0.31429; 
as[1][4]=   0.00000; ac[1][4]=   0.00482; 
as[1][5]=   0.00000; ac[1][5]=   0.03898; 
as[1][7]=   0.00000; ac[1][7]=  -0.00090; 
as[1][8]=  -0.00000; ac[1][8]=   0.01027; 
as[1][10]=   0.00000; ac[1][10]=   0.00018; 
as[1][11]=   0.00000; ac[1][11]=   0.00337; 
as[1][14]=  -0.00000; ac[1][14]=   0.00127; 
as[1][17]=   0.00000; ac[1][17]=   0.00051; 
as[1][20]=  -0.00000; ac[1][20]=   0.00022; 
as[2][1]=   0.00000; ac[2][1]=   0.55418; 
as[2][2]=  -0.00000; ac[2][2]=   0.31429; 
as[2][4]=  -0.00000; ac[2][4]=   0.00482; 
as[2][5]=   0.00000; ac[2][5]=  -0.03898; 
as[2][7]=   0.00000; ac[2][7]=   0.00090; 
as[2][8]=   0.00000; ac[2][8]=   0.01027; 
as[2][10]=  -0.00000; ac[2][10]=   0.00018; 
as[2][11]=   0.00000; ac[2][11]=  -0.00337; 
as[2][14]=   0.00000; ac[2][14]=   0.00127; 
as[2][17]=   0.00000; ac[2][17]=  -0.00051; 
as[2][20]=   0.00000; ac[2][20]=   0.00022; 
bs[0][2]=  -0.62857; bc[0][2]=  -0.00000; 
bs[0][4]=   0.00964; bc[0][4]=  -0.00000; 
bs[0][8]=  -0.02054; bc[0][8]=  -0.00000; 
bs[0][10]=   0.00036; bc[0][10]=   0.00000; 
bs[0][14]=  -0.00253; bc[0][14]=   0.00000; 
bs[0][20]=  -0.00043; bc[0][20]=   0.00000; 
bs[1][1]=   0.55419; bc[1][1]=  -0.00000; 
bs[1][2]=   0.31429; bc[1][2]=   0.00000; 
bs[1][4]=  -0.00482; bc[1][4]=  -0.00000; 
bs[1][5]=   0.03898; bc[1][5]=   0.00000; 
bs[1][7]=   0.00090; bc[1][7]=  -0.00000; 
bs[1][8]=   0.01027; bc[1][8]=  -0.00000; 
bs[1][10]=  -0.00018; bc[1][10]=  -0.00000; 
bs[1][11]=   0.00337; bc[1][11]=   0.00000; 
bs[1][14]=   0.00127; bc[1][14]=  -0.00000; 
bs[1][17]=   0.00051; bc[1][17]=   0.00000; 
bs[1][20]=   0.00022; bc[1][20]=  -0.00000; 
bs[2][1]=  -0.55418; bc[2][1]=  -0.00000; 
bs[2][2]=   0.31429; bc[2][2]=  -0.00000; 
bs[2][4]=  -0.00482; bc[2][4]=   0.00000; 
bs[2][5]=  -0.03898; bc[2][5]=   0.00000; 
bs[2][7]=  -0.00090; bc[2][7]=  -0.00000; 
bs[2][8]=   0.01027; bc[2][8]=   0.00000; 
bs[2][10]=  -0.00018; bc[2][10]=   0.00000; 
bs[2][11]=  -0.00337; bc[2][11]=   0.00000; 
bs[2][14]=   0.00127; bc[2][14]=   0.00000; 
bs[2][17]=  -0.00051; bc[2][17]=   0.00000; 
bs[2][20]=   0.00022; bc[2][20]=   0.00000; 
		} else
		if (name.equals("Hill3_3")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    M=80; N=3; K=32;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
x[0]=   -0.97225651599086; y[0]=   -0.00000000005484; 
x[1]=    0.18954134893937; y[1]=    0.00000000273666; 
x[2]=    0.78271516695564; y[2]=   -0.00000000268182; 
vx[0]=    0.03347471765456; vy[0]=   -0.97682648430202; 
vx[1]=    0.06601011090532; vy[1]=    1.42032894168400; 
vx[2]=   -0.09948482853364; vy[2]=   -0.44350245710859; 
as[0][1]=  -0.00000; ac[0][1]=  -0.97099; 
as[0][7]=   0.00000; ac[0][7]=  -0.00029; 
as[0][9]=   0.00000; ac[0][9]=  -0.00096; 
as[1][1]=   0.00000; ac[1][1]=   0.48549; 
as[1][3]=  -0.00000; ac[1][3]=  -0.30139; 
as[1][5]=   0.00000; ac[1][5]=   0.00535; 
as[1][7]=  -0.00000; ac[1][7]=   0.00014; 
as[1][9]=   0.00000; ac[1][9]=   0.00048; 
as[1][11]=  -0.00000; ac[1][11]=  -0.00060; 
as[2][1]=  -0.00000; ac[2][1]=   0.48549; 
as[2][3]=   0.00000; ac[2][3]=   0.30139; 
as[2][5]=  -0.00000; ac[2][5]=  -0.00535; 
as[2][7]=  -0.00000; ac[2][7]=   0.00014; 
as[2][9]=  -0.00000; ac[2][9]=   0.00048; 
as[2][11]=  -0.00000; ac[2][11]=   0.00060; 
bs[0][1]=  -0.97099; bc[0][1]=  -0.00000; 
bs[0][7]=   0.00029; bc[0][7]=  -0.00000; 
bs[0][9]=  -0.00096; bc[0][9]=   0.00000; 
bs[1][1]=   0.48549; bc[1][1]=   0.00000; 
bs[1][3]=   0.30139; bc[1][3]=   0.00000; 
bs[1][5]=   0.00535; bc[1][5]=   0.00000; 
bs[1][7]=  -0.00014; bc[1][7]=   0.00000; 
bs[1][9]=   0.00048; bc[1][9]=   0.00000; 
bs[1][11]=   0.00060; bc[1][11]=   0.00000; 
bs[2][1]=   0.48549; bc[2][1]=  -0.00000; 
bs[2][3]=  -0.30139; bc[2][3]=  -0.00000; 
bs[2][5]=  -0.00535; bc[2][5]=  -0.00000; 
bs[2][7]=  -0.00014; bc[2][7]=   0.00000; 
bs[2][9]=   0.00048; bc[2][9]=  -0.00000; 
bs[2][11]=  -0.00060; bc[2][11]=   0.00000; 
		} else

	        /* 3 equal masses, sun, earth, and moon.
		 * Earth/moon orbit each other twice per year.
		 */
		if (name.equals("Hill3_2")) { // BorderCollie is Hill3_1
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    M=80; N=3; K=18;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
x[0]=   -0.97979191100525; y[0]=    0.00000000019830; 
x[1]=    0.10799162884718; y[1]=   -0.00000003618930; 
x[2]=    0.87180028215768; y[2]=    0.00000003599100; 
vx[0]=    0.03567419178898; vy[0]=   -0.98972768995942; 
vx[1]=    0.02737497524415; vy[1]=    1.33042789782613; 
vx[2]=   -0.06304916703308; vy[2]=   -0.34070020786593; 
as[0][1]=  -0.00000; ac[0][1]=  -0.97625; 
as[0][5]=   0.00000; ac[0][5]=  -0.00086; 
as[0][7]=  -0.00000; ac[0][7]=  -0.00260; 
as[1][1]=  -0.00000; ac[1][1]=   0.48813; 
as[1][2]=  -0.00000; ac[1][2]=  -0.39208; 
as[1][4]=   0.00000; ac[1][4]=   0.01102; 
as[1][5]=   0.00000; ac[1][5]=   0.00043; 
as[1][7]=  -0.00000; ac[1][7]=   0.00130; 
as[1][8]=   0.00000; ac[1][8]=  -0.00102; 
as[1][10]=   0.00000; ac[1][10]=   0.00020; 
as[2][1]=   0.00000; ac[2][1]=   0.48813; 
as[2][2]=  -0.00000; ac[2][2]=   0.39208; 
as[2][4]=   0.00000; ac[2][4]=  -0.01102; 
as[2][5]=  -0.00000; ac[2][5]=   0.00043; 
as[2][7]=   0.00000; ac[2][7]=   0.00130; 
as[2][8]=   0.00000; ac[2][8]=   0.00102; 
as[2][10]=   0.00000; ac[2][10]=  -0.00020; 
bs[0][1]=  -0.97625; bc[0][1]=   0.00000; 
bs[0][5]=   0.00086; bc[0][5]=   0.00000; 
bs[0][7]=  -0.00260; bc[0][7]=   0.00000; 
bs[1][1]=   0.48813; bc[1][1]=  -0.00000; 
bs[1][2]=   0.39208; bc[1][2]=   0.00000; 
bs[1][4]=   0.01102; bc[1][4]=   0.00000; 
bs[1][5]=  -0.00043; bc[1][5]=  -0.00000; 
bs[1][7]=   0.00130; bc[1][7]=  -0.00000; 
bs[1][8]=   0.00102; bc[1][8]=  -0.00000; 
bs[1][10]=   0.00020; bc[1][10]=   0.00000; 
bs[2][1]=   0.48813; bc[2][1]=   0.00000; 
bs[2][2]=  -0.39208; bc[2][2]=   0.00000; 
bs[2][4]=  -0.01102; bc[2][4]=   0.00000; 
bs[2][5]=  -0.00043; bc[2][5]=   0.00000; 
bs[2][7]=   0.00130; bc[2][7]=   0.00000; 
bs[2][8]=  -0.00102; bc[2][8]=  -0.00000; 
bs[2][10]=  -0.00020; bc[2][10]=   0.00000; 
		} else
		if (name.equals("Hill15")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    M=300; N=3; K=22;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
		    x[0]=    0.58501847204871; y[0]=   -0.00000000789290; 
		    x[1]=    0.37772268396161; y[1]=    0.00000000788560; 
		    x[2]=   -0.96274754536021; y[2]=    0.00000000000730; 
		    vx[0]=   -0.24736039203787; vy[0]=    2.01572248105017; 
		    vx[1]=    0.23727895876078; vy[1]=   -1.05305170775055; 
		    vx[2]=    0.01008150018376; vy[2]=   -0.96267716218242; 
		    as[0][1]=   0.00000; ac[0][1]=   0.48137; 
		    as[0][13]=  -0.00000; ac[0][13]=  -0.00019; 
		    as[0][14]=  -0.00000; ac[0][14]=  -0.00000; 
		    as[0][15]=   0.00000; ac[0][15]=   0.10383; 
		    as[1][1]=   0.00000; ac[1][1]=   0.48137; 
		    as[1][13]=   0.00000; ac[1][13]=   0.00019; 
		    as[1][14]=   0.00000; ac[1][14]=  -0.00000; 
		    as[1][15]=  -0.00000; ac[1][15]=  -0.10383; 
		    as[2][1]=  -0.00000; ac[2][1]=  -0.96275; 
		    bs[0][1]=   0.48137; bc[0][1]=  -0.00000; 
		    bs[0][15]=   0.10383; bc[0][15]=  -0.00000; 
		    bs[1][1]=   0.48137; bc[1][1]=  -0.00000; 
		    bs[1][13]=  -0.00019; bc[1][13]=   0.00000; 
		    bs[1][15]=  -0.10383; bc[1][15]=   0.00000; 
		    bs[2][1]=  -0.96275; bc[2][1]=   0.00000; 
		} else
		if (name.equals("DoubleDouble5")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    M=300; N=4; K=22;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
		    x[0]=    1.01596309261827; y[0]=   -0.00000000008221; 
		    x[1]=    0.58657173154520; y[1]=    0.00000000008216; 
		    x[2]=   -0.58657173154519; y[2]=    0.00000000008221; 
		    x[3]=   -1.01596309261827; y[3]=   -0.00000000008216; 
		    vx[0]=   -0.01937801962993; vy[0]=    1.87456397124724; 
		    vx[1]=    0.01434352456210; vy[1]=   -0.27203969131947; 
		    vx[2]=   -0.01434352538380; vy[2]=    0.27203969130656; 
		    vx[3]=    0.01937802045163; vy[3]=   -1.87456397123433; 
		    as[0][1]=  -0.00000; ac[0][1]=   0.80127; 
		    as[0][5]=   0.00000; ac[0][5]=   0.21470; 
		    as[1][1]=  -0.00000; ac[1][1]=   0.80127; 
		    as[1][5]=  -0.00000; ac[1][5]=  -0.21470; 
		    as[2][1]=   0.00000; ac[2][1]=  -0.80127; 
		    as[2][5]=  -0.00000; ac[2][5]=   0.21470; 
		    as[3][1]=   0.00000; ac[3][1]=  -0.80127; 
		    as[3][5]=   0.00000; ac[3][5]=  -0.21470; 
		    bs[0][1]=   0.80127; bc[0][1]=   0.00000; 
		    bs[0][5]=   0.21470; bc[0][5]=  -0.00000; 
		    bs[1][1]=   0.80127; bc[1][1]=   0.00000; 
		    bs[1][5]=  -0.21470; bc[1][5]=   0.00000; 
		    bs[2][1]=  -0.80127; bc[2][1]=  -0.00000; 
		    bs[2][5]=   0.21470; bc[2][5]=   0.00000; 
		    bs[3][1]=  -0.80127; bc[3][1]=  -0.00000; 
		    bs[3][5]=  -0.21470; bc[3][5]=  -0.00000; 
		} else
		if (name.equals("DoubleDouble10")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    M=300; N=4; K=22;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
		    x[0]=    0.93226647492271; y[0]=   -0.00000000003687; 
		    x[1]=    0.66102587375825; y[1]=    0.00000000003677; 
		    x[2]=   -0.66102587375825; y[2]=    0.00000000003691; 
		    x[3]=   -0.93226647492271; y[3]=   -0.00000000003682; 
		    vx[0]=   -0.04509508807505; vy[0]=    2.15195176877809; 
		    vx[1]=    0.04008962900481; vy[1]=   -0.55866990351989; 
		    vx[2]=   -0.04008962974117; vy[2]=    0.55866990349675; 
		    vx[3]=    0.04509508881141; vy[3]=   -2.15195176875495; 
		    ac[0][0]=   0.00000; bc[0][0]=   0.00000; 
		    ac[1][0]=   0.00000; bc[1][0]=   0.00000; 
		    ac[2][0]=   0.00000; bc[2][0]=   0.00000; 
		    ac[3][0]=   0.00000; bc[3][0]=   0.00000; 
		    as[0][1]=   0.00000; ac[0][1]=   0.79665; 
		    as[0][10]=   0.00000; ac[0][10]=   0.13562; 
		    as[1][1]=   0.00000; ac[1][1]=   0.79665; 
		    as[1][10]=  -0.00000; ac[1][10]=  -0.13562; 
		    as[2][1]=  -0.00000; ac[2][1]=  -0.79665; 
		    as[2][10]=  -0.00000; ac[2][10]=   0.13562; 
		    as[3][1]=  -0.00000; ac[3][1]=  -0.79665; 
		    as[3][10]=   0.00000; ac[3][10]=  -0.13562; 
		    bs[0][1]=   0.79665; bc[0][1]=  -0.00000; 
		    bs[0][10]=   0.13562; bc[0][10]=  -0.00000; 
		    bs[1][1]=   0.79665; bc[1][1]=  -0.00000; 
		    bs[1][10]=  -0.13562; bc[1][10]=   0.00000; 
		    bs[2][1]=  -0.79665; bc[2][1]=   0.00000; 
		    bs[2][10]=   0.13562; bc[2][10]=   0.00000; 
		    bs[3][1]=  -0.79665; bc[3][1]=   0.00000; 
		    bs[3][10]=  -0.13562; bc[3][10]=  -0.00000; 
		} else
		if (name.equals("DoubleDouble20")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    M=300; N=4; K=22;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
		    x[0]=    0.88037908528378; y[0]=   -0.00000000001214; 
		    x[1]=    0.70934244056923; y[1]=    0.00000000001221; 
		    x[2]=   -0.70934244056923; y[2]=    0.00000000001187; 
		    x[3]=   -0.88037908528379; y[3]=   -0.00000000001193; 
		    vx[0]=   -0.10982126907957; vy[0]=    2.50072402897539; 
		    vx[1]=    0.10482702804303; vy[1]=   -0.91101296305011; 
		    vx[2]=   -0.10482702852327; vy[2]=    0.91101296301989; 
		    vx[3]=    0.10982126955981; vy[3]=   -2.50072402894518; 
		    as[0][1]=  -0.00000; ac[0][1]=   0.79486; 
		    as[0][20]=   0.00000; ac[0][20]=   0.08552; 
		    as[1][1]=  -0.00000; ac[1][1]=   0.79486; 
		    as[1][20]=  -0.00000; ac[1][20]=  -0.08552; 
		    as[2][1]=   0.00000; ac[2][1]=  -0.79486; 
		    as[2][20]=  -0.00000; ac[2][20]=   0.08552; 
		    as[3][1]=   0.00000; ac[3][1]=  -0.79486; 
		    as[3][20]=   0.00000; ac[3][20]=  -0.08552; 
		    bs[0][1]=   0.79486; bc[0][1]=   0.00000; 
		    bs[0][20]=   0.08552; bc[0][20]=  -0.00000; 
		    bs[1][1]=   0.79486; bc[1][1]=   0.00000; 
		    bs[1][20]=  -0.08552; bc[1][20]=   0.00000; 
		    bs[2][1]=  -0.79486; bc[2][1]=  -0.00000; 
		    bs[2][20]=   0.08552; bc[2][20]=   0.00000; 
		    bs[3][1]=  -0.79486; bc[3][1]=  -0.00000; 
		    bs[3][20]=  -0.08552; bc[3][20]=  -0.00000; 
		} else
		if (name.equals("BorderCollie4")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    M=80; N=4; K=12;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
		    x[0]=   -1.25697252090716; y[0]=    0.02509889235199; 
		    x[1]=   -0.34268241269845; y[1]=   -0.04637708923925; 
		    x[2]=    0.77080325104245; y[2]=    0.72882686022197; 
		    x[3]=    0.82885183279544; y[3]=   -0.70754866333471; 
		    /*
		    vx[0]=    0.01563776934108; vy[0]=   -1.28457164569553; 
		    vx[1]=   -0.01625007172457; vy[1]=    1.35122919476817; 
		    vx[2]=    0.65272962947865; vy[2]=   -0.02562047572471; 
		    vx[3]=   -0.65211731816480; vy[3]=   -0.04103693686309; 
		    */
		    vx[0]=   -0.01563776934108; vy[0]=    1.28457164569553; 
		    vx[1]=    0.01625007172457; vy[1]=   -1.35122919476817; 
		    vx[2]=   -0.65272962947865; vy[2]=    0.02562047572471; 
		    vx[3]=    0.65211731816480; vy[3]=    0.04103693686309; 
		    as[0][1]=  -0.02410; ac[0][1]=  -1.25002; 
		    as[0][5]=  -0.00022; ac[0][5]=  -0.00129; 
		    as[0][7]=  -0.00114; ac[0][7]=  -0.00543; 
		    as[1][1]=  -0.02669; ac[1][1]=  -0.38534; 
		    as[1][3]=   0.00280; ac[1][3]=   0.03317; 
		    as[1][5]=   0.00101; ac[1][5]=   0.00690; 
		    as[1][7]=   0.00035; ac[1][7]=   0.00169; 
		    as[1][9]=   0.00017; ac[1][9]=   0.00062; 
		    as[2][1]=   0.71996; ac[2][1]=   0.78761; 
		    as[2][3]=  -0.02439; ac[2][3]=  -0.01477; 
		    as[2][5]=   0.00520; ac[2][5]=  -0.00362; 
		    as[2][7]=   0.00050; ac[2][7]=   0.00185; 
		    as[2][9]=  -0.00053; ac[2][9]=  -0.00019; 
		    as[3][1]=  -0.66917; ac[3][1]=   0.84776; 
		    as[3][3]=   0.02159; ac[3][3]=  -0.01840; 
		    as[3][5]=  -0.00600; ac[3][5]=  -0.00199; 
		    as[3][7]=   0.00029; ac[3][7]=   0.00190; 
		    as[3][9]=   0.00036; ac[3][9]=  -0.00043; 
		    bs[0][1]=  -1.25002; bc[0][1]=   0.02410; 
		    bs[0][5]=   0.00129; bc[0][5]=  -0.00022; 
		    bs[0][7]=  -0.00543; bc[0][7]=   0.00114; 
		    bs[1][1]=   1.21869; bc[1][1]=  -0.04276; 
		    bs[1][3]=   0.02654; bc[1][3]=  -0.00210; 
		    bs[1][5]=   0.00603; bc[1][5]=  -0.00087; 
		    bs[1][7]=   0.00193; bc[1][7]=  -0.00041; 
		    bs[1][9]=   0.00051; bc[1][9]=  -0.00014; 
		    bs[2][1]=   0.04574; bc[2][1]=   0.70390; 
		    bs[2][3]=  -0.01085; bc[2][3]=   0.02977; 
		    bs[2][5]=  -0.00448; bc[2][5]=  -0.00506; 
		    bs[2][7]=   0.00177; bc[2][7]=  -0.00026; 
		    bs[2][9]=  -0.00011; bc[2][9]=   0.00060; 
		    bs[3][1]=  -0.01441; bc[3][1]=  -0.68524; 
		    bs[3][3]=  -0.01570; bc[3][3]=  -0.02767; 
		    bs[3][5]=  -0.00285; bc[3][5]=   0.00614; 
		    bs[3][7]=   0.00172; bc[3][7]=  -0.00047; 
		    bs[3][9]=  -0.00040; bc[3][9]=  -0.00046; 
		} else
		if (name.equals("Lagrange3a")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    M=480; N=3; K=47;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
		    x[0]=  0.087927; y[0]=  0.560990; 
		    x[1]=  0.441868; y[1]= -0.356642; 
		    x[2]= -0.529795; y[2]= -0.204348; 
		    vx[0]=  0.663505; vy[0]=  0.401219; 
		    vx[1]=  0.015714; vy[1]= -0.775222; 
		    vx[2]= -0.679219; vy[2]=  0.374002; 
		    as[0][1 ]=  0.490521   ; bc[0][ 1]= 0.490521;
		    as[0][5 ]= -0.0306453  ; bc[0][ 5]= 0.0306453;
		    as[0][7 ]=  0.0622122  ; bc[0][ 7]= 0.0622122;
		    as[0][11]=  0.00592451 ; bc[0][11]= -0.00592451;
		    as[0][13]= -0.02955    ; bc[0][13]= -0.02955;
		    as[0][17]= -0.00202766 ; bc[0][17]= 0.00202766;
		    as[0][19]=  0.0167348  ; bc[0][19]= 0.0167348;
		    as[0][23]=  0.000814572; bc[0][23]= -0.000814572;
		    as[0][25]= -0.0103064  ; bc[0][25]= -0.0103064;
		    as[0][29]= -0.000329243; bc[0][29]= 0.000329244;
		    as[0][31]=  0.00664728 ; bc[0][31]= 0.00664728;
		    as[0][35]=  0.000120832; bc[0][35]= -0.000120832;
		    as[0][37]= -0.00440981 ; bc[0][37]= -0.00440981;
		    as[0][43]=  0.00298052 ; bc[0][43]= 0.00298052;
		    ac[0][2 ]= 0.16897     ; bs[0][ 2]=   0.16897    ;
		    ac[0][4 ]= -0.100262   ; bs[0][ 4]=   0.100262   ;
		    ac[0][8 ]= -0.0118675  ; bs[0][ 8]=  -0.0118675  ;
		    ac[0][10]= 0.041535    ; bs[0][10]=  -0.041535   ;
		    ac[0][14]= 0.00334927  ; bs[0][14]=   0.00334927 ;
		    ac[0][16]= -0.0219229  ; bs[0][16]=   0.0219229  ;
		    ac[0][20]= -0.0012744  ; bs[0][20]=  -0.0012744  ;
		    ac[0][22]= 0.0130357   ; bs[0][22]=  -0.0130357  ;
		    ac[0][26]= 0.000521167 ; bs[0][26]=   0.000521167;
		    ac[0][28]= -0.00824097 ; bs[0][28]=   0.00824097 ;
		    ac[0][32]= -0.000202875; bs[0][32]=  -0.000202875;
		    ac[0][34]= 0.00539905  ; bs[0][34]=  -0.00539905 ;
		    ac[0][40]= -0.00361843 ; bs[0][40]=   0.00361843 ;
		    ac[0][46]= 0.00246312  ; bs[0][46]=  -0.00246312 ;
		    for (j=1; j<N; j++) {
			for (k=0; k<K; k++) {
			    as[j][k] = as[0][k]*Math.cos(2*Math.PI*k*j/N)
				     - ac[0][k]*Math.sin(2*Math.PI*k*j/N);
			    ac[j][k] = ac[0][k]*Math.cos(2*Math.PI*k*j/N)
				     + as[0][k]*Math.sin(2*Math.PI*k*j/N);
			    bs[j][k] = bs[0][k]*Math.cos(2*Math.PI*k*j/N)
				     - bc[0][k]*Math.sin(2*Math.PI*k*j/N);
			    bc[j][k] = bc[0][k]*Math.cos(2*Math.PI*k*j/N)
				     + bs[0][k]*Math.sin(2*Math.PI*k*j/N);
			}
		    }

		} else
		if (name.equals("Test2")) {
		    m[0] = 1/G;
		    m[1] = 3.0034495182e-6/G;
		    m[2] = 0/G;
		    m[3] = 0/G;
		    rad[0] = 2e-1;
		    rad[1] = 1e-1;
		    rad[2] = 0.5e-1;
		    rad[3] = 0.5e-1;
		    col[0]=Color.yellow;
		    col[1]=Color.blue;
		    col[2]=Color.orange;
		    col[3]=Color.orange;
		    M=80; N=4; K=0;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
		    /*
		    x[0]=   -0.00000300347161; y[0]=   -0.00000000004340; 
		    x[1]=    1.00000238505843; y[1]=    0.00000000006899; 
		    x[2]=    0.99528674238119; y[2]=    0.00000000000003; 
		    vx[0]=    0.00000001887103; vy[0]=   -0.00000300337863; 
		    vx[1]=   -0.00628311792466; vy[1]=    0.99997605916196; 
		    vx[2]=   -0.00483562497236; vy[2]=    1.03385305640020; 
		    */
		    /*
		    // this one is nice
		    x[0]=   -0.00000355696914; y[0]=   -0.00000001260530; 
		    x[1]=    1.18429430031603; y[1]=    0.00419694589742; 
		    x[2]=    1.33616573469558; y[2]=   -0.01429773312446; 
		    vx[0]=    0.00000002428329; vy[0]=   -0.00000249247082; 
		    vx[1]=   -0.00808513863552; vy[1]=    0.82986911919042; 
		    vx[2]=    0.00622544971083; vy[2]=    0.70499520041427; 
		    */
		    // this one is great
		    x[0]=   -0.00000305362226; y[0]=    0.00000000006084; 
		    x[1]=    1.01670470604894; y[1]=   -0.00002025658753; 
		    x[2]=    0.95028015381917; y[2]=   -0.00019550950867; 
		    vx[0]=    0.00000001819623; vy[0]=   -0.00000295362770; 
		    vx[1]=   -0.00605844253290; vy[1]=    0.98341147443581; 
		    vx[2]=   -0.00675986913750; vy[2]=    1.05074551611996; 
		    x[3]=    0.77119741551444; y[3]=   -0.00345190432922; 
		    vx[3]=   -0.02180242399325; vy[3]=    1.26192836811418; 
		} else
		if (name.equals("Test3")) {
		    m[0] = 0.012277471/G;
		    m[1] = (1-0.012277471)/G;
		    m[2] = 0/G;
		    rad[0] = 2e-1;
		    rad[1] = 1e-1;
		    rad[2] = 0.5e-1;
		    col[0]=Color.yellow;
		    col[1]=Color.blue;
		    col[2]=Color.orange;
		    M=80; N=3; K=0;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
		    x[0]=   -0.00000305362226; y[0]=    0.00000000006084; 
		    x[1]=    1.01670470604894; y[1]=   -0.00002025658753; 
		    x[2]=    0.95028015381917; y[2]=   -0.00019550950867; 
		    vx[0]=    0.00000001819623; vy[0]=   -0.00000295362770; 
		    vx[1]=   -0.00605844253290; vy[1]=    0.98341147443581; 
		    vx[2]=   -0.00675986913750; vy[2]=    1.05074551611996; 
		} else
		if (name.equals("L5")) {
		    m[0] = 1/G;
		    m[1] = 3.0034495182e-6/G;
		    m[2] = 0/G;
		    m[3] = 0/G;
		    rad[0] = 2e-1;
		    rad[1] = 1e-1;
		    rad[2] = 0.5e-1;
		    rad[3] = 0.5e-1;
		    col[0]=Color.yellow;
		    col[1]=Color.blue;
		    col[2]=Color.orange;
		    col[3]=Color.orange;
		    M=80; N=4; K=0;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
		    x[0]=   -0.00000305362178; y[0]=    0.00000000003167; 
		    x[1]=    1.01670470323982; y[1]=   -0.00001082901117; 
		    x[2]=    0.48715846284775; y[2]=    0.86546296502841; 
		    vx[0]=    0.00000001822494; vy[0]=   -0.00000295362906; 
		    vx[1]=   -0.00606772437697; vy[1]=    0.98341208480921; 
		    vx[2]=   -0.87002050044904; vy[2]=    0.50658221004099; 
		    x[3]=    0.48715846284775; y[3]=   -0.86546296502841; 
		    vx[3]=   0.87002050044904; vy[3]=    0.50658221004099; 
		} else
		if (name.equals("Test")) {
		    m[0] = 1.33e+11/G;
		    m[1] = 398600/G;
		    m[2] = 4978/G;
		    m[3] = 0/G;
		    rad[0] =2e+6;
		    rad[1] =1e+6;
		    rad[2] =0.5e+6;
		    rad[3] =0.25e+6;
		    col[0]=Color.yellow;
		    col[1]=Color.blue;
		    col[2]=Color.orange;
		    col[3]=Color.white;
		    M=80; N=4; K=0;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
		    x[0]=  0.000000; y[0]=  0.000000; 
		    x[1]=148800000.000000; y[1]=  0.000000; 
		    x[2]=148800000.000000; y[2]=400000.000000; 
		    x[3]=148800000.000000; y[3]=6578.000000; 
		    vx[0]=  0.000000; vy[0]=  0.000000; 
		    vx[1]= -0.000000; vy[1]= 29.646767; 
		    vx[2]= -1.039353; vy[2]= 29.646754; 
		    vx[3]= -0.735840; vy[3]=  7.263530; 
		} else
		if (name.equals("Ducati3_2")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    m[n-1] *= 2; rad[n-1] *= Math.pow(2,1.0/3.0);
		    M=80; N=3; K=40;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
	    as[0][1]=   1.402716819428644; bc[0][1]=   0.213564231020102; 
	    as[0][3]=   0.041810689051194; bc[0][3]=   0.047346035951955; 
	    as[0][5]=   0.010110200768179; bc[0][5]=   0.004558967395620; 
	    as[0][7]=   0.001473869909513; bc[0][7]=   0.001581999750768; 
	    as[0][9]=   0.000497800842500; bc[0][9]=   0.000274711917643; 
	    as[0][11]=   0.000099491918249; bc[0][11]=   0.000104597977357; 
	    as[0][13]=   0.000036661915448; bc[0][13]=   0.000022368476791; 
	    as[0][15]=   0.000008457396395; bc[0][15]=   0.000008793196337; 
	    as[0][17]=   0.000003230841857; bc[0][17]=   0.000002099576339; 
	    as[0][19]=   0.000000810512755; bc[0][19]=   0.000000836846572; 
	    as[0][21]=   0.000000315283443; bc[0][21]=   0.000000214187397; 
	    as[0][23]=   0.000000083632215; bc[0][23]=   0.000000085936777; 
	    as[0][25]=   0.000000032862902; bc[0][25]=   0.000000023083943; 
	    as[0][27]=   0.000000009074471; bc[0][27]=   0.000000009291889; 
	    as[0][29]=   0.000000003587305; bc[0][29]=   0.000000002587012; 
	    as[0][31]=   0.000000001021252; bc[0][31]=   0.000000001042921; 
	    as[0][33]=   0.000000000405227; bc[0][33]=   0.000000000298555; 
	    as[0][35]=   0.000000000118172; bc[0][35]=   0.000000000120431; 
	    as[0][37]=   0.000000000046998; bc[0][37]=   0.000000000035258; 
	    as[0][39]=   0.000000000013973; bc[0][39]=   0.000000000014227; 
	    as[1][1]=   0.213564231020102; bc[1][1]=   1.402716819428644; 
	    as[1][3]=  -0.047346035951955; bc[1][3]=  -0.041810689051194; 
	    as[1][5]=   0.004558967395620; bc[1][5]=   0.010110200768179; 
	    as[1][7]=  -0.001581999750768; bc[1][7]=  -0.001473869909513; 
	    as[1][9]=   0.000274711917643; bc[1][9]=   0.000497800842500; 
	    as[1][11]=  -0.000104597977357; bc[1][11]=  -0.000099491918249; 
	    as[1][13]=   0.000022368476791; bc[1][13]=   0.000036661915448; 
	    as[1][15]=  -0.000008793196337; bc[1][15]=  -0.000008457396395; 
	    as[1][17]=   0.000002099576339; bc[1][17]=   0.000003230841857; 
	    as[1][19]=  -0.000000836846572; bc[1][19]=  -0.000000810512755; 
	    as[1][21]=   0.000000214187397; bc[1][21]=   0.000000315283443; 
	    as[1][23]=  -0.000000085936777; bc[1][23]=  -0.000000083632215; 
	    as[1][25]=   0.000000023083943; bc[1][25]=   0.000000032862902; 
	    as[1][27]=  -0.000000009291889; bc[1][27]=  -0.000000009074471; 
	    as[1][29]=   0.000000002587012; bc[1][29]=   0.000000003587305; 
	    as[1][31]=  -0.000000001042921; bc[1][31]=  -0.000000001021251; 
	    as[1][33]=   0.000000000298555; bc[1][33]=   0.000000000405227; 
	    as[1][35]=  -0.000000000120431; bc[1][35]=  -0.000000000118172; 
	    as[1][37]=   0.000000000035258; bc[1][37]=   0.000000000046998; 
	    as[1][39]=  -0.000000000014227; bc[1][39]=  -0.000000000013973; 
	    as[2][1]=  -0.808140525224626; bc[2][1]=  -0.808140525224625; 
	    as[2][3]=   0.002767673450380; bc[2][3]=  -0.002767673450380; 
	    as[2][5]=  -0.007334584081899; bc[2][5]=  -0.007334584081899; 
	    as[2][7]=   0.000054064920627; bc[2][7]=  -0.000054064920627; 
	    as[2][9]=  -0.000386256380071; bc[2][9]=  -0.000386256380071; 
	    as[2][11]=   0.000002553029554; bc[2][11]=  -0.000002553029554; 
	    as[2][13]=  -0.000029515196119; bc[2][13]=  -0.000029515196119; 
	    as[2][15]=   0.000000167899971; bc[2][15]=  -0.000000167899971; 
	    as[2][17]=  -0.000002665209098; bc[2][17]=  -0.000002665209098; 
	    as[2][19]=   0.000000013166909; bc[2][19]=  -0.000000013166909; 
	    as[2][21]=  -0.000000264735420; bc[2][21]=  -0.000000264735420; 
	    as[2][23]=   0.000000001152281; bc[2][23]=  -0.000000001152281; 
	    as[2][25]=  -0.000000027973423; bc[2][25]=  -0.000000027973423; 
	    as[2][27]=   0.000000000108709; bc[2][27]=  -0.000000000108709; 
	    as[2][29]=  -0.000000003087159; bc[2][29]=  -0.000000003087159; 
	    as[2][31]=   0.000000000010835; bc[2][31]=  -0.000000000010835; 
	    as[2][33]=  -0.000000000351891; bc[2][33]=  -0.000000000351891; 
	    as[2][35]=   0.000000000001130; bc[2][35]=  -0.000000000001130; 
	    as[2][37]=  -0.000000000041128; bc[2][37]=  -0.000000000041128; 
	    as[2][39]=   0.000000000000127; bc[2][39]=  -0.000000000000127; 
	    x[0]=   0.000000000000000; y[0]=   0.267464978707892; 
	    x[1]=  -0.000000000000000; y[1]=   1.369971653334316; 
	    x[2]=  -0.000000000000000; y[2]=  -0.818718316021356; 
	    /*
	    vx[0]=   1.595024893244509; vy[0]=  -0.007614159935814; 
	    vx[1]=   0.084798653316563; vy[1]=  -0.010822763417961; 
	    vx[2]=  -0.839911773280789; vy[2]=   0.009218461676890; 
	    */
	    vx[0]=  -1.595024893244509; vy[0]=   0.007614159935814; 
	    vx[1]=  -0.084798653316563; vy[1]=   0.010822763417961; 
	    vx[2]=   0.839911773280789; vy[2]=  -0.009218461676890; 

		} else
		if (name.equals("Ducati3_0.5")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    m[n-1] /= 2; rad[n-1] *= Math.pow(0.5,1.0/3.0);
		    M=80; N=3; K=40;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
	    as[0][1]=   0.894087552790666; bc[0][1]=  -0.330544168677541; 
	    as[0][3]=   0.017018826977239; bc[0][3]=   0.020060381612618; 
	    as[0][5]=   0.004570582781581; bc[0][5]=   0.002614540402871; 
	    as[0][7]=   0.000840818621854; bc[0][7]=   0.000915561746402; 
	    as[0][9]=   0.000275421194616; bc[0][9]=   0.000230458052854; 
	    as[0][11]=   0.000075210038110; bc[0][11]=   0.000079678364048; 
	    as[0][13]=   0.000026304105059; bc[0][13]=   0.000024357120353; 
	    as[0][15]=   0.000008331720106; bc[0][15]=   0.000008703838454; 
	    as[0][17]=   0.000003023703343; bc[0][17]=   0.000002912471129; 
	    as[0][19]=   0.000001032206334; bc[0][19]=   0.000001069150073; 
	    as[0][21]=   0.000000383924347; bc[0][21]=   0.000000376710286; 
	    as[0][23]=   0.000000137027853; bc[0][23]=   0.000000141119644; 
	    as[0][25]=   0.000000051878057; bc[0][25]=   0.000000051401995; 
	    as[0][27]=   0.000000019067857; bc[0][27]=   0.000000019556295; 
	    as[0][29]=   0.000000007316752; bc[0][29]=   0.000000007289892; 
	    as[0][31]=   0.000000002745804; bc[0][31]=   0.000000002807401; 
	    as[0][33]=   0.000000001064838; bc[0][33]=   0.000000001064483; 
	    as[0][35]=   0.000000000405841; bc[0][35]=   0.000000000413985; 
	    as[0][37]=   0.000000000158727; bc[0][37]=   0.000000000159060; 
	    as[0][39]=   0.000000000061217; bc[0][39]=   0.000000000062379; 
	    as[1][1]=  -0.330544168677542; bc[1][1]=   0.894087552790666; 
	    as[1][3]=  -0.020060381612618; bc[1][3]=  -0.017018826977239; 
	    as[1][5]=   0.002614540402871; bc[1][5]=   0.004570582781581; 
	    as[1][7]=  -0.000915561746402; bc[1][7]=  -0.000840818621854; 
	    as[1][9]=   0.000230458052854; bc[1][9]=   0.000275421194616; 
	    as[1][11]=  -0.000079678364048; bc[1][11]=  -0.000075210038110; 
	    as[1][13]=   0.000024357120353; bc[1][13]=   0.000026304105059; 
	    as[1][15]=  -0.000008703838454; bc[1][15]=  -0.000008331720106; 
	    as[1][17]=   0.000002912471129; bc[1][17]=   0.000003023703343; 
	    as[1][19]=  -0.000001069150073; bc[1][19]=  -0.000001032206334; 
	    as[1][21]=   0.000000376710286; bc[1][21]=   0.000000383924347; 
	    as[1][23]=  -0.000000141119644; bc[1][23]=  -0.000000137027853; 
	    as[1][25]=   0.000000051401995; bc[1][25]=   0.000000051878057; 
	    as[1][27]=  -0.000000019556295; bc[1][27]=  -0.000000019067857; 
	    as[1][29]=   0.000000007289892; bc[1][29]=   0.000000007316752; 
	    as[1][31]=  -0.000000002807401; bc[1][31]=  -0.000000002745804; 
	    as[1][33]=   0.000000001064483; bc[1][33]=   0.000000001064838; 
	    as[1][35]=  -0.000000000413985; bc[1][35]=  -0.000000000405841; 
	    as[1][37]=   0.000000000159060; bc[1][37]=   0.000000000158727; 
	    as[1][39]=  -0.000000000062379; bc[1][39]=  -0.000000000061217; 
	    as[2][1]=  -1.127086768226172; bc[2][1]=  -1.127086768226171; 
	    as[2][3]=   0.006083109270757; bc[2][3]=  -0.006083109270757; 
	    as[2][5]=  -0.014370246368902; bc[2][5]=  -0.014370246368903; 
	    as[2][7]=   0.000149486249096; bc[2][7]=  -0.000149486249096; 
	    as[2][9]=  -0.001011758494941; bc[2][9]=  -0.001011758494941; 
	    as[2][11]=   0.000008936651876; bc[2][11]=  -0.000008936651876; 
	    as[2][13]=  -0.000101322450824; bc[2][13]=  -0.000101322450824; 
	    as[2][15]=   0.000000744236696; bc[2][15]=  -0.000000744236696; 
	    as[2][17]=  -0.000011872348945; bc[2][17]=  -0.000011872348945; 
	    as[2][19]=   0.000000073887478; bc[2][19]=  -0.000000073887478; 
	    as[2][21]=  -0.000001521269267; bc[2][21]=  -0.000001521269267; 
	    as[2][23]=   0.000000008183582; bc[2][23]=  -0.000000008183582; 
	    as[2][25]=  -0.000000206560104; bc[2][25]=  -0.000000206560104; 
	    as[2][27]=   0.000000000976877; bc[2][27]=  -0.000000000976877; 
	    as[2][29]=  -0.000000029213288; bc[2][29]=  -0.000000029213288; 
	    as[2][31]=   0.000000000123194; bc[2][31]=  -0.000000000123194; 
	    as[2][33]=  -0.000000004258642; bc[2][33]=  -0.000000004258642; 
	    as[2][35]=   0.000000000016288; bc[2][35]=  -0.000000000016288; 
	    as[2][37]=  -0.000000000635575; bc[2][37]=  -0.000000000635575; 
	    as[2][39]=   0.000000000002324; bc[2][39]=  -0.000000000002324; 
	    x[0]=  -0.000000000000000; y[0]=  -0.306605905333318; 
	    x[1]=  -0.000000000000000; y[1]=   0.881018950045771; 
	    x[2]=  -0.000000000000000; y[2]=  -1.148826089424827; 
	    /*
	    vx[0]=   0.977606709756782; vy[0]=   0.000035836180345; 
	    vx[1]=  -0.382604245677866; vy[1]=  -0.007246113081741; 
	    vx[2]=  -1.190004928157757; vy[2]=   0.014420553802792; 
	    */
	    vx[0]=  -0.977606709756782; vy[0]=  -0.000035836180345; 
	    vx[1]=   0.382604245677866; vy[1]=   0.007246113081741; 
	    vx[2]=   1.190004928157757; vy[2]=  -0.014420553802792; 

		} else
		if (name.equals("Ducati3_0.1")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    m[n-1] /= 10; rad[n-1] *= Math.pow(0.1,1.0/3.0);
		    M=80; N=3; K=40;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
	    as[0][1]=   0.689366551618465; bc[0][1]=  -0.561132745088638; 
	    as[0][3]=   0.004109581016852; bc[0][3]=   0.004958008939741; 
	    as[0][5]=   0.001149587098495; bc[0][5]=   0.000726335184293; 
	    as[0][7]=   0.000233486623707; bc[0][7]=   0.000256651510264; 
	    as[0][9]=   0.000077158861235; bc[0][9]=   0.000072261676088; 
	    as[0][11]=   0.000023566254773; bc[0][11]=   0.000025115004840; 
	    as[0][13]=   0.000008453583416; bc[0][13]=   0.000008397913115; 
	    as[0][15]=   0.000002936851350; bc[0][15]=   0.000003081300118; 
	    as[0][17]=   0.000001105020875; bc[0][17]=   0.000001114770696; 
	    as[0][19]=   0.000000408830602; bc[0][19]=   0.000000424900336; 
	    as[0][21]=   0.000000158664329; bc[0][21]=   0.000000160826245; 
	    as[0][23]=   0.000000060947521; bc[0][23]=   0.000000062942482; 
	    as[0][25]=   0.000000024171498; bc[0][25]=   0.000000024533090; 
	    as[0][27]=   0.000000009520610; bc[0][27]=   0.000000009787577; 
	    as[0][29]=   0.000000003836838; bc[0][29]=   0.000000003894402; 
	    as[0][31]=   0.000000001538674; bc[0][31]=   0.000000001576427; 
	    as[0][33]=   0.000000000627793; bc[0][33]=   0.000000000636949; 
	    as[0][35]=   0.000000000255194; bc[0][35]=   0.000000000260798; 
	    as[0][37]=   0.000000000105139; bc[0][37]=   0.000000000106651; 
	    as[0][39]=   0.000000000043191; bc[0][39]=   0.000000000044085; 
	    as[1][1]=  -0.561132745088637; bc[1][1]=   0.689366551618465; 
	    as[1][3]=  -0.004958008939740; bc[1][3]=  -0.004109581016852; 
	    as[1][5]=   0.000726335184293; bc[1][5]=   0.001149587098495; 
	    as[1][7]=  -0.000256651510264; bc[1][7]=  -0.000233486623707; 
	    as[1][9]=   0.000072261676088; bc[1][9]=   0.000077158861235; 
	    as[1][11]=  -0.000025115004840; bc[1][11]=  -0.000023566254773; 
	    as[1][13]=   0.000008397913115; bc[1][13]=   0.000008453583416; 
	    as[1][15]=  -0.000003081300118; bc[1][15]=  -0.000002936851350; 
	    as[1][17]=   0.000001114770696; bc[1][17]=   0.000001105020875; 
	    as[1][19]=  -0.000000424900336; bc[1][19]=  -0.000000408830602; 
	    as[1][21]=   0.000000160826245; bc[1][21]=   0.000000158664329; 
	    as[1][23]=  -0.000000062942482; bc[1][23]=  -0.000000060947521; 
	    as[1][25]=   0.000000024533090; bc[1][25]=   0.000000024171498; 
	    as[1][27]=  -0.000000009787577; bc[1][27]=  -0.000000009520610; 
	    as[1][29]=   0.000000003894402; bc[1][29]=   0.000000003836838; 
	    as[1][31]=  -0.000000001576427; bc[1][31]=  -0.000000001538674; 
	    as[1][33]=   0.000000000636949; bc[1][33]=   0.000000000627793; 
	    as[1][35]=  -0.000000000260798; bc[1][35]=  -0.000000000255194; 
	    as[1][37]=   0.000000000106651; bc[1][37]=   0.000000000105139; 
	    as[1][39]=  -0.000000000044085; bc[1][39]=  -0.000000000043191; 
	    as[2][1]=  -1.282338065297950; bc[2][1]=  -1.282338065297950; 
	    as[2][3]=   0.008484279228886; bc[2][3]=  -0.008484279228886; 
	    as[2][5]=  -0.018759222827880; bc[2][5]=  -0.018759222827880; 
	    as[2][7]=   0.000231648865568; bc[2][7]=  -0.000231648865568; 
	    as[2][9]=  -0.001494205373233; bc[2][9]=  -0.001494205373233; 
	    as[2][11]=   0.000015487500670; bc[2][11]=  -0.000015487500670; 
	    as[2][13]=  -0.000168514965313; bc[2][13]=  -0.000168514965313; 
	    as[2][15]=   0.000001444487671; bc[2][15]=  -0.000001444487671; 
	    as[2][17]=  -0.000022197915712; bc[2][17]=  -0.000022197915712; 
	    as[2][19]=   0.000000160697336; bc[2][19]=  -0.000000160697336; 
	    as[2][21]=  -0.000003194905739; bc[2][21]=  -0.000003194905739; 
	    as[2][23]=   0.000000019949610; bc[2][23]=  -0.000000019949610; 
	    as[2][25]=  -0.000000487045879; bc[2][25]=  -0.000000487045879; 
	    as[2][27]=   0.000000002669670; bc[2][27]=  -0.000000002669670; 
	    as[2][29]=  -0.000000077312394; bc[2][29]=  -0.000000077312394; 
	    as[2][31]=   0.000000000377533; bc[2][31]=  -0.000000000377533; 
	    as[2][33]=  -0.000000012647422; bc[2][33]=  -0.000000012647422; 
	    as[2][35]=   0.000000000056040; bc[2][35]=  -0.000000000056040; 
	    as[2][37]=  -0.000000002117895; bc[2][37]=  -0.000000002117895; 
	    as[2][39]=   0.000000000008940; bc[2][39]=  -0.000000000008940; 
	    x[0]=   0.000000000000000; y[0]=  -0.555081089280441; 
	    x[1]=  -0.000000000000000; y[1]=   0.686232991705607; 
	    x[2]=  -0.000000000000000; y[2]=  -1.311519024251342; 
	    /*
	    vx[0]=   0.710159151725980; vy[0]=   0.004138380172086; 
	    vx[1]=  -0.573687800810489; vy[1]=  -0.005882145551611; 
	    vx[2]=  -1.364713509154591; vy[2]=   0.017437653795249; 
	    */
	    vx[0]=  -0.710159151725980; vy[0]=  -0.004138380172086; 
	    vx[1]=   0.573687800810489; vy[1]=   0.005882145551611; 
	    vx[2]=   1.364713509154591; vy[2]=  -0.017437653795249; 

		} else
		if (name.equals("Ducati3_10")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    m[n-1] *= 10; rad[n-1] *= Math.pow(10,1.0/3.0);
		    M=80; N=3; K=40;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
	    as[0][1]=   2.500218116871762; bc[0][1]=   1.331619454241626; 
	    as[0][3]=   0.069804699947113; bc[0][3]=   0.075460805028738; 
	    as[0][5]=   0.013381760228298; bc[0][5]=   0.003255706763703; 
	    as[0][7]=   0.001209039706950; bc[0][7]=   0.001282199048642; 
	    as[0][9]=   0.000521928011746; bc[0][9]=  -0.000030586079365; 
	    as[0][11]=   0.000046431298701; bc[0][11]=   0.000048903207620; 
	    as[0][13]=   0.000033616143996; bc[0][13]=  -0.000010486416504; 
	    as[0][15]=   0.000002434124034; bc[0][15]=   0.000002557688526; 
	    as[0][17]=   0.000002773306471; bc[0][17]=  -0.000001396210839; 
	    as[0][19]=   0.000000152291529; bc[0][19]=   0.000000160017210; 
	    as[0][21]=   0.000000264558248; bc[0][21]=  -0.000000169355830; 
	    as[0][23]=   0.000000010719762; bc[0][23]=   0.000000011279974; 
	    as[0][25]=   0.000000027683292; bc[0][25]=  -0.000000020358981; 
	    as[0][27]=   0.000000000821915; bc[0][27]=   0.000000000867105; 
	    as[0][29]=   0.000000003083182; bc[0][29]=  -0.000000002471610; 
	    as[0][31]=   0.000000000067261; bc[0][31]=   0.000000000071212; 
	    as[0][33]=   0.000000000358861; bc[0][33]=  -0.000000000304344; 
	    as[0][35]=   0.000000000005790; bc[0][35]=   0.000000000006158; 
	    as[0][37]=   0.000000000043148; bc[0][37]=  -0.000000000038021; 
	    as[0][39]=   0.000000000000513; bc[0][39]=   0.000000000000549; 
	    as[1][1]=   1.331619454241627; bc[1][1]=   2.500218116871761; 
	    as[1][3]=  -0.075460805028738; bc[1][3]=  -0.069804699947113; 
	    as[1][5]=   0.003255706763703; bc[1][5]=   0.013381760228298; 
	    as[1][7]=  -0.001282199048642; bc[1][7]=  -0.001209039706950; 
	    as[1][9]=  -0.000030586079365; bc[1][9]=   0.000521928011746; 
	    as[1][11]=  -0.000048903207620; bc[1][11]=  -0.000046431298701; 
	    as[1][13]=  -0.000010486416504; bc[1][13]=   0.000033616143996; 
	    as[1][15]=  -0.000002557688526; bc[1][15]=  -0.000002434124034; 
	    as[1][17]=  -0.000001396210839; bc[1][17]=   0.000002773306471; 
	    as[1][19]=  -0.000000160017210; bc[1][19]=  -0.000000152291529; 
	    as[1][21]=  -0.000000169355830; bc[1][21]=   0.000000264558248; 
	    as[1][23]=  -0.000000011279974; bc[1][23]=  -0.000000010719762; 
	    as[1][25]=  -0.000000020358981; bc[1][25]=   0.000000027683292; 
	    as[1][27]=  -0.000000000867105; bc[1][27]=  -0.000000000821915; 
	    as[1][29]=  -0.000000002471610; bc[1][29]=   0.000000003083182; 
	    as[1][31]=  -0.000000000071212; bc[1][31]=  -0.000000000067261; 
	    as[1][33]=  -0.000000000304344; bc[1][33]=   0.000000000358861; 
	    as[1][35]=  -0.000000000006158; bc[1][35]=  -0.000000000005790; 
	    as[1][37]=  -0.000000000038021; bc[1][37]=   0.000000000043148; 
	    as[1][39]=  -0.000000000000549; bc[1][39]=  -0.000000000000513; 
	    as[2][1]=  -0.383183757111363; bc[2][1]=  -0.383183757111363; 
	    as[2][3]=   0.000565610508162; bc[2][3]=  -0.000565610508162; 
	    as[2][5]=  -0.001663746699200; bc[2][5]=  -0.001663746699200; 
	    as[2][7]=   0.000007315934169; bc[2][7]=  -0.000007315934169; 
	    as[2][9]=  -0.000049134193238; bc[2][9]=  -0.000049134193238; 
	    as[2][11]=   0.000000247190892; bc[2][11]=  -0.000000247190892; 
	    as[2][13]=  -0.000002312972749; bc[2][13]=  -0.000002312972749; 
	    as[2][15]=   0.000000012356449; bc[2][15]=  -0.000000012356449; 
	    as[2][17]=  -0.000000137709563; bc[2][17]=  -0.000000137709563; 
	    as[2][19]=   0.000000000772568; bc[2][19]=  -0.000000000772568; 
	    as[2][21]=  -0.000000009520242; bc[2][21]=  -0.000000009520242; 
	    as[2][23]=   0.000000000056021; bc[2][23]=  -0.000000000056021; 
	    as[2][25]=  -0.000000000732431; bc[2][25]=  -0.000000000732431; 
	    as[2][27]=   0.000000000004519; bc[2][27]=  -0.000000000004519; 
	    as[2][29]=  -0.000000000061157; bc[2][29]=  -0.000000000061157; 
	    as[2][31]=   0.000000000000395; bc[2][31]=  -0.000000000000395; 
	    as[2][33]=  -0.000000000005452; bc[2][33]=  -0.000000000005452; 
	    as[2][35]=   0.000000000000037; bc[2][35]=  -0.000000000000037; 
	    as[2][37]=  -0.000000000000513; bc[2][37]=  -0.000000000000513; 
	    as[2][39]=   0.000000000000004; bc[2][39]=  -0.000000000000004; 
	    x[0]=   0.000000000000000; y[0]=   1.411627136985569; 
	    x[1]=   0.000000000000000; y[1]=   2.443095721305436; 
	    x[2]=   0.000000000000000; y[2]=  -0.385472285829125; 
	    /*
	    vx[0]=   2.790386559838153; vy[0]=  -0.018818814090859; 
	    vx[1]=   1.111565091617343; vy[1]=  -0.019109515708231; 
	    vx[2]=  -0.390195165145574; vy[2]=   0.003792832979909; 
	    */
	    vx[0]=  -2.790386559838153; vy[0]=   0.018818814090859; 
	    vx[1]=  -1.111565091617343; vy[1]=   0.019109515708231; 
	    vx[2]=   0.390195165145574; vy[2]=  -0.003792832979909; 

		} else
		if (name.equals("Ducati3_1.2")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    m[n-1] *= 1.2; rad[n-1] *= Math.pow(1.2,1.0/3.0);
		    M=80; N=3; K=40;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
	    as[0][1]=   1.167595786002363; bc[0][1]=  -0.034155345532810; 
	    as[0][3]=   0.031532534014338; bc[0][3]=   0.036287269992928; 
	    as[0][5]=   0.008009814620430; bc[0][5]=   0.004031308768587; 
	    as[0][7]=   0.001298260947658; bc[0][7]=   0.001401113829413; 
	    as[0][9]=   0.000428792925759; bc[0][9]=   0.000295037641445; 
	    as[0][11]=   0.000099604449745; bc[0][11]=   0.000104984384699; 
	    as[0][13]=   0.000035038757947; bc[0][13]=   0.000027418402705; 
	    as[0][15]=   0.000009532586192; bc[0][15]=   0.000009924071796; 
	    as[0][17]=   0.000003443601140; bc[0][17]=   0.000002888279807; 
	    as[0][19]=   0.000001023463579; bc[0][19]=   0.000001057404751; 
	    as[0][21]=   0.000000375172211; bc[0][21]=   0.000000328401676; 
	    as[0][23]=   0.000000117947076; bc[0][23]=   0.000000121228866; 
	    as[0][25]=   0.000000043647890; bc[0][25]=   0.000000039311176; 
	    as[0][27]=   0.000000014263121; bc[0][27]=   0.000000014605043; 
	    as[0][29]=   0.000000005314391; bc[0][29]=   0.000000004883567; 
	    as[0][31]=   0.000000001786194; bc[0][31]=   0.000000001823819; 
	    as[0][33]=   0.000000000669042; bc[0][33]=   0.000000000623962; 
	    as[0][35]=   0.000000000229715; bc[0][35]=   0.000000000234049; 
	    as[0][37]=   0.000000000086403; bc[0][37]=   0.000000000081514; 
	    as[0][39]=   0.000000000030160; bc[0][39]=   0.000000000030700; 
	    as[1][1]=  -0.034155345532810; bc[1][1]=   1.167595786002364; 
	    as[1][3]=  -0.036287269992928; bc[1][3]=  -0.031532534014338; 
	    as[1][5]=   0.004031308768587; bc[1][5]=   0.008009814620430; 
	    as[1][7]=  -0.001401113829413; bc[1][7]=  -0.001298260947658; 
	    as[1][9]=   0.000295037641445; bc[1][9]=   0.000428792925759; 
	    as[1][11]=  -0.000104984384699; bc[1][11]=  -0.000099604449745; 
	    as[1][13]=   0.000027418402705; bc[1][13]=   0.000035038757947; 
	    as[1][15]=  -0.000009924071796; bc[1][15]=  -0.000009532586192; 
	    as[1][17]=   0.000002888279807; bc[1][17]=   0.000003443601140; 
	    as[1][19]=  -0.000001057404751; bc[1][19]=  -0.000001023463579; 
	    as[1][21]=   0.000000328401676; bc[1][21]=   0.000000375172211; 
	    as[1][23]=  -0.000000121228866; bc[1][23]=  -0.000000117947076; 
	    as[1][25]=   0.000000039311176; bc[1][25]=   0.000000043647890; 
	    as[1][27]=  -0.000000014605043; bc[1][27]=  -0.000000014263121; 
	    as[1][29]=   0.000000004883567; bc[1][29]=   0.000000005314391; 
	    as[1][31]=  -0.000000001823819; bc[1][31]=  -0.000000001786194; 
	    as[1][33]=   0.000000000623962; bc[1][33]=   0.000000000669042; 
	    as[1][35]=  -0.000000000234049; bc[1][35]=  -0.000000000229715; 
	    as[1][37]=   0.000000000081514; bc[1][37]=   0.000000000086403; 
	    as[1][39]=  -0.000000000030700; bc[1][39]=  -0.000000000030160; 
	    as[2][1]=  -0.944533700391305; bc[2][1]=  -0.944533700391304; 
	    as[2][3]=   0.003962279982159; bc[2][3]=  -0.003962279982159; 
	    as[2][5]=  -0.010034269490848; bc[2][5]=  -0.010034269490848; 
	    as[2][7]=   0.000085710734796; bc[2][7]=  -0.000085710734796; 
	    as[2][9]=  -0.000603192139337; bc[2][9]=  -0.000603192139337; 
	    as[2][11]=   0.000004483279128; bc[2][11]=  -0.000004483279128; 
	    as[2][13]=  -0.000052047633876; bc[2][13]=  -0.000052047633876; 
	    as[2][15]=   0.000000326238004; bc[2][15]=  -0.000000326238004; 
	    as[2][17]=  -0.000005276567456; bc[2][17]=  -0.000005276567456; 
	    as[2][19]=   0.000000028284309; bc[2][19]=  -0.000000028284309; 
	    as[2][21]=  -0.000000586311573; bc[2][21]=  -0.000000586311573; 
	    as[2][23]=   0.000000002734825; bc[2][23]=  -0.000000002734825; 
	    as[2][25]=  -0.000000069132555; bc[2][25]=  -0.000000069132555; 
	    as[2][27]=   0.000000000284934; bc[2][27]=  -0.000000000284934; 
	    as[2][29]=  -0.000000008498298; bc[2][29]=  -0.000000008498298; 
	    as[2][31]=   0.000000000031354; bc[2][31]=  -0.000000000031354; 
	    as[2][33]=  -0.000000001077504; bc[2][33]=  -0.000000001077504; 
	    as[2][35]=   0.000000000003612; bc[2][35]=  -0.000000000003612; 
	    as[2][37]=  -0.000000000139931; bc[2][37]=  -0.000000000139931; 
	    as[2][39]=   0.000000000000450; bc[2][39]=  -0.000000000000450; 
	    x[0]=   0.000000000000000; y[0]=   0.008006168467695; 
	    x[1]=   0.000000000000000; y[1]=   1.143132211079798; 
	    x[2]=   0.000000000000000; y[2]=  -0.959281982956253; 
	    /*
	    vx[0]=   1.316768353837866; vy[0]=  -0.004419577531932; 
	    vx[1]=  -0.130881869125341; vy[1]=  -0.009141805952915; 
	    vx[2]=  -0.988238737260448; vy[2]=   0.011301152904040; 
	    */
	    vx[0]=  -1.316768353837866; vy[0]=   0.004419577531932; 
	    vx[1]=   0.130881869125341; vy[1]=   0.009141805952915; 
	    vx[2]=   0.988238737260448; vy[2]=  -0.011301152904040; 
		} else
		if (name.equals("Ducati3_1.3")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    m[n-1] *= 1.3; rad[n-1] *= Math.pow(1.3,1.0/3.0);
		    M=80; N=3; K=40;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
	    as[0][1]=   1.200607556358488; bc[0][1]=   0.000953088124940; 
	    as[0][3]=   0.033091761653644; bc[0][3]=   0.037987408256145; 
	    as[0][5]=   0.008348586912074; bc[0][5]=   0.004137844595024; 
	    as[0][7]=   0.001333022505006; bc[0][7]=   0.001437381205791; 
	    as[0][9]=   0.000441308286451; bc[0][9]=   0.000295358255962; 
	    as[0][11]=   0.000100433041677; bc[0][11]=   0.000105810296229; 
	    as[0][13]=   0.000035475956287; bc[0][13]=   0.000026992872406; 
	    as[0][15]=   0.000009449870250; bc[0][15]=   0.000009835316246; 
	    as[0][17]=   0.000003427491382; bc[0][17]=   0.000002800582291; 
	    as[0][19]=   0.000000998029226; bc[0][19]=   0.000001030947803; 
	    as[0][21]=   0.000000367084936; bc[0][21]=   0.000000313752415; 
	    as[0][23]=   0.000000113175783; bc[0][23]=   0.000000116311188; 
	    as[0][25]=   0.000000041990032; bc[0][25]=   0.000000037009179; 
	    as[0][27]=   0.000000013470002; bc[0][27]=   0.000000013791797; 
	    as[0][29]=   0.000000005027753; bc[0][29]=   0.000000004530438; 
	    as[0][31]=   0.000000001660476; bc[0][31]=   0.000000001695357; 
	    as[0][33]=   0.000000000622584; bc[0][33]=   0.000000000570369; 
	    as[0][35]=   0.000000000210229; bc[0][35]=   0.000000000214186; 
	    as[0][37]=   0.000000000079100; bc[0][37]=   0.000000000073418; 
	    as[0][39]=   0.000000000027175; bc[0][39]=   0.000000000027661; 
	    as[1][1]=   0.000953088124940; bc[1][1]=   1.200607556358487; 
	    as[1][3]=  -0.037987408256145; bc[1][3]=  -0.033091761653644; 
	    as[1][5]=   0.004137844595024; bc[1][5]=   0.008348586912074; 
	    as[1][7]=  -0.001437381205791; bc[1][7]=  -0.001333022505006; 
	    as[1][9]=   0.000295358255962; bc[1][9]=   0.000441308286451; 
	    as[1][11]=  -0.000105810296229; bc[1][11]=  -0.000100433041677; 
	    as[1][13]=   0.000026992872406; bc[1][13]=   0.000035475956287; 
	    as[1][15]=  -0.000009835316246; bc[1][15]=  -0.000009449870250; 
	    as[1][17]=   0.000002800582291; bc[1][17]=   0.000003427491382; 
	    as[1][19]=  -0.000001030947803; bc[1][19]=  -0.000000998029226; 
	    as[1][21]=   0.000000313752415; bc[1][21]=   0.000000367084936; 
	    as[1][23]=  -0.000000116311188; bc[1][23]=  -0.000000113175783; 
	    as[1][25]=   0.000000037009179; bc[1][25]=   0.000000041990032; 
	    as[1][27]=  -0.000000013791797; bc[1][27]=  -0.000000013470002; 
	    as[1][29]=   0.000000004530438; bc[1][29]=   0.000000005027753; 
	    as[1][31]=  -0.000000001695357; bc[1][31]=  -0.000000001660476; 
	    as[1][33]=   0.000000000570369; bc[1][33]=   0.000000000622584; 
	    as[1][35]=  -0.000000000214186; bc[1][35]=  -0.000000000210229; 
	    as[1][37]=   0.000000000073418; bc[1][37]=   0.000000000079100; 
	    as[1][39]=  -0.000000000027661; bc[1][39]=  -0.000000000027175; 
	    as[2][1]=  -0.924277418833618; bc[2][1]=  -0.924277418833618; 
	    as[2][3]=   0.003765882001923; bc[2][3]=  -0.003765882001923; 
	    as[2][5]=  -0.009604947313152; bc[2][5]=  -0.009604947313152; 
	    as[2][7]=   0.000080275923681; bc[2][7]=  -0.000080275923681; 
	    as[2][9]=  -0.000566666571087; bc[2][9]=  -0.000566666571087; 
	    as[2][11]=   0.000004136349656; bc[2][11]=  -0.000004136349656; 
	    as[2][13]=  -0.000048052945148; bc[2][13]=  -0.000048052945148; 
	    as[2][15]=   0.000000296496920; bc[2][15]=  -0.000000296496920; 
	    as[2][17]=  -0.000004790825902; bc[2][17]=  -0.000004790825902; 
	    as[2][19]=   0.000000025321982; bc[2][19]=  -0.000000025321982; 
	    as[2][21]=  -0.000000523721040; bc[2][21]=  -0.000000523721040; 
	    as[2][23]=   0.000000002411850; bc[2][23]=  -0.000000002411850; 
	    as[2][25]=  -0.000000060768623; bc[2][25]=  -0.000000060768623; 
	    as[2][27]=   0.000000000247534; bc[2][27]=  -0.000000000247534; 
	    as[2][29]=  -0.000000007352455; bc[2][29]=  -0.000000007352455; 
	    as[2][31]=   0.000000000026831; bc[2][31]=  -0.000000000026831; 
	    as[2][33]=  -0.000000000917657; bc[2][33]=  -0.000000000917657; 
	    as[2][35]=   0.000000000003044; bc[2][35]=  -0.000000000003044; 
	    as[2][37]=  -0.000000000117322; bc[2][37]=  -0.000000000117322; 
	    as[2][39]=   0.000000000000373; bc[2][39]=  -0.000000000000373; 
	    x[0]=   0.000000000000000; y[0]=   0.044958038428844; 
	    x[1]=  -0.000000000000000; y[1]=   1.174900976165617; 
	    x[2]=  -0.000000000000000; y[2]=  -0.938353088149799; 
	    /*
	    vx[0]=   1.356512982939901; vy[0]=  -0.004898058577053; 
	    vx[1]=  -0.100606071800509; vy[1]=  -0.009375247643666; 
	    vx[2]=  -0.966082239338206; vy[2]=   0.010979466323631; 
	    */
	    vx[0]=  -1.356512982939901; vy[0]=   0.004898058577053; 
	    vx[1]=   0.100606071800509; vy[1]=   0.009375247643666; 
	    vx[2]=   0.966082239338206; vy[2]=  -0.010979466323631; 
		} else
		if (name.equals("Ducati3_alluneq")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    m[0]*=  0.709209; rad[n-1] *= Math.pow(0.709209,1.0/3.0);
		    m[1]*=  0.289873; rad[n-1] *= Math.pow(0.289873,1.0/3.0);
		    m[2]*=  1.42189 ; rad[n-1] *= Math.pow(1.421899,1.0/3.0);
		    M=80; N=3; K=40;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
	    as[0][1]=   1.066715861635849; bc[0][1]=   0.517787371379908; 
	    as[0][3]=   0.022671153111922; bc[0][3]=   0.026213314961083; 
	    as[0][5]=   0.003910758444749; bc[0][5]=   0.000611835889285; 
	    as[0][7]=   0.000092578597617; bc[0][7]=   0.000407007347104; 
	    as[0][9]=   0.000128995668928; bc[0][9]=  -0.000039900555274; 
	    as[0][11]=  -0.000013094355661; bc[0][11]=   0.000022537115780; 
	    as[0][13]=   0.000009246437659; bc[0][13]=  -0.000006048073069; 
	    as[0][15]=  -0.000002266738781; bc[0][15]=   0.000002321641590; 
	    as[0][17]=   0.000001018013311; bc[0][17]=  -0.000000823684077; 
	    as[0][19]=  -0.000000340091879; bc[0][19]=   0.000000311560734; 
	    as[0][21]=   0.000000139511392; bc[0][21]=  -0.000000120570609; 
	    as[0][23]=  -0.000000052702095; bc[0][23]=   0.000000047222079; 
	    as[0][25]=   0.000000021466812; bc[0][25]=  -0.000000018990071; 
	    as[0][27]=  -0.000000008583350; bc[0][27]=   0.000000007685006; 
	    as[0][29]=   0.000000003542001; bc[0][29]=  -0.000000003167430; 
	    as[0][31]=  -0.000000001461856; bc[0][31]=   0.000000001314134; 
	    as[0][33]=   0.000000000612994; bc[0][33]=  -0.000000000551902; 
	    as[0][35]=  -0.000000000258357; bc[0][35]=   0.000000000233339; 
	    as[0][37]=   0.000000000109874; bc[0][37]=  -0.000000000099451; 
	    as[0][39]=  -0.000000000047127; bc[0][39]=   0.000000000042753; 
	    as[1][1]=   0.149621824779950; bc[1][1]=   1.472967108604428; 
	    as[1][3]=  -0.060189213117593; bc[1][3]=  -0.055229210902209; 
	    as[1][5]=   0.008823155636588; bc[1][5]=   0.015657308175209; 
	    as[1][7]=  -0.003167692180573; bc[1][7]=  -0.003576806881987; 
	    as[1][9]=   0.000923302427373; bc[1][9]=   0.001251544443551; 
	    as[1][11]=  -0.000348592748112; bc[1][11]=  -0.000408738682415; 
	    as[1][13]=   0.000125406123882; bc[1][13]=   0.000154348196592; 
	    as[1][15]=  -0.000049513699362; bc[1][15]=  -0.000057798689087; 
	    as[1][17]=   0.000019484467591; bc[1][17]=   0.000022949852069; 
	    as[1][19]=  -0.000007975617781; bc[1][19]=  -0.000009185059487; 
	    as[1][21]=   0.000003289393721; bc[1][21]=   0.000003780335622; 
	    as[1][23]=  -0.000001383966417; bc[1][23]=  -0.000001572147180; 
	    as[1][25]=   0.000000587427895; bc[1][25]=   0.000000664240419; 
	    as[1][27]=  -0.000000252333194; bc[1][27]=  -0.000000283259921; 
	    as[1][29]=   0.000000109230591; bc[1][29]=   0.000000122045580; 
	    as[1][31]=  -0.000000047672139; bc[1][31]=  -0.000000052985325; 
	    as[1][33]=   0.000000020938526; bc[1][33]=   0.000000023174452; 
	    as[1][35]=  -0.000000009253271; bc[1][35]=  -0.000000010197803; 
	    as[1][37]=   0.000000004110995; bc[1][37]=   0.000000004512852; 
	    as[1][39]=  -0.000000001835544; bc[1][39]=  -0.000000002006828; 
	    as[2][1]=  -0.562557625614731; bc[2][1]=  -0.558547011669947; 
	    as[2][3]=   0.000962551771670; bc[2][3]=  -0.001815369602247; 
	    as[2][5]=  -0.003749328535137; bc[2][5]=  -0.003497137056480; 
	    as[2][7]=   0.000599603207628; bc[2][7]=   0.000526176826398; 
	    as[2][9]=  -0.000252568675774; bc[2][9]=  -0.000235243720857; 
	    as[2][11]=   0.000077596805090; bc[2][11]=   0.000072086206362; 
	    as[2][13]=  -0.000030177757486; bc[2][13]=  -0.000028449441681; 
	    as[2][15]=   0.000011224676942; bc[2][15]=   0.000010625105662; 
	    as[2][17]=  -0.000004479950765; bc[2][17]=  -0.000004267820115; 
	    as[2][19]=   0.000001795574368; bc[2][19]=   0.000001717106123; 
	    as[2][21]=  -0.000000740175279; bc[2][21]=  -0.000000710537527; 
	    as[2][23]=   0.000000308428037; bc[2][23]=   0.000000296951360; 
	    as[2][25]=  -0.000000130462769; bc[2][25]=  -0.000000125943083; 
	    as[2][27]=   0.000000055722924; bc[2][27]=   0.000000053913476; 
	    as[2][29]=  -0.000000024034894; bc[2][29]=  -0.000000023300893; 
	    as[2][31]=   0.000000010447790; bc[2][31]=   0.000000010146356; 
	    as[2][33]=  -0.000000004574366; bc[2][33]=  -0.000000004449167; 
	    as[2][35]=   0.000000002015275; bc[2][35]=   0.000000001962584; 
	    as[2][37]=  -0.000000000892888; bc[2][37]=  -0.000000000870406; 
	    as[2][39]=   0.000000000397708; bc[2][39]=   0.000000000387796; 
	    x[0]=  -0.000000000000000; y[0]=   0.544997840700911; 
	    x[1]=  -0.000000000000000; y[1]=   1.430774192768532; 
	    x[2]=   0.000000000000000; y[2]=  -0.563517375806289; 
	    /*
	    vx[0]=   1.155930872295865; vy[0]=  -0.006873797478103; 
	    vx[1]=  -0.003359279913766; vy[1]=  -0.011007978258256; 
	    vx[2]=  -0.575868756629479; vy[2]=   0.005672639347606; 
	    */
	    vx[0]=  -1.155930872295865; vy[0]=   0.006873797478103; 
	    vx[1]=   0.003359279913766; vy[1]=   0.011007978258256; 
	    vx[2]=   0.575868756629479; vy[2]=  -0.005672639347606; 
			    } else
		if (name.equals("Ducati3_alluneq2")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    m[0]*=  1.10921 ; rad[n-1] *= Math.pow(1.10921 ,1.0/3.0);
		    m[1]*=  0.689873; rad[n-1] *= Math.pow(0.689873,1.0/3.0);
		    m[2]*=  1.42189 ; rad[n-1] *= Math.pow(1.421899,1.0/3.0);
		    M=80; N=3; K=40;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
	    as[0][1]=   1.103669284058645; bc[0][1]=   0.214881840619989; 
	    as[0][3]=   0.027171550422745; bc[0][3]=   0.032134566739029; 
	    as[0][5]=   0.005849607571281; bc[0][5]=   0.002080729528329; 
	    as[0][7]=   0.000521560159544; bc[0][7]=   0.000782501177522; 
	    as[0][9]=   0.000215606737970; bc[0][9]=   0.000061429454396; 
	    as[0][11]=   0.000016382386769; bc[0][11]=   0.000039791139664; 
	    as[0][13]=   0.000012609118432; bc[0][13]=   0.000001559695586; 
	    as[0][15]=   0.000000259357766; bc[0][15]=   0.000002780366890; 
	    as[0][17]=   0.000000961514889; bc[0][17]=  -0.000000092704349; 
	    as[0][19]=  -0.000000060296763; bc[0][19]=   0.000000242345738; 
	    as[0][21]=   0.000000089917868; bc[0][21]=  -0.000000029864860; 
	    as[0][23]=  -0.000000014037025; bc[0][23]=   0.000000025249604; 
	    as[0][25]=   0.000000009932807; bc[0][25]=  -0.000000005331335; 
	    as[0][27]=  -0.000000002386519; bc[0][27]=   0.000000003033423; 
	    as[0][29]=   0.000000001249314; bc[0][29]=  -0.000000000856113; 
	    as[0][31]=  -0.000000000379721; bc[0][31]=   0.000000000405636; 
	    as[0][33]=   0.000000000172801; bc[0][33]=  -0.000000000134708; 
	    as[0][35]=  -0.000000000059924; bc[0][35]=   0.000000000058606; 
	    as[0][37]=   0.000000000025526; bc[0][37]=  -0.000000000021295; 
	    as[0][39]=  -0.000000000009588; bc[0][39]=   0.000000000008992; 
	    as[1][1]=  -0.033029354905445; bc[1][1]=   1.383083962238464; 
	    as[1][3]=  -0.049185727005706; bc[1][3]=  -0.044573744196884; 
	    as[1][5]=   0.007007090245610; bc[1][5]=   0.012342884271501; 
	    as[1][7]=  -0.002505414145732; bc[1][7]=  -0.002624823493864; 
	    as[1][9]=   0.000699238513784; bc[1][9]=   0.000895521594853; 
	    as[1][11]=  -0.000257154393584; bc[1][11]=  -0.000275515737381; 
	    as[1][13]=   0.000088080483269; bc[1][13]=   0.000100838815504; 
	    as[1][15]=  -0.000033593520807; bc[1][15]=  -0.000035818451018; 
	    as[1][17]=   0.000012589180781; bc[1][17]=   0.000013707530145; 
	    as[1][19]=  -0.000004951295118; bc[1][19]=  -0.000005225170193; 
	    as[1][21]=   0.000001946653598; bc[1][21]=   0.000002064978828; 
	    as[1][23]=  -0.000000784417438; bc[1][23]=  -0.000000819616882; 
	    as[1][25]=   0.000000317587119; bc[1][25]=   0.000000331776791; 
	    as[1][27]=  -0.000000130418186; bc[1][27]=  -0.000000135147625; 
	    as[1][29]=   0.000000053862855; bc[1][29]=   0.000000055717161; 
	    as[1][31]=  -0.000000022450093; bc[1][31]=  -0.000000023111479; 
	    as[1][33]=   0.000000009407533; bc[1][33]=   0.000000009664731; 
	    as[1][35]=  -0.000000003967972; bc[1][35]=  -0.000000004063267; 
	    as[1][37]=   0.000000001681659; bc[1][37]=   0.000000001718355; 
	    as[1][39]=  -0.000000000716363; bc[1][39]=  -0.000000000729957; 
	    as[2][1]=  -0.844940337212697; bc[2][1]=  -0.838672223065054; 
	    as[2][3]=   0.002667556983486; bc[2][3]=  -0.003441700983526; 
	    as[2][5]=  -0.007962939549718; bc[2][5]=  -0.007611680017556; 
	    as[2][7]=   0.000808710094794; bc[2][7]=   0.000663086628883; 
	    as[2][9]=  -0.000507449742209; bc[2][9]=  -0.000482409563764; 
	    as[2][11]=   0.000111986275430; bc[2][11]=   0.000102633804277; 
	    as[2][13]=  -0.000052571149965; bc[2][13]=  -0.000050141647480; 
	    as[2][15]=   0.000016096568482; bc[2][15]=   0.000015209435910; 
	    as[2][17]=  -0.000006858085741; bc[2][17]=  -0.000006578296262; 
	    as[2][19]=   0.000002449304637; bc[2][19]=   0.000002346093730; 
	    as[2][21]=  -0.000001014621027; bc[2][21]=  -0.000000978588265; 
	    as[2][23]=   0.000000391533554; bc[2][23]=   0.000000377964354; 
	    as[2][25]=  -0.000000161835320; bc[2][25]=  -0.000000156812398; 
	    as[2][27]=   0.000000065137955; bc[2][27]=   0.000000063204520; 
	    as[2][29]=  -0.000000027107741; bc[2][29]=  -0.000000026364983; 
	    as[2][31]=   0.000000011188546; bc[2][31]=   0.000000010896785; 
	    as[2][33]=  -0.000000004699145; bc[2][33]=  -0.000000004584046; 
	    as[2][35]=   0.000000001971926; bc[2][35]=   0.000000001925697; 
	    as[2][37]=  -0.000000000835819; bc[2][37]=  -0.000000000817099; 
	    as[2][39]=   0.000000000355044; bc[2][39]=   0.000000000347146; 
	    x[0]=   0.000000000000000; y[0]=   0.249985340910743; 
	    x[1]=  -0.000000000000000; y[1]=   1.348923268587781; 
	    x[2]=  -0.000000000000000; y[2]=  -0.849482170439130; 
	    /*
	    vx[0]=   1.220238110857629; vy[0]=  -0.005279710577315; 
	    vx[1]=  -0.158781455601828; vy[1]=  -0.010582626434576; 
	    vx[2]=  -0.874862589120183; vy[2]=   0.009253144563406; 
	    */
	    vx[0]=  -1.220238110857629; vy[0]=   0.005279710577315; 
	    vx[1]=   0.158781455601828; vy[1]=   0.010582626434576; 
	    vx[2]=   0.874862589120183; vy[2]=  -0.009253144563406; 
					} else
		if (name.equals("FigureEight4")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    M=80; N=n; K=12;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
		    x[0]= -0.053403; y[0]= -0.475557; 
		    x[1]=  1.308483; y[1]=  0.000000; 
		    x[2]= -0.053403; y[2]=  0.475557; 
		    x[3]= -1.201677; y[3]=  0.000000; 
		    vx[0]=  1.090166; vy[0]=  0.649463; 
		    vx[1]=  0.000000; vy[1]= -0.412297; 
		    vx[2]= -1.090166; vy[2]=  0.649463; 
		    vx[3]= -0.000000; vy[3]= -0.886630; 
		    as[0][1]=  1.22554   ; bc[0][1]= -0.36094;
		    as[0][3]= -0.0476103 ; bc[0][3]= -0.0947144;
		    as[0][5]= -0.0132978 ; bc[0][5]= -0.029182;
		    as[0][7]=  0.00904298; bc[0][7]=  0.00767837;
		    as[0][9]=  0.00356273; bc[0][9]=  0.00367382;
		    ac[0][2 ]= -0.0686594; bs[0][ 2]=  0.368149  ;
		    ac[0][6 ]=  0.018014 ; bs[0][ 6]= -0.0189613 ;
		    ac[0][10]= -0.003399 ; bs[0][10]=  0.00356253;
		    for (j=1; j<N; j++) {
			for (k=0; k<K; k++) {
			    as[j][k] = as[0][k]*Math.cos(2*Math.PI*k*j/N)
				     - ac[0][k]*Math.sin(2*Math.PI*k*j/N);
			    ac[j][k] = ac[0][k]*Math.cos(2*Math.PI*k*j/N)
				     + as[0][k]*Math.sin(2*Math.PI*k*j/N);
			    bs[j][k] = bs[0][k]*Math.cos(2*Math.PI*k*j/N)
				     - bc[0][k]*Math.sin(2*Math.PI*k*j/N);
			    bc[j][k] = bc[0][k]*Math.cos(2*Math.PI*k*j/N)
				     + bs[0][k]*Math.sin(2*Math.PI*k*j/N);
			}
		    }
		} else
		if (name.equals("PlateSaucer4")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    M=80; N=4; K=12;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
		    x[0]=   0.190377281458587; y[0]=  -0.579443000126026; 
		    x[1]=  -0.899755923113152; y[1]=  -0.961172971669757; 
		    x[2]=  -0.190377279665413; y[2]=   0.579443000482841; 
		    x[3]=   0.899755922943973; y[3]=   0.961172972059574; 
		    vx[0]=   0.510147268363174; vy[0]=  -0.077596581636764; 
		    vx[1]=  -0.902343124626421; vy[1]=   0.939430273710724; 
		    vx[2]=  -0.510147270525570; vy[2]=   0.077596586252894; 
		    vx[3]=   0.902343124929559; vy[3]=  -0.939430275595325; 
		    as[0][1]=   0.57027; ac[0][1]=   0.11313; 
		    as[0][3]=  -0.00860; ac[0][3]=   0.07435; 
		    as[0][5]=  -0.00235; ac[0][5]=   0.00422; 
		    as[0][7]=  -0.00326; ac[0][7]=  -0.00148; 
		    as[0][9]=  -0.00022; ac[0][9]=  -0.00023; 
		    as[0][11]=   0.00023; ac[0][11]=  -0.00025; 
		    as[1][1]=  -0.97294; ac[1][1]=  -0.89743; 
		    as[1][3]=   0.00226; ac[1][3]=   0.00487; 
		    as[1][5]=   0.01263; ac[1][5]=  -0.00714; 
		    as[1][7]=  -0.00015; ac[1][7]=   0.00002; 
		    as[1][9]=   0.00022; ac[1][9]=   0.00106; 
		    bs[0][1]=   0.11313; bc[0][1]=  -0.57027; 
		    bs[0][3]=  -0.07435; bc[0][3]=  -0.00860; 
		    bs[0][5]=   0.00422; bc[0][5]=   0.00235; 
		    bs[0][7]=   0.00148; bc[0][7]=  -0.00326; 
		    bs[0][9]=  -0.00023; bc[0][9]=   0.00022; 
		    bs[0][11]=   0.00025; bc[0][11]=   0.00023; 
		    bs[1][1]=   0.89743; bc[1][1]=  -0.97294; 
		    bs[1][3]=   0.00487; bc[1][3]=  -0.00226; 
		    bs[1][5]=   0.00714; bc[1][5]=   0.01263; 
		    bs[1][7]=   0.00002; bc[1][7]=   0.00015; 
		    bs[1][9]=  -0.00106; bc[1][9]=   0.00022; 
		    for (j=1; j<N/2; j++) {
			for (k=0; k<K; k++) {
			    as[2*j  ][k] = as[0][k]*Math.cos(4*Math.PI*k*j/N)
				         - ac[0][k]*Math.sin(4*Math.PI*k*j/N);
			    ac[2*j  ][k] = ac[0][k]*Math.cos(4*Math.PI*k*j/N)
				         + as[0][k]*Math.sin(4*Math.PI*k*j/N);
			    bs[2*j  ][k] = bs[0][k]*Math.cos(4*Math.PI*k*j/N)
				         - bc[0][k]*Math.sin(4*Math.PI*k*j/N);
			    bc[2*j  ][k] = bc[0][k]*Math.cos(4*Math.PI*k*j/N)
				         + bs[0][k]*Math.sin(4*Math.PI*k*j/N);
			    as[2*j+1][k] = as[1][k]*Math.cos(4*Math.PI*k*j/N)
				         - ac[1][k]*Math.sin(4*Math.PI*k*j/N);
			    ac[2*j+1][k] = ac[1][k]*Math.cos(4*Math.PI*k*j/N)
				         + as[1][k]*Math.sin(4*Math.PI*k*j/N);
			    bs[2*j+1][k] = bs[1][k]*Math.cos(4*Math.PI*k*j/N)
				         - bc[1][k]*Math.sin(4*Math.PI*k*j/N);
			    bc[2*j+1][k] = bc[1][k]*Math.cos(4*Math.PI*k*j/N)
				         + bs[1][k]*Math.sin(4*Math.PI*k*j/N);
			}
		    }
		} else
		if (name.equals("OrthQuasiEllipse4")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    x[0]=   0.000142836666237; y[0]=   1.387894345301505; 
		    x[1]=   0.002064792013622; y[1]=   0.316077076241605; 
		    x[2]=  -0.000142836666237; y[2]=  -1.387894345301505; 
		    x[3]=  -0.002064792013621; y[3]=  -0.316077076241605; 
		    vx[0]=   0.117159944384558; vy[0]=  -0.001639970899285; 
		    vx[1]=   1.691036997916313; vy[1]=  -0.002402332788673; 
		    vx[2]=  -0.117159944384557; vy[2]=   0.001639970899289; 
		    vx[3]=  -1.691036997916309; vy[3]=   0.002402332788674; 
		    M=240; N=n; K=12;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
		    as[0][1 ]=  1.41373   ; ac[0][ 1]=  -0.00216618;
		    as[0][3 ]=  0.0378243 ; ac[0][ 3]=  -0.000190166;
		    as[0][5 ]=  0.0151218 ; ac[0][ 5]=  -0.000115328;
		    as[0][7 ]=  0.00477286; ac[0][ 7]=  -5.12234e-05;
		    as[0][9 ]=  0.00241482; ac[0][ 9]=  -3.26323e-05;
		    as[0][11]=  0.00107303; ac[0][11]=  -1.7701e-05 ;
		    as[1][1 ]=  0.242059  ; ac[1][ 1]=  -0.00105679 ;
		    as[1][3 ]= -0.0520339 ; ac[1][ 3]=   0.000244734;
		    as[1][5 ]=  0.011894  ; ac[1][ 5]=  -9.35933e-05;
		    as[1][7 ]= -0.00577896; ac[1][ 7]=   6.09092e-05;
		    as[1][9 ]=  0.00237984; ac[1][ 9]=  -3.21944e-05;
		    as[1][11]= -0.00124824; ac[1][11]=   2.04017e-05;

		    bs[0][1 ]=  0.00105679;  bc[0][ 1]=   0.242059;
		    bs[0][3 ]=  0.000244734; bc[0][ 3]=   0.0520339;
		    bs[0][5 ]=  9.35933e-05; bc[0][ 5]=   0.011894;
		    bs[0][7 ]=  6.09092e-05; bc[0][ 7]=   0.00577896;
		    bs[0][9 ]=  3.21944e-05; bc[0][ 9]=   0.00237984;
		    bs[0][11]=  2.04017e-05; bc[0][11]=   0.00124824;
		    bs[1][1 ]=  0.00216618 ; bc[1][ 1]=   1.41373;
		    bs[1][3 ]= -0.000190166; bc[1][ 3]=  -0.0378243;
		    bs[1][5 ]=  0.000115328; bc[1][ 5]=   0.0151218;
		    bs[1][7 ]= -5.12234e-05; bc[1][ 7]=  -0.00477286;
		    bs[1][9 ]=  3.26323e-05; bc[1][ 9]=   0.00241482;
		    bs[1][11]= -1.7701e-05 ; bc[1][11]=  -0.00107303;
		    for (j=1; j<N/2; j++) {
			for (k=0; k<K; k++) {
			    as[2*j  ][k] = as[0][k]*Math.cos(4*Math.PI*k*j/N)
				         - ac[0][k]*Math.sin(4*Math.PI*k*j/N);
			    ac[2*j  ][k] = ac[0][k]*Math.cos(4*Math.PI*k*j/N)
				         + as[0][k]*Math.sin(4*Math.PI*k*j/N);
			    bs[2*j  ][k] = bs[0][k]*Math.cos(4*Math.PI*k*j/N)
				         - bc[0][k]*Math.sin(4*Math.PI*k*j/N);
			    bc[2*j  ][k] = bc[0][k]*Math.cos(4*Math.PI*k*j/N)
				         + bs[0][k]*Math.sin(4*Math.PI*k*j/N);
			    as[2*j+1][k] = as[1][k]*Math.cos(4*Math.PI*k*j/N)
				         - ac[1][k]*Math.sin(4*Math.PI*k*j/N);
			    ac[2*j+1][k] = ac[1][k]*Math.cos(4*Math.PI*k*j/N)
				         + as[1][k]*Math.sin(4*Math.PI*k*j/N);
			    bs[2*j+1][k] = bs[1][k]*Math.cos(4*Math.PI*k*j/N)
				         - bc[1][k]*Math.sin(4*Math.PI*k*j/N);
			    bc[2*j+1][k] = bc[1][k]*Math.cos(4*Math.PI*k*j/N)
				         + bs[1][k]*Math.sin(4*Math.PI*k*j/N);
			}
		    }

		} else
		if (name.equals("Ducati3")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    /*
		    x[0]=  0.000000; y[0]= -0.070947; 
		    x[1]= -0.000000; y[1]=  1.075903; 
		    x[2]= -0.000000; y[2]= -1.004957; 
		    vx[0]=  1.231869; vy[0]=  0.000000; 
		    vx[1]= -0.195090; vy[1]=  0.000000; 
		    vx[2]= -1.036779; vy[2]= -0.000000; 
		    */
		    x[0]=   0.000000000000000; y[0]=  -0.070946548000287; 
		    x[1]=  -0.000000000000000; y[1]=   1.075903540029147; 
		    x[2]=  -0.000000000000000; y[2]=  -1.004956992028862; 
		    /*
		    vx[0]=   1.231870762148251; vy[0]=  -0.000000000000000; 
		    vx[1]=  -0.195089035291636; vy[1]=   0.000000000000000; 
		    vx[2]=  -1.036781726856616; vy[2]=  -0.000000000000000; 
		    */
		    vx[0]=  -1.231870762148251; vy[0]=   0.000000000000000; 
		    vx[1]=   0.195089035291636; vy[1]=  -0.000000000000000; 
		    vx[2]=   1.036781726856616; vy[2]=   0.000000000000000; 
    
		    M=60; N=3; K=12;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
		    as[0][ 1]= 1.097640051076110; bc[0][ 1]=-0.108957450183278;
		    as[0][ 3]= 0.028090182332076; bc[0][ 3]= 0.032506096677853;
		    as[0][ 5]= 0.007235718283341; bc[0][ 5]= 0.003762748006231;
		    as[0][ 7]= 0.001210181377945; bc[0][ 7]= 0.001308690517932;
		    as[0][ 9]= 0.000397710135507; bc[0][ 9]= 0.000289369558798;
		    as[0][11]= 0.000096328571545; bc[0][11]= 0.000101643863860;
    
		    as[1][ 1]=-0.108957450183278; bc[1][ 1]= 1.097640051076109;
		    as[1][ 3]=-0.032506096677853; bc[1][ 3]=-0.028090182332075;
		    as[1][ 5]= 0.003762748006231; bc[1][ 5]= 0.007235718283341;
		    as[1][ 7]=-0.001308690517932; bc[1][ 7]=-0.001210181377945;
		    as[1][ 9]= 0.000289369558798; bc[1][ 9]= 0.000397710135507;
		    as[1][11]=-0.000101643863860; bc[1][11]=-0.000096328571545;
    
		    as[2][ 1]=-0.988682600892831; bc[2][ 1]=-0.988682600892830;
		    as[2][ 3]= 0.004415914345778; bc[2][ 3]=-0.004415914345778;
		    as[2][ 5]=-0.010998466289572; bc[2][ 5]=-0.010998466289572;
		    as[2][ 7]= 0.000098509139987; bc[2][ 7]=-0.000098509139987;
		    as[2][ 9]=-0.000687079694304; bc[2][ 9]=-0.000687079694304;
		    as[2][11]= 0.000005315292314; bc[2][11]=-0.000005315292314;
    
		} else
		if (name.equals("FoldedTriLoop4")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    x[0]=  0.539823; y[0]=  0.265293; 
		    x[1]=  0.094823; y[1]= -0.000000; 
		    x[2]=  0.539823; y[2]= -0.265293; 
		    x[3]= -1.174468; y[3]= -0.000000; 
		    vx[0]=  1.604690; vy[0]= -0.635285; 
		    vx[1]= -0.000000; vy[1]=  1.689055; 
		    vx[2]= -1.604690; vy[2]= -0.635285; 
		    vx[3]= -0.000000; vy[3]= -0.418485; 
		    M=60; N=n; K=19;
		    as = new double[N][K]; 
		    ac = new double[N][K];
		    bc = new double[N][K]; 
		    bs = new double[N][K];
		    as[0][1 ]=  0.821499   ; bc[0][ 1]=  -0.113489;
		    as[0][3 ]=  0.158004   ; bc[0][ 3]=  0.316925;
		    as[0][5 ]=  0.000138446; bc[0][ 5]=  0.0404517;
		    as[0][7 ]=  0.0310197  ; bc[0][ 7]=  0.0157086;
		    as[0][9 ]=  0.00556824 ; bc[0][ 9]=  -0.00231044;
		    as[0][11]=  0.00182091 ; bc[0][11]=  0.00600803;
		    as[0][13]= -0.000615311; bc[0][13]=  0.00141225;
		    as[0][15]=  0.00133346 ; bc[0][15]=  0.000283035;
		    as[0][17]=  0.000389685; bc[0][17]=  -0.00021534;

		    ac[0][2 ]=  0.518392   ; bs[0][ 2]= -0.15166    ;
		    ac[0][6 ]=  0.0124713  ; bs[0][ 6]= -0.0492401  ;
		    ac[0][10]=  0.00854337 ; bs[0][10]= -0.000727563;
		    ac[0][14]= -8.99605e-05; bs[0][14]= -0.00201866 ;
		    ac[0][18]=  0.000519525; bs[0][18]=  0.000119832;

		    for (j=1; j<N; j++) {
			for (k=0; k<K; k++) {
			    as[j][k] = as[0][k]*Math.cos(2*Math.PI*k*j/N)
				     - ac[0][k]*Math.sin(2*Math.PI*k*j/N);
			    ac[j][k] = ac[0][k]*Math.cos(2*Math.PI*k*j/N)
				     + as[0][k]*Math.sin(2*Math.PI*k*j/N);
			    bs[j][k] = bs[0][k]*Math.cos(2*Math.PI*k*j/N)
				     - bc[0][k]*Math.sin(2*Math.PI*k*j/N);
			    bc[j][k] = bc[0][k]*Math.cos(2*Math.PI*k*j/N)
				     + bs[0][k]*Math.sin(2*Math.PI*k*j/N);
			}
		    }
		} else
		if (name.equals("Trefoil4")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    x[0]=  0.528694; y[0]=  0.567126; 
		    x[1]= -0.055073; y[1]=  0.000000; 
		    x[2]=  0.528694; y[2]= -0.567126; 
		    x[3]= -1.002315; y[3]= -0.000000; 
		    vx[0]=  0.175828; vy[0]=  1.023611; 
		    vx[1]=  0.000000; vy[1]= -1.753434; 
		    vx[2]= -0.175828; vy[2]=  1.023611; 
		    vx[3]=  0.000000; vy[3]= -0.293789; 

		    M=240; N=n; K=28;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
		    as[0][1 ]=  0.514114   ; bc[0][ 1]=  0.514114;
		    as[0][5 ]= -0.0527294  ; bc[0][ 5]=  0.0527294;
		    as[0][7 ]= -0.00481555 ; bc[0][ 7]=  -0.00481555;
		    as[0][11]= -0.00637452 ; bc[0][11]=  0.00637452;
		    as[0][13]=  0.000186782; bc[0][13]=  0.000186781;
		    as[0][17]=  0.00128426 ; bc[0][17]=  -0.00128426;
		    as[0][19]=  4.39804e-05; bc[0][19]=  4.39799e-05;
		    as[0][23]=  0.000319484; bc[0][23]=  -0.000319484;
		    as[0][25]= -6.22123e-06; bc[0][25]=  -6.22101e-06;

		    ac[0][2 ]=  0.533093   ; bs[0][ 2]= 0.533093   ;
		    ac[0][10]= -0.000716251; bs[0][10]= 0.000716249;
		    ac[0][14]= -0.00391518 ; bs[0][14]=-0.00391518 ;
		    ac[0][22]=  1.71084e-05; bs[0][22]=-1.71079e-05;
		    ac[0][26]=  0.000234809; bs[0][26]= 0.000234809;

		    for (j=1; j<N; j++) {
			for (k=0; k<K; k++) {
			    as[j][k] = as[0][k]*Math.cos(2*Math.PI*k*j/N)
				     - ac[0][k]*Math.sin(2*Math.PI*k*j/N);
			    ac[j][k] = ac[0][k]*Math.cos(2*Math.PI*k*j/N)
				     + as[0][k]*Math.sin(2*Math.PI*k*j/N);
			    bs[j][k] = bs[0][k]*Math.cos(2*Math.PI*k*j/N)
				     - bc[0][k]*Math.sin(2*Math.PI*k*j/N);
			    bc[j][k] = bc[0][k]*Math.cos(2*Math.PI*k*j/N)
				     + bs[0][k]*Math.sin(2*Math.PI*k*j/N);
			}
		    }
		} else
		if (name.equals("Lagrange2")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    M=240; N=2; K=4;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
		    x[0]=    0.62996820478109; y[0]=   -0.00000000001876; 
		    x[1]=   -0.62996820478109; y[1]=    0.00000000001876; 
		    vx[0]=   -0.00659677045079; vy[0]=    0.62992214996201; 
		    vx[1]=    0.00659677045079; vy[1]=   -0.62992214996201; 
		    as[0][1]=   0.00000; ac[0][1]=   0.62997; 
		    as[1][1]=  -0.00000; ac[1][1]=  -0.62997; 
		    bs[0][1]=   0.62997; bc[0][1]=  -0.00000; 
		    bs[1][1]=  -0.62997; bc[1][1]=   0.00000; 
		} else
		if (name.equals("Lagrange3")) {
		    /*
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    x[0]=  0.524558; y[0]= -0.000000;
		    x[1]= -0.262279; y[1]= -0.454280;
		    x[2]= -0.262279; y[2]=  0.454280;
		    vx[0]=  0.000000; vy[0]=  1.049115;
		    vx[1]=  0.908560; vy[1]= -0.524558;
		    vx[2]= -0.908560; vy[2]= -0.524558;
		    as = ac = bs = bc = null;
		    */
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    M=240; N=3; K=4;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
		    x[0]=    0.83269332620584; y[0]=    0.00000000000632; 
		    x[1]=   -0.41634666324539; y[1]=    0.72113357413130; 
		    x[2]=   -0.41634666323445; y[2]=   -0.72113357413762; 
		    vx[0]=   -0.00871962538236; vy[0]=    0.83263245084674; 
		    vx[1]=   -0.71672104183439; vy[1]=   -0.42386764265378; 
		    vx[2]=    0.72544066721962; vy[2]=   -0.40876480846693; 
		    as[0][1]=  -0.00000; ac[0][1]=   0.83269; 
		    as[1][1]=  -0.72113; ac[1][1]=  -0.41635; 
		    as[2][1]=   0.72113; ac[2][1]=  -0.41635; 
		    bs[0][1]=   0.83269; bc[0][1]=   0.00000; 
		    bs[1][1]=  -0.41635; bc[1][1]=   0.72113; 
		    bs[2][1]=  -0.41635; bc[2][1]=  -0.72113; 
		} else
		if (name.equals("Lagrange20")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    as = ac = bs = bc = null;
		    x[0]= -0.0000000000000000; y[0]=  2.1497222072178794;
		    x[1]=  0.6643006952155469; y[1]=  2.0445073133989644;
		    x[2]=  1.2635750099282934; y[2]=  1.7391617988244870;
		    x[3]=  1.7391617988244867; y[3]=  1.2635750099282934;
		    x[4]=  2.0445073133989644; y[4]=  0.6643006952155470;
		    x[5]=  2.1497222072178794; y[5]=  0.0000000000000001;
		    x[6]=  2.0445073133989649; y[6]= -0.6643006952155468;
		    x[7]=  1.7391617988244867; y[7]= -1.2635750099282932;
		    x[8]=  1.2635750099282936; y[8]= -1.7391617988244867;
		    x[9]=  0.6643006952155471; y[9]= -2.0445073133989644;
		    x[10]=  0.0000000000000002; y[10]= -2.1497222072178794;
		    x[11]= -0.6643006952155459; y[11]= -2.0445073133989649;
		    x[12]= -1.2635750099282932; y[12]= -1.7391617988244867;
		    x[13]= -1.7391617988244865; y[13]= -1.2635750099282936;
		    x[14]= -2.0445073133989644; y[14]= -0.6643006952155474;
		    x[15]= -2.1497222072178794; y[15]= -0.0000000000000004;
		    x[16]= -2.0445073133989649; y[16]=  0.6643006952155466;
		    x[17]= -1.7391617988244870; y[17]=  1.2635750099282930;
		    x[18]= -1.2635750099282939; y[18]=  1.7391617988244863;
		    x[19]= -0.6643006952155475; y[19]=  2.0445073133989644;
		    vx[0]=  2.1497222072178794; vy[0]=  0.0000000000000001;
		    vx[1]=  2.0445073133989644; vy[1]= -0.6643006952155462;
		    vx[2]=  1.7391617988244867; vy[2]= -1.2635750099282936;
		    vx[3]=  1.2635750099282932; vy[3]= -1.7391617988244872;
		    vx[4]=  0.6643006952155467; vy[4]= -2.0445073133989644;
		    vx[5]=  0.0000000000000001; vy[5]= -2.1497222072178794;
		    vx[6]= -0.6643006952155465; vy[6]= -2.0445073133989649;
		    vx[7]= -1.2635750099282930; vy[7]= -1.7391617988244872;
		    vx[8]= -1.7391617988244865; vy[8]= -1.2635750099282939;
		    vx[9]= -2.0445073133989644; vy[9]= -0.6643006952155465;
		    vx[10]= -2.1497222072178794; vy[10]= -0.0000000000000001;
		    vx[11]= -2.0445073133989649; vy[11]=  0.6643006952155457;
		    vx[12]= -1.7391617988244870; vy[12]=  1.2635750099282934;
		    vx[13]= -1.2635750099282934; vy[13]=  1.7391617988244865;
		    vx[14]= -0.6643006952155472; vy[14]=  2.0445073133989644;
		    vx[15]= -0.0000000000000004; vy[15]=  2.1497222072178794;
		    vx[16]=  0.6643006952155465; vy[16]=  2.0445073133989649;
		    vx[17]=  1.2635750099282927; vy[17]=  1.7391617988244870;
		    vx[18]=  1.7391617988244865; vy[18]=  1.2635750099282941;
		    vx[19]=  2.0445073133989644; vy[19]=  0.6643006952155472;
		} else
		if (name.equals("FigureEight3")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    x[0]= -0.000000; y[0]=  0.000000;
		    x[1]=  0.965502; y[1]= -0.241961;
		    x[2]= -0.965502; y[2]=  0.241961;
		    vx[0]=  0.933232; vy[0]=  0.865461;
		    vx[1]= -0.466616; vy[1]= -0.432730;
		    vx[2]= -0.466616; vy[2]= -0.432730;
		    
		    M=60; N=n; K=12;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
		    as[0][ 1]= 1.09586; 
		    bs[0][ 2]= 0.337301;
		    bs[0][ 4]= 0.0557184;
		    as[0][ 5]= -0.0252682; 
		    as[0][ 7]= -0.0058423; 
		    bs[0][ 8]= -0.00299532;
		    bs[0][10]= -0.000805329;
		    as[0][11]= 0.000419168; 
    
		    for (j=1; j<N; j++) {
			for (k=0; k<K; k++) {
			    as[j][k] = as[0][k]*Math.cos(2*Math.PI*k*j/N)
				     - ac[0][k]*Math.sin(2*Math.PI*k*j/N);
			    ac[j][k] = ac[0][k]*Math.cos(2*Math.PI*k*j/N)
				     + as[0][k]*Math.sin(2*Math.PI*k*j/N);
			    bs[j][k] = bs[0][k]*Math.cos(2*Math.PI*k*j/N)
				     - bc[0][k]*Math.sin(2*Math.PI*k*j/N);
			    bc[j][k] = bc[0][k]*Math.cos(2*Math.PI*k*j/N)
				     + bs[0][k]*Math.sin(2*Math.PI*k*j/N);
			}
		    }
		} else
		if (name.equals("FigureEight5")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    x[0]= -0.000000; y[0]= -0.000000;
		    x[1]=  1.360275; y[1]=  0.271516;
		    x[2]=  0.798117; y[2]= -0.484380;
		    x[3]= -0.798117; y[3]=  0.484380;
		    x[4]= -1.360275; y[4]= -0.271516;
		    vx[0]=  1.219124; vy[0]=  1.250650;
		    vx[1]=  0.545471; vy[1]= -0.753552;
		    vx[2]= -1.155034; vy[2]=  0.128226;
		    vx[3]= -1.155034; vy[3]=  0.128226;
		    vx[4]=  0.545471; vy[4]= -0.753552; 

		    M=60; N=n; K=13;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
		    as[0][ 1]=  1.40826; 
		    as[0][ 3]= -0.0373261; 
		    as[0][ 7]= -0.0051869; 
		    as[0][ 9]= -0.00300556; 
		    as[0][11]= -0.00107705; 

		    bs[0][ 2]=  0.499831;
		    bs[0][ 4]=  0.0364537;
		    bs[0][ 6]=  0.0150979;
		    bs[0][ 8]=  0.00318318;
		    bs[0][12]= -0.000417088;
    
		    for (j=1; j<N; j++) {
			for (k=0; k<K; k++) {
			    as[j][k] = as[0][k]*Math.cos(2*Math.PI*k*j/N)
				     - ac[0][k]*Math.sin(2*Math.PI*k*j/N);
			    ac[j][k] = ac[0][k]*Math.cos(2*Math.PI*k*j/N)
				     + as[0][k]*Math.sin(2*Math.PI*k*j/N);
			    bs[j][k] = bs[0][k]*Math.cos(2*Math.PI*k*j/N)
				     - bc[0][k]*Math.sin(2*Math.PI*k*j/N);
			    bc[j][k] = bc[0][k]*Math.cos(2*Math.PI*k*j/N)
				     + bs[0][k]*Math.sin(2*Math.PI*k*j/N);
			}
		    }
		} else
		if (name.equals("FigureEight21")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    ac = as = bc = bs = null;
		    x[0]= -0.000000; y[0]=  0.000000;
		    x[1]=  0.693626; y[1]=  0.566822;
		    x[2]=  1.344507; y[2]=  0.872421;
		    x[3]=  1.998686; y[3]=  0.922288;
		    x[4]=  2.561524; y[4]=  0.666798;
		    x[5]=  2.856775; y[5]=  0.149360;
		    x[6]=  2.754042; y[6]= -0.432179;
		    x[7]=  2.301777; y[7]= -0.833035;
		    x[8]=  1.674675; y[8]= -0.932312;
		    x[9]=  1.017246; y[9]= -0.748207;
		    x[10]=  0.363023; y[10]= -0.333940;
		    x[11]= -0.363023; y[11]=  0.333940;
		    x[12]= -1.017246; y[12]=  0.748207;
		    x[13]= -1.674675; y[13]=  0.932312;
		    x[14]= -2.301777; y[14]=  0.833035;
		    x[15]= -2.754042; y[15]=  0.432179;
		    x[16]= -2.856775; y[16]= -0.149360;
		    x[17]= -2.561524; y[17]= -0.666798;
		    x[18]= -1.998686; y[18]= -0.922288;
		    x[19]= -1.344507; y[19]= -0.872421;
		    x[20]= -0.693626; y[20]= -0.566822;
		    vx[0]=  2.493336; vy[0]=  2.527415;
		    vx[1]=  2.152379; vy[1]=  1.406752;
		    vx[2]=  2.198561; vy[2]=  0.638598;
		    vx[3]=  2.116286; vy[3]= -0.331912;
		    vx[4]=  1.540184; vy[4]= -1.367219;
		    vx[5]=  0.353393; vy[5]= -1.984221;
		    vx[6]= -1.005987; vy[6]= -1.750830;
		    vx[7]= -1.907531; vy[7]= -0.848925;
		    vx[8]= -2.202054; vy[8]=  0.171649;
		    vx[9]= -2.179169; vy[9]=  1.003520;
		    vx[10]= -2.312731; vy[10]=  1.798881;
		    vx[11]= -2.312731; vy[11]=  1.798881;
		    vx[12]= -2.179169; vy[12]=  1.003520;
		    vx[13]= -2.202054; vy[13]=  0.171649;
		    vx[14]= -1.907531; vy[14]= -0.848925;
		    vx[15]= -1.005987; vy[15]= -1.750830;
		    vx[16]=  0.353393; vy[16]= -1.984221;
		    vx[17]=  1.540184; vy[17]= -1.367219;
		    vx[18]=  2.116286; vy[18]= -0.331912;
		    vx[19]=  2.198561; vy[19]=  0.638598;
		    vx[20]=  2.152379; vy[20]=  1.406752;
		} else
		if (name.equals("Trefoil5")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    x[0]=  0.000000; y[0]=  0.000000;
		    x[1]=  0.425476; y[1]=  0.520799;
		    x[2]=  0.720075; y[2]= -0.596513;
		    x[3]= -0.720075; y[3]=  0.596513;
		    x[4]= -0.425476; y[4]= -0.520799;
		    vx[0]=  2.255045; vy[0]=  0.619970;
		    vx[1]= -1.020932; vy[1]= -0.673295;
		    vx[2]= -0.106591; vy[2]=  0.363310;
		    vx[3]= -0.106591; vy[3]=  0.363310;
		    vx[4]= -1.020932; vy[4]= -0.673295; 

		    M=240; N=n; K=50;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
		    as[0][1 ] = 0.643827    ;
		    bs[0][2 ] = 0.675539    ;
		    as[0][3 ] = 0.413228    ;
		    bs[0][4 ] = -0.0761439  ;
		    bs[0][6 ] = 0.0430496   ;
		    as[0][7 ] = 0.0639452   ;
		    bs[0][8 ] = -0.0415704  ;
		    as[0][9 ] =-0.00032588  ;
		    as[0][11] = 0.0159929   ;
		    bs[0][12] = -0.0192728  ;
		    as[0][13] =-0.00149213  ;
		    bs[0][14] = -0.00432179 ;
		    bs[0][16] = -0.00572237 ;
		    as[0][17] =-0.00144205  ;
		    bs[0][18] = -0.0033671  ;
		    as[0][19] =-0.00435731  ;
		    as[0][21] =-0.000879326 ;
		    bs[0][22] = -0.00176042 ;
		    as[0][23] =-0.0039462   ;
		    bs[0][24] = 0.0019967   ;
		    bs[0][26] = -0.000578541;
		    as[0][27] =-0.00230343  ;
		    bs[0][28] = 0.00207375  ;
		    as[0][29] = 0.000479664 ;
		    as[0][31] =-0.000868299 ;
		    bs[0][32] = 0.00135601  ;
		    as[0][33] = 0.000576605 ;
		    bs[0][34] = 0.000128186 ;
		    bs[0][36] = 0.000547986 ;
		    as[0][37] = 0.000428867 ;
		    bs[0][38] = 8.86675e-05 ;
		    as[0][39] = 0.000348361 ;
		    as[0][41] = 0.000200918 ;
		    bs[0][42] = 2.84153e-05 ;
		    as[0][43] = 0.000357392 ;
		    bs[0][44] = -0.000248437;
		    as[0][47] = 0.000227683 ;
		    bs[0][48] = -0.000276428;

		    for (j=1; j<N; j++) {
			for (k=0; k<K; k++) {
			    as[j][k] = as[0][k]*Math.cos(2*Math.PI*k*j/N)
				     - ac[0][k]*Math.sin(2*Math.PI*k*j/N);
			    ac[j][k] = ac[0][k]*Math.cos(2*Math.PI*k*j/N)
				     + as[0][k]*Math.sin(2*Math.PI*k*j/N);
			    bs[j][k] = bs[0][k]*Math.cos(2*Math.PI*k*j/N)
				     - bc[0][k]*Math.sin(2*Math.PI*k*j/N);
			    bc[j][k] = bc[0][k]*Math.cos(2*Math.PI*k*j/N)
				     + bs[0][k]*Math.sin(2*Math.PI*k*j/N);
			}
		    }
		} else
		if (name.equals("Braid4")) {
		    for (j=0; j<n; j++) {m[j] = 1/G; rad[j] = 5*r0;}
		    x[0]= -0.000000; y[0]=  0.157025;
		    x[1]= -1.382856; y[1]= -0.000000;
		    x[2]= -0.000000; y[2]= -0.157025;
		    x[3]=  1.382856; y[3]= -0.000000;
		    vx[0]= -1.871553; vy[0]=  0.000000;
		    vx[1]=  0.000000; vy[1]=  0.584916;
		    vx[2]=  1.871553; vy[2]=  0.000000;
		    vx[3]=  0.000000; vy[3]= -0.584916;

		    M=240; N=n; K=34;
		    as = new double[N][K]; 
		    ac = new double[N][K]; 
		    bs = new double[N][K];
		    bc = new double[N][K];
		    as[0][1 ]= -1.36932    ; bc[0][ 1]= -0.182739;
		    as[0][3 ]=  0.0206285  ; bc[0][ 3]=  0.228366;
		    as[0][5 ]= -0.00846933 ; bc[0][ 5]=  0.0679047;
		    as[0][7 ]= -0.0231746  ; bc[0][ 7]=  0.0142275;
		    as[0][9 ]= -0.0109996  ; bc[0][ 9]=  0.00962889;
		    as[0][11]= -0.00537014 ; bc[0][11]=  0.00700937;
		    as[0][13]= -0.00347894 ; bc[0][13]=  0.00425501;
		    as[0][15]= -0.00245484 ; bc[0][15]=  0.0025754;
		    as[0][17]= -0.00161298 ; bc[0][17]=  0.00173591;
		    as[0][19]= -0.00109424 ; bc[0][19]=  0.00118847;
		    as[0][21]= -0.000757246; bc[0][21]=  0.000820975;
		    as[0][23]= -0.000537598; bc[0][23]=  0.000570515;
		    as[0][25]= -0.000380553; bc[0][25]=  0.000405137;
		    as[0][27]= -0.000273892; bc[0][27]=  0.000289283;
		    as[0][29]= -0.000197625; bc[0][29]=  0.000208928;
		    as[0][31]= -0.000144401; bc[0][31]=  0.000151469;
		    as[0][33]= -0.000105603; bc[0][33]=  0.000110881;

		    for (j=1; j<N; j++) {
			for (k=0; k<K; k++) {
			    as[j][k] = as[0][k]*Math.cos(2*Math.PI*k*j/N)
				     - ac[0][k]*Math.sin(2*Math.PI*k*j/N);
			    ac[j][k] = ac[0][k]*Math.cos(2*Math.PI*k*j/N)
				     + as[0][k]*Math.sin(2*Math.PI*k*j/N);
			    bs[j][k] = bs[0][k]*Math.cos(2*Math.PI*k*j/N)
				     - bc[0][k]*Math.sin(2*Math.PI*k*j/N);
			    bc[j][k] = bc[0][k]*Math.cos(2*Math.PI*k*j/N)
				     + bs[0][k]*Math.sin(2*Math.PI*k*j/N);
			}
		    }
		} else

		if (name.equals("L1-L5 M/m=40")) {
		  double a;
		  r0 = 0.002;
		  if (n>7) {n = 7;}
		  m[0] = 200000; m[1] = 5000; 
		  m[2]=m[3]=m[4]=m[5]=m[6]=0;
		  x[0] = -0.00125; y[0] = 0; 
		  x[1] =  0.05;  y[1] = 0; 
		  //x[0] = -0.5; y[0] = 0; 
		  //x[1] =  1;   y[1] = 0; 
		  v = Math.sqrt(G*m[0]/x[1])/(x[1]-x[0]); 
		  vx[1] = 0; vy[1] = v*x[1]; 

		  as = ac = bs = bc = null;

		  // L4 and L5 form an equilateral triangle with two big masses
		  x[2] = (x[0]+x[1])/2; y[2] = (x[1]-x[2])*Math.sqrt(3); 
		  x[3] = (x[0]+x[1])/2; y[3] =-(x[1]-x[2])*Math.sqrt(3); 
		  th = Math.atan(y[2]/x[2]);
		  a = Math.sqrt(x[2]*x[2]+y[2]*y[2]);
		  vx[2] = -v*a*Math.sin(th); vy[2] = v*a*Math.cos(th); 
		  vx[3] = v*a*Math.sin(th); vy[3] = v*a*Math.cos(th); 

		  // L2, L3, and L4
		  x[4] = 6.1003489037265997e-02; y[4] = 0;
		  x[5] = 4.0373442924202974e-02; y[5] = 0;
		  x[6] =-5.1770792333589014e-02; y[6] = 0;
		  vx[4] = 0; vy[4] = v*x[4]; 
		  vx[5] = 0; vy[5] = v*x[5]; 
		  vx[6] = 0; vy[6] = v*x[6]; 

		  rad[0]=4*r0; rad[1]=2*r0;
		  col[0]=Color.yellow;
		  col[1]=Color.blue;
		  col[2]=col[3]=Color.orange;
		  col[4]=col[5]=col[6]=Color.white;
		  for (j=2; j<n; j++) { rad[j] = r0/2; }
		} else

		if (name.equals("L1-L5 M/m=2")) {
		  double a;
		  r0 = 0.002;
		  if (n>7) {n = 7;}
		  m[0] = 200000; m[1] = 100000; 
		  m[2]=m[3]=m[4]=m[5]=m[6]=0;
		  x[0] = -0.025; y[0] = 0; 
		  x[1] =  0.05;  y[1] = 0; 
		  //x[0] = -0.5; y[0] = 0; 
		  //x[1] =  1;   y[1] = 0; 
		  v = Math.sqrt(G*m[0]/x[1])/(x[1]-x[0]); 
		  vx[1] = 0; vy[1] = v*x[1]; 

		  as = ac = bs = bc = null;

		  // L4 and L5 form an equilateral triangle with two big masses
		  x[2] = (x[0]+x[1])/2; y[2] = (x[1]-x[2])*Math.sqrt(3); 
		  x[3] = (x[0]+x[1])/2; y[3] =-(x[1]-x[2])*Math.sqrt(3); 
		  th = Math.atan(y[2]/x[2]);
		  a = Math.sqrt(x[2]*x[2]+y[2]*y[2]);
		  vx[2] = -v*a*Math.sin(th); vy[2] = v*a*Math.cos(th); 
		  vx[3] = v*a*Math.sin(th); vy[3] = v*a*Math.cos(th); 

		  // L2, L3, and L4
		  x[4] = 9.3678554166024677e-02; y[4] = 0;
		  x[5] = 1.7806367863889522e-02; y[5] = 0;
		  x[6] =-8.5227097049376421e-02; y[6] = 0;
		  vx[4] = 0; vy[4] = v*x[4]; 
		  vx[5] = 0; vy[5] = v*x[5]; 
		  vx[6] = 0; vy[6] = v*x[6]; 

		  rad[0]=4*r0; rad[1]=2*r0;
		  col[0]=Color.yellow;
		  col[1]=Color.blue;
		  col[2]=col[3]=Color.orange;
		  col[4]=col[5]=col[6]=Color.white;
		  for (j=2; j<n; j++) { rad[j] = r0/2; }
		}
		if (name.equals("Belbruno")) {
		  double a;
		  r0 = 0.002;
		  m[0] = 200000; m[1] = 0.000003*m[0]; 
		  m[2]= 0.0300*m[0];
		  m[3]=m[4]=m[5]=m[6]=0;
		  x[0] = -0.05*0.0000003;  y[0] = 0; 
		  x[1] =  0.05;  y[1] = 0; 
		  v = Math.sqrt(G*m[0]/x[1])/(x[1]-x[0]); 
		  vx[1] = 0; vy[1] = v*x[1]; 

		  as = ac = bs = bc = null;

		  // L4 and L5 form an equilateral triangle with two big masses
		  x[2] = (x[0]+x[1])/2; y[2] = (x[1]-x[2])*Math.sqrt(3); 
		  x[3] = (x[0]+x[1])/2; y[3] =-(x[1]-x[2])*Math.sqrt(3); 
		  th = Math.atan(y[2]/x[2]);
		  a = Math.sqrt(x[2]*x[2]+y[2]*y[2]);
		  vx[2] = -v*a*Math.sin(th); vy[2] = v*a*Math.cos(th); 
		  vx[3] = v*a*Math.sin(th); vy[3] = v*a*Math.cos(th); 

		  // L2, L3, and L4
		  x[4] = 9.3678554166024677e-02; y[4] = 0;
		  x[5] = 1.7806367863889522e-02; y[5] = 0;
		  x[6] =-8.5227097049376421e-02; y[6] = 0;
		  vx[4] = 0; vy[4] = v*x[4]; 
		  vx[5] = 0; vy[5] = v*x[5]; 
		  vx[6] = 0; vy[6] = v*x[6]; 

		  rad[0]=4*r0; rad[1]=2*r0;
		  col[0]=Color.yellow;
		  col[1]=Color.blue;
		  col[2]=col[3]=Color.orange;
		  col[3]=Color.white;
		  col[4]=col[5]=col[6]=Color.white;
		  for (j=2; j<n; j++) { rad[j] = r0/2; }
		}
		if (name.equals("BelbrunoVenus")) {
		  /* strange units.  too be fixed 
		   * 0.05 distance = 1 au
		   * 200000 mass = 1 solar mass 
		   * G = 1 */
		  double a;
		  r0 = 0.002;
		  m[0] = 200000;   // sun
		  m[1] = 0.000003001*m[0];  // earth
		  m[2]= 0*m[1];           // lagrange
		  m[3]= 0.000002443*m[0]; // venus
		  // m[3]= 0;
		  m[4]=m[5]=m[6]=0;
		  x[0] = -0.05*0.0000003;  y[0] = 0; 
		  x[1] =  0.05;  y[1] = 0; 
		  x[3] =  0.723*0.05;  y[3] = 0; 
		  v = Math.sqrt(G*m[0]/x[1])/(x[1]-x[0]); 
		  vx[1] = 0; vy[1] = v*x[1]; 

		  as = ac = bs = bc = null;

		  // L4 and L5 form an equilateral triangle with two big masses
		  x[2] = (x[0]+x[1])/2; y[2] = (x[1]-x[2])*Math.sqrt(3); 
		  // x[2] += 0.01*x[1];
		  th = Math.atan(y[2]/x[2]);
		  a = Math.sqrt(x[2]*x[2]+y[2]*y[2]);
		  vx[2] = -v*a*Math.sin(th); vy[2] = v*a*Math.cos(th); 
		  v = Math.sqrt(G*m[0]/x[3])/(x[3]-x[0]); 
		  a = Math.sqrt(x[3]*x[3]+y[3]*y[3]);
		  vx[3] = 0; vy[3] = v*a;

		  // L2, L3, and L4
		  x[4] = 9.3678554166024677e-02; y[4] = 0;
		  x[5] = 1.7806367863889522e-02; y[5] = 0;
		  x[6] =-8.5227097049376421e-02; y[6] = 0;
		  vx[4] = 0; vy[4] = v*x[4]; 
		  vx[5] = 0; vy[5] = v*x[5]; 
		  vx[6] = 0; vy[6] = v*x[6]; 

		  rad[0]=4*r0; rad[1]=2*r0;
		  rad[2]=r0;
		  rad[3]=2*r0;
		  col[0]=Color.yellow;
		  col[1]=Color.blue;
		  col[2]=col[3]=Color.orange;
		  col[3]=Color.white;
		  col[4]=col[5]=col[6]=Color.white;
		  for (j=4; j<n; j++) { rad[j] = r0/2; }
		}
		if (name.equals("Saturn")) {
		  double a;
		  r0 = 0.002;
		  m[0] = 200000;  // Saturn
		  // m[1] = 0.00000003*m[0];  // Ring Dust
		  m[1] = 0.01*m[0]/n;  // Ring Dust
		  for (j=2; j<n; j++) {
		      m[j] = m[1];
		  }
		  x[0] = 0;  y[0] = 0; 
		  r = 0.05;
		  v = Math.sqrt(G*m[0]/r)/r; 
		  vx[1] = 0; vy[1] = 0;

		  as = ac = bs = bc = null;

		  for (j=1; j<n; j++) {
		    th = j/(n-1.0);
		    x[j] = r*Math.cos(th*2*Math.PI);
		    y[j] = r*Math.sin(th*2*Math.PI);
		    vx[j] = -v*r*Math.sin(th*2*Math.PI);
		    vy[j] = v*r*Math.cos(th*2*Math.PI);
		  }

		  rad[0]=4*r0; 
		  col[0]=Color.yellow;
		  for (j=1; j<n; j++) {rad[j]=r0; col[j]=Color.white;}
		}
		if (name.equals("Saturn2")) {
		  double a;
		  r0 = 0.002;
		  m[0] = 200000;  // Saturn
		  // m[1] = 0.00000003*m[0];  // Ring Dust
		  m[1] = 0.01*m[0];  // Ring Dust
		  m[2] = 0.01*m[0]/n;  // Ring Dust
		  for (j=3; j<n; j++) {
		      m[j] = m[2];
		  }
		  r = 0.1;
		  x[0] = -0.01*r;  y[0] = 0; 
		  x[1] =       r;  y[0] = 0; 
		  v = Math.sqrt(G*m[0]/x[1])/(x[1]-x[0]); 
		  vx[0] = 0; vy[0] = -0.01*v*x[1]; 
		  vx[1] = 0; vy[1] = v*x[1]; 

		  as = ac = bs = bc = null;

		  r = 0.05;
		  v = Math.sqrt(G*m[0]/r)/r; 
		  for (j=2; j<n; j++) {
		    th = j/(n-2.0);
		    x[j] = r*Math.cos(th*2*Math.PI);
		    y[j] = r*Math.sin(th*2*Math.PI);
		    vx[j] = -v*r*Math.sin(th*2*Math.PI);
		    vy[j] = v*r*Math.cos(th*2*Math.PI);
		  }

		  rad[0]=4*r0; 
		  col[0]=Color.yellow;
		  rad[1]=r0; 
		  col[1]=Color.red;
		  for (j=2; j<n; j++) {rad[j]=r0; col[j]=Color.white;}
		}
        }
}