void make_parent(int t) {
  //if(!automatic) { p++; pmax=p; }
  if(!automatic || !run) {
    if(oldparents) {
      for(int i=0; i<pmax; i++) forrest[i].parent = false;
      fillparent = 0;
      p = 1;
      oldparents = false;
    } else {
      p++;
    }
    pmax = p;
    if(parentlist[fillparent]!=-1) {
      boolean replace = true;
      
      if(replace) {
        forrest[parentlist[fillparent]].elite = false;
        forrest[parentlist[fillparent]].parent = false;
      }
    }
  }
  forrest[t].parent = true;
  parentlist[fillparent] = t;
 
  for(int i=0; i<forrest[t].genotype.length; i++) parents[fillparent].genotype[i] = forrest[t].genotype[i];
  for(int i=0; i<forrest[t].substitute.length; i++) parents[fillparent].substitute[i] = forrest[t].substitute[i];
  parents[fillparent].deg = forrest[t].deg;
  
  fillparent++;
  if(run && automatic) {
    if(fillparent >= pmax) {
      fillparent = 0;
    }
  } if(fillparent >= maxparents) {
    fillparent = 0; p = maxparents; pmax = p;
  }
}

void evolve() {
  println(pmax);
  displfitness=false;
  clear();
  int j;
  //save untouched or manuell-parents
  if(!automatic) {
    oldparents = true;
    for(j=0; j<pmax; j++) {
      for(int i=0; i<parents[j].genotype.length; i++) forrest[j].genotype[i] = parents[j].genotype[i];
      for(int i=0; i<parents[j].substitute.length; i++) forrest[j].substitute[i] = parents[j].substitute[i];
      forrest[j].deg = parents[j].deg;
    }
  } else {
    for(j=0; j<top; j++) {
      for(int i=0; i<parents[j].genotype.length; i++) forrest[j].genotype[i] = parents[j].genotype[i];
      for(int i=0; i<parents[j].substitute.length; i++) forrest[j].substitute[i] = parents[j].substitute[i];
      forrest[j].deg = parents[j].deg;
    }
  }

  for( ; j<t; j++) {
 
    int p1 = int(random(pmax));
    int p2 = p1;
    int rnum;

     // recombination
    if(random(1) < prob_crossover && pmax>1) {
      while(p2==p1) p2 = int(random(pmax));
      rnum = int(random(parents[p1].symbols));
      forrest[j].genotype[0] = parents[p1].genotype[0];
      forrest[j].deg = ((int) random(2)==0) ? parents[p1].deg : parents[p2].deg;
      for(int i=0; i<rnum; i++) forrest[j].substitute[i] = parents[p1].substitute[i];
      for(int i=rnum; i<parents[p2].substitute.length; i++) forrest[j].substitute[i] = parents[p2].substitute[i];
    } else {
      forrest[j].genotype[0] = parents[p1].genotype[0];
      forrest[j].deg = parents[p1].deg;
      for(int i=0; i<parents[p1].substitute.length; i++) forrest[j].substitute[i] = parents[p1].substitute[i];
    }
  
    // mutation
    if(random(1) < prob_mutation) {
      rnum = int(random(parents[p1].symbols+5));
      
      if(rnum<=2) {
        //println("mutate genotype");
        forrest[j].genotype[0] = mutate(forrest[j].genotype[0]);
      } else if(rnum<=4) {
        //println("mutate deg");
        forrest[j].deg += random(-PI/4,PI/4);
        if(forrest[j].deg<PI*0.1) forrest[j].deg=PI*0.1;
        else if(forrest[j].deg>PI) forrest[j].deg=PI;
        //mutate deg
      } else {
        //println("mutate rule "+(rnum-5));
        if((abc[rnum-5]>='A'&&abc[rnum-5]<='Z') || evolve_lowercase) forrest[j].substitute[rnum-5] = mutate(forrest[j].substitute[rnum-5]);
      }
    }

  }
  for(int i=0; i<t; i++) {
    for(j=0; j<abc.length; j++) forrest[i].substitute[j] = clean(forrest[i].substitute[j]);
    forrest[i].genotype[0] = clean(forrest[i].genotype[0]);
    forrest[i].build_gen();
    if(!automatic && i<pmax) {
      forrest[i].parent = true;
      parentlist[i] = i;
    } else if(automatic && i<top) {
      forrest[i].parent = true;
      forrest[i].elite = true;
      parentlist[i] = i;
    }
  }
  measure_fitness();
  //fillparent = 0;
  //if(!automatic) { p=0; pmax=0; }
}

String mutate(String input) {
  String output;
  int rnum, r1 = 0, r2 = 0, v1, v2;
  rnum = int(random(6));
  
  if(rnum == 0) {
  //randomly change 1 letter
    r1 = int(random(input.length()));
    output = input.substring(0,r1) + abc[int(random(abc.length))]+ input.substring(r1+1); 
  } else if (rnum==1) {
  //randomize whole part
    r1 = int(random(input.length())); r2 = int(random(input.length()));
    v1 = min(r1,r2); v2 = max(r1,r2);
    output = input.substring(0,v1);
    for(int k=v1; k<v2; k++) output+=abc[int(random(abc.length))];
    output +=input.substring(v2); 
  } else if (rnum==2) {
  //randomize whole 
    v1 = input.length(); output = "";
    for(int k=0; k<v1; k++) output+=abc[int(random(abc.length))];
  } else if (rnum==3) {
  // minus 1 random letter
    if(input.length()>1) {
      r1 = int(random(input.length()));
      output = input.substring(0,r1) + input.substring(r1+1); 
    } else output=input;
  } else if (rnum==4) {
  // add 1 random letter
    r1 = int(random(input.length()));
    output = input.substring(0,r1) + abc[int(random(abc.length))]+ input.substring(r1); 
  } else if(rnum==5) {
    // add push/pop
    r1 = int(random(input.length()-1));
    r2 = int(random(r1+1,input.length()));
    output = input.substring(0,r1)+'['+input.substring(r1,r2)+']'+input.substring(r2);
  } else output = input;
  return output;
}




void measure_fitness() {

  fitness_checker = true;
}



void clear() {
  for(int i=0; i<t; i++) {
      forrest[i].parent = false;
      forrest[i].elite = false;
  }
}




void rank() {
  println("rank "+pmax);
  oldparents = false; 

  for(int i=0; i<t; i++) {
      forrest[i].parent = false;
      forrest[i].elite = false;
  }
  for(int i=0; i<t; i++) {
    ranking[i] = i;
  }
  int dummy;
  //shuffling
  while(descending()==false) {
    for(int i=t-1; i>0; i--) {
      if(forrest[ranking[i]].fitness > forrest[ranking[i-1]].fitness) {
        dummy = ranking[i];
        ranking[i] = ranking[i-1];
        ranking[i-1] = dummy;
      }
    }
  }

  add_fitness();
  
  // eliminate duplicates of best one
  for(int i=t-1; i>0; i--) {
    for(int j=i-1; j>0; j--) {
      if(forrest[ranking[i]].genotype[g].equals(forrest[ranking[j]].genotype[g])) {
        //println(ranking[i]+"="+ranking[j]);
        forrest[ranking[i]].fitness = f_worst;
        forrest[ranking[i]].startup();
        forrest[ranking[i]].fitness_check = false;
      }
    }
  }
  
  //shuffling again
  while(descending()==false) {
    for(int i=t-1; i>0; i--) {
      if(forrest[ranking[i]].fitness > forrest[ranking[i-1]].fitness) {
        dummy = ranking[i];
        ranking[i] = ranking[i-1];
        ranking[i-1] = dummy;
      }
    }
  }
  
  //clean out parent list
  fillparent = 0;

  for(int i=0; i<pmax; i++) {
    if(!automatic && !run) {
      println(1);
      automatic = true; run=true;
      make_parent(ranking[i]);
      automatic = false;  run=false;
    } else if(!automatic) {
      println(2);
      automatic = true; 
      make_parent(ranking[i]);
      automatic = false;
    } else {
      if(!run) {
        print(3);
        run=true;
        make_parent(ranking[i]);
        run = false;
      } else make_parent(ranking[i]);
      if(i<top) forrest[ranking[i]].elite = true;
    }

  }
  
  displfitness=true;
}

void add_fitness() {
   //compute best and average values
  f_best = 0; f_av = 0; f_worst = 0;
  for(int i=0; i<t; i++) f_av+=forrest[i].fitness;
  f_best = forrest[ranking[0]].fitness;
  f_worst = forrest[ranking[t-1]].fitness;
  f_av /= float(t); 
}

boolean descending() {
  boolean inorder = true;
  for(int i=0; i<t-1; i++) {
    if(forrest[ranking[i]].fitness < forrest[ranking[i+1]].fitness) {
      inorder = false; break;
    }
  }
  return inorder;
}


void randomize() {
  displfitness=false;
  generation=0;
  oldparents = false;
  counter = 0;
  for(int i=0; i<t; i++) {
    forrest[i].startup();
  }
  measure_fitness();
  //p=pp; pmax=p;
}


void switch_mode() {
  automatic = (automatic) ? false : true;
  if(automatic) {
    pmax=pp; fillparent = 0; run = true; bt_run.n = "STOP";
  } else {
    pmax=p; fillparent = 0;
  }
  if(automatic) println("automatic"); else println("manuell");
}


void testplant() {
  
  forrest[0].genotype[0] = "A";
  for(int i=0; i<abc.length; i++) forrest[0].substitute[i]=str(abc[i]);
  forrest[0].substitute[1] = "ARA";
  forrest[0].deg=1.3;
  forrest[0].build_gen();
  
  forrest[1].genotype[0] = "AAAF[R]";
  for(int i=1; i<abc.length; i++) forrest[1].substitute[i]=str(abc[i]);
  forrest[1].substitute[0] = "FF";
  forrest[1].deg=3;
  forrest[1].build_gen();
  
   forrest[2].genotype[0] = "R";
  for(int i=0; i<abc.length; i++) forrest[2].substitute[i]=str(abc[i]);
  forrest[2].substitute[2] = "RAR";
  forrest[2].deg=1.3;
  forrest[2].build_gen();
}


void plant() {
  clear();
  // char[] abc = { 'r','b','j','F','A','J','C','R','K' };
  oldparents = false;
  
  forrest[0].genotype[0] = "FFF[J]F[bFFC]F";
  for(int i=0; i<abc.length; i++) forrest[0].substitute[i]=str(abc[i]);
  forrest[0].substitute[4] = "A[J]F";
  forrest[0].substitute[5] = "F[bF]F[rFA]";
  forrest[0].deg=PI/4;
  forrest[0].build_gen();

  // koch
  /*
  forrest[1].genotype[0] = "F";
  for(int i=0; i<abc.length; i++) forrest[1].substitute[i]=str(abc[i]);
  forrest[1].substitute[3] = "FrrFbbFbbFrrF";
  forrest[1].build_gen();
  */
  forrest[1].genotype[0] = "FFCJK";
  for(int i=0; i<abc.length; i++) forrest[1].substitute[i]=str(abc[i]);
  forrest[1].substitute[5] = "+CJK";
  forrest[1].substitute[8] = "[bF-C]";
  forrest[1].build_gen();
  
  forrest[2].genotype[0] = "F";
  for(int i=0; i<abc.length; i++) forrest[2].substitute[i]=str(abc[i]);
  forrest[2].substitute[3] = "FF[bFJ][rFJ]";
  forrest[2].substitute[5] = "bArrA";
  forrest[2].deg=PI/4;
  forrest[2].build_gen();
  
  forrest[3].genotype[0] = "J";
  for(int i=0; i<abc.length; i++) forrest[3].substitute[i]=str(abc[i]);
  forrest[3].substitute[3] = "FF";
  forrest[3].substitute[5] = "F[rFjC][bFjC]J";
  forrest[3].deg=PI/4;
  forrest[3].build_gen();
  
  forrest[4].genotype[0] = "FFFCJ";
  for(int i=0; i<abc.length; i++) forrest[4].substitute[i]=str(abc[i]);
  forrest[4].substitute[3] = "F";
  forrest[4].substitute[5] = "FRRRRF";
  forrest[4].deg=PI/4;
  forrest[4].build_gen();
 
  forrest[5].genotype[0] = "FJC";
  for(int i=0; i<abc.length; i++) forrest[5].substitute[i]=str(abc[i]);
  forrest[5].substitute[5] = "[bFR][rFA]jJ";
  forrest[5].deg=PI/4;
  forrest[5].build_gen();
  
  forrest[6].genotype[0] = "FF[bJ][rJ]C";
  for(int i=0; i<abc.length; i++) forrest[6].substitute[i]=str(abc[i]);
  forrest[6].substitute[5] = "CbAArrAA";
  forrest[6].substitute[6] = "[RRrJ][AAlJ]";
  forrest[6].deg=PI/4;
  forrest[6].build_gen();
  
   forrest[7].genotype[0] = "F[bFFAA]F[rFJ]";
  for(int i=0; i<abc.length; i++) forrest[7].substitute[i]=str(abc[i]);
  forrest[7].substitute[4] = "F";
  forrest[7].substitute[5] = "bRR[FC]lAA";
  forrest[7].substitute[6] = "[RRrJ][AAl]";
  forrest[7].deg=PI/4;
  forrest[7].build_gen();
  
   forrest[8].genotype[0] = "FFFAAKAAJ";
  for(int i=0; i<abc.length; i++) forrest[8].substitute[i]=str(abc[i]);
  forrest[8].substitute[3] = "F";
  forrest[8].substitute[5] = "[bbAA][rrRR]jC";
  forrest[8].substitute[6] = "[bAAbJ][rRR]";
  forrest[8].substitute[8] = "[bbFRRbK]";
  forrest[8].deg=PI/4;
  forrest[8].build_gen();
  
  forrest[9].genotype[0] = "FFCAARRF[RR][AA]CrrJbK";
  for(int i=0; i<abc.length; i++) forrest[9].substitute[i]=str(abc[i]);
  forrest[9].substitute[5] = "RRKRRbbRR";
  forrest[9].substitute[6] = "[bbAAJ]F";
  forrest[9].substitute[8] = "[rrFJ]F";
  forrest[9].deg=PI/4;
  forrest[9].build_gen();
  
  forrest[10].genotype[0] = "bbFjFrrJ";
  for(int i=0; i<abc.length; i++) forrest[10].substitute[i]=str(abc[i]);
  forrest[10].substitute[5] = "[FRRFKFRRCRRbb]";
  forrest[10].substitute[6] = "C";
  forrest[10].substitute[2] = "Fj";
  forrest[10].substitute[7] = "+CR";
  forrest[10].substitute[8] = "[rrFAAFJFAACj]";
  forrest[10].deg=PI/4;
  forrest[10].build_gen();
  
  forrest[11].genotype[0] = "JC";
  for(int i=0; i<abc.length; i++) forrest[11].substitute[i]=str(abc[i]);
  forrest[11].substitute[5] = "FRRFAAF";
  forrest[11].substitute[7] = "bbAAAArrAAAA";
  forrest[11].substitute[6] = "CFFJ";
  forrest[11].deg=PI/4;
  forrest[11].build_gen();
  
  forrest[12].genotype[0] = "FF[JbbAArrAA][KAArrAA]C";
  for(int i=0; i<abc.length; i++) forrest[12].substitute[i]=str(abc[i]);
  forrest[12].substitute[5] = "JbbAA";
  forrest[12].substitute[4] = "A";
  forrest[12].substitute[7] = "R";
  forrest[12].substitute[8] = "KrrRR";
  forrest[12].substitute[6] = "F[JbbAArrAA][KAArrAA]C";
  forrest[12].deg=PI/4;
  forrest[12].build_gen();
 
  forrest[13].genotype[0] = "FF[J]FFK";
  for(int i=0; i<abc.length; i++) forrest[13].substitute[i]=str(abc[i]);
  forrest[13].substitute[5] = "FRRFK";
  forrest[13].substitute[8] = "AAJ";
  forrest[13].substitute[3] = "F";
  forrest[13].build_gen();
  
  forrest[14].genotype[0] = "F[ArJ]RbK";
  for(int i=0; i<abc.length; i++) forrest[14].substitute[i]=str(abc[i]);
  forrest[14].substitute[5] = "bbAJ";
  forrest[14].substitute[8] = "rrRK";
  forrest[14].deg=PI/4;
  forrest[14].build_gen();
   
  forrest[15].genotype[0] = "FFKFK";
  for(int i=0; i<abc.length; i++) forrest[15].substitute[i]=str(abc[i]);
  forrest[15].substitute[5] = "AAbbbbAA";
  forrest[15].substitute[8] ="AAAAJAAAAJFbbbbK";
  forrest[15].deg=PI/4;
  forrest[15].build_gen();
  
  //char[] abc = { 'r','b','j','F','A','J','C','R','K','+','-' };
  
  forrest[16].genotype[0] = "F";
  for(int i=0; i<abc.length; i++) forrest[16].substitute[i]=str(abc[i]);
  forrest[16].substitute[3] = "F[bFF]f[rFF]F";
  forrest[16].deg=PI/6;
  forrest[16].build_gen();
  
  forrest[17].genotype[0] = "F";
  for(int i=0; i<abc.length; i++) forrest[17].substitute[i]=str(abc[i]);
  forrest[17].substitute[3] = "FFr[rFbFbF]b[bFrFrF]";
  forrest[17].deg=PI/6;
  forrest[17].build_gen();
  
  forrest[18].genotype[0] = "CFJ";
  for(int i=0; i<abc.length; i++) forrest[18].substitute[i]=str(abc[i]);
  forrest[18].substitute[5] = "[C][bjFJ][rjFJ]";
  forrest[18].deg=PI/6;
  forrest[18].build_gen();
  
  forrest[19].genotype[0] = "F+F[bJ][rJ]C";
  for(int i=0; i<abc.length; i++) forrest[19].substitute[i]=str(abc[i]);
  forrest[19].substitute[5] = "Rr-FKACb";
  forrest[19].substitute[6] = "[RRrJ][AAlJ]";
  forrest[19].substitute[7] = "FR";
  forrest[19].build_gen();

  forrest[20].genotype[0] = "FFCAARRF[RR][AA]CrrJbK";
  for(int i=0; i<abc.length; i++) forrest[20].substitute[i]=str(abc[i]);
  forrest[20].substitute[5] = "RRKRRbbRR";
  forrest[20].substitute[6] = "[bbAAJ]F";
  forrest[20].substitute[8] = "[rrFJ]AF";
  forrest[20].build_gen();

  forrest[21].genotype[0] = "FF[J]FFK";
  for(int i=0; i<abc.length; i++) forrest[21].substitute[i]=str(abc[i]);
  forrest[21].substitute[2] = "Fj";
  forrest[21].substitute[5] = "-jCJb-CrrRK+R+C";
  forrest[21].substitute[8] = "[bF-C]";
  forrest[21].build_gen();

  forrest[22].genotype[0] = "FFCAARRF[RR][AA]CrrJbK";
  for(int i=0; i<abc.length; i++) forrest[22].substitute[i]=str(abc[i]);
  forrest[22].substitute[5] = "RRKRRbbRR";
  forrest[22].substitute[8] = "[rrFJ]AF";
  forrest[22].build_gen();

  forrest[23].genotype[0] = "FFCAARRF[RR][AA]C[rrJ]bK";
  for(int i=0; i<abc.length; i++) forrest[23].substitute[i]=str(abc[i]);
  forrest[23].substitute[5] = "RRKRRbbRR";
  forrest[23].substitute[6] = "[bbAAJ]F";
  forrest[23].substitute[7] = "FR";
  forrest[23].substitute[8] = "[rrFJ]F";
  forrest[23].build_gen();

  forrest[24].genotype[0] = "CK-";
  for(int i=0; i<abc.length; i++) forrest[24].substitute[i]=str(abc[i]);
  forrest[24].substitute[4] = "ArJ";
  forrest[24].substitute[6] = "F[JbbAArrAA][KAArrAA]C";
  forrest[24].build_gen();
  
  forrest[25].genotype[0] = "rKFA+b+";
  for(int i=0; i<abc.length; i++) forrest[25].substitute[i]=str(abc[i]);
  forrest[25].substitute[3] = "AF";
  forrest[25].substitute[4] = "[+RFF-]";
  forrest[25].substitute[6] = "bjC";
  forrest[25].substitute[7] = "rF";
  forrest[25].substitute[8] = "+C";
  forrest[25].deg = PI/8;
  forrest[25].build_gen();

  forrest[26].genotype[0] = "ARb";
  for(int i=0; i<abc.length; i++) forrest[26].substitute[i]=str(abc[i]);
  forrest[26].substitute[3] = "FC";
  forrest[26].substitute[4] = "F";
  forrest[26].substitute[5] = "bbAC";
  forrest[26].substitute[6] = "-F[R]";
  forrest[26].substitute[7] = "rrK";
  forrest[26].substitute[8] = "RF";
  forrest[26].deg = PI/4;
  forrest[26].build_gen();
  
  forrest[27].genotype[0] = "Kb+JrJb";
for(int i=0; i<abc.length; i++) forrest[27].substitute[i]=str(abc[i]);
forrest[27].substitute[3] = "Ab+";
forrest[27].substitute[4] = "+A[K]";
forrest[27].substitute[5] = "rJ+";
forrest[27].substitute[6] = "jFRC";
forrest[27].substitute[7] = "-FA";
forrest[27].substitute[8] = "F+Rb";
forrest[27].deg = 1.04954;
forrest[27].build_gen();

forrest[28].genotype[0] = "RrCFACC++JjKrJ";
for(int i=0; i<abc.length; i++) forrest[28].substitute[i]=str(abc[i]);
forrest[28].substitute[4] = "A[J]F";
forrest[28].substitute[5] = "F[bF]F[rFA]";
forrest[28].deg = 0.31415927;
forrest[28].build_gen();


  forrest[29].genotype[0] = "FF[J]FFK";
  for(int i=0; i<abc.length; i++) forrest[29].substitute[i]=str(abc[i]);
  forrest[29].substitute[4] = "[A]";
  forrest[29].substitute[5] = "RrJCCr";
  forrest[29].substitute[6] = "[RRrJ][AAlJ]";
  forrest[29].substitute[7] = "FR";
  forrest[29].build_gen();
  
  forrest[30].genotype[0] = "rKFAJ+b+";
  for(int i=0; i<abc.length; i++) forrest[30].substitute[i]=str(abc[i]);
  forrest[30].substitute[3] = "AF";
  forrest[30].substitute[4] = "[+RFF-]";
  forrest[30].substitute[6] = "Aj";
  forrest[30].substitute[7] = "JRJ";
  forrest[30].substitute[8] = "+[J]";
  forrest[30].build_gen();

  forrest[31].genotype[0] = "F[ArJ]RbK";
  for(int i=0; i<abc.length; i++) forrest[31].substitute[i]=str(abc[i]);
  forrest[31].substitute[4] = "J";
  forrest[31].substitute[5] = "+CJK";
  forrest[31].substitute[8] = "[bF-C]";
  forrest[31].build_gen();

  forrest[32].genotype[0] = "J";
  for(int i=0; i<abc.length; i++) forrest[32].substitute[i]=str(abc[i]); 
  forrest[32].substitute[3] = "FF";
  forrest[32].substitute[4] = "J";
  forrest[32].substitute[5] = "F[rFjC][bFjJ]J";
  forrest[32].substitute[6] = "R";
  forrest[32].substitute[7] = "RCF";
  forrest[32].substitute[8] = "Kj-r";
  forrest[32].build_gen();
  
  forrest[33].genotype[0] = "FF[J]FFK";
  for(int i=0; i<abc.length; i++) forrest[33].substitute[i]=str(abc[i]);
  forrest[33].substitute[5] = "FRRFK";
  forrest[33].substitute[8] = "ACJ";
  forrest[33].deg = 0.39314348;
  forrest[33].build_gen();

  forrest[34].genotype[0] = "j[KJKCj]";
  for(int i=0; i<abc.length; i++) forrest[34].substitute[i]=str(abc[i]);
  forrest[34].substitute[3] = "K+K";
  forrest[34].substitute[4] = "+";
  forrest[34].substitute[5] = "[J]j";
  forrest[34].substitute[6] = "AFC";
  forrest[34].substitute[7] = "[C]";
  forrest[34].substitute[8] = "F[R]";
  forrest[34].deg = 0.5744916;
  forrest[34].build_gen();
  
  forrest[35].genotype[0] = "[R]";
for(int i=0; i<abc.length; i++) forrest[35].substitute[i]=str(abc[i]);
forrest[35].substitute[3] = "FF[bF[j][rF]J]";
forrest[35].substitute[5] = "RbJbb-C-JCR-+";
forrest[35].substitute[6] = "KJr+rjJKCKjr";
forrest[35].substitute[7] = "[K]";
forrest[35].substitute[8] = "[[F]]A";
forrest[35].deg = 1.6615319;
forrest[35].build_gen();


forrest[36].genotype[0] = "KFR";
for(int i=0; i<abc.length; i++) forrest[36].substitute[i]=str(abc[i]);
forrest[36].substitute[3] = "AC";
forrest[36].substitute[5] = "bCbK";
forrest[36].substitute[7] = "C";
forrest[36].substitute[8] = "RKF";
forrest[36].deg = 1.4021204;
forrest[36].build_gen();


forrest[37].genotype[0] = "C+KAAFAJ-bKjrK";
for(int i=0; i<abc.length; i++) forrest[37].substitute[i]=str(abc[i]);
forrest[37].substitute[5] = "C";
forrest[37].substitute[6] = "FC";
forrest[37].substitute[7] = "C";
forrest[37].substitute[8] = "[[bC]]";
forrest[37].deg = 1.9473474;
forrest[37].build_gen();



  forrest[41].genotype[0] = "brA+JFFbCbA";
  for(int i=0; i<abc.length; i++) forrest[41].substitute[i]=str(abc[i]);
  forrest[41].substitute[3] = "C";
  forrest[41].substitute[4] = "C+";
  forrest[41].substitute[5] = "+C";
  forrest[41].substitute[6] = "[-b-F]F";
  forrest[41].substitute[7] = "RR-+";
  forrest[41].substitute[8] = "b";
  forrest[41].deg = 0.5014323;
  forrest[41].build_gen();
  
  forrest[38].genotype[0] = "[KF-A+bK-+j]";
for(int i=0; i<abc.length; i++) forrest[38].substitute[i]=str(abc[i]);
forrest[38].substitute[6] = "FC";
forrest[38].substitute[7] = "F";
forrest[38].substitute[8] = "[C]";
forrest[38].deg = 2.6293314;
forrest[38].build_gen();

forrest[39].genotype[0] = "KrR+bFRj";
for(int i=0; i<abc.length; i++) forrest[39].substitute[i]=str(abc[i]);
forrest[39].substitute[4] = "R++";
forrest[39].substitute[5] = "[Fjj]";
forrest[39].substitute[6] = "AFC";
forrest[39].substitute[7] = "[C]";
forrest[39].substitute[8] = "F[R]";
forrest[39].deg = 1.4766102;
forrest[39].build_gen();

forrest[40].genotype[0] = "jAA[FKbF]";
for(int i=0; i<abc.length; i++) forrest[40].substitute[i]=str(abc[i]);
forrest[40].substitute[3] = "A-A";
forrest[40].substitute[4] = "[R]+";
forrest[40].substitute[5] = "A+-";
forrest[40].substitute[6] = "jRK";
forrest[40].substitute[7] = "A-F";
forrest[40].substitute[8] = "RA";
forrest[40].deg = 1.0405445;
forrest[40].build_gen();



forrest[42].genotype[0] = "FFCCJ";
for(int i=0; i<abc.length; i++) forrest[42].substitute[i]=str(abc[i]);
forrest[42].substitute[3] = "[RF]";
forrest[42].substitute[4] = "[R]+";
forrest[42].substitute[5] = "A-";
forrest[42].substitute[7] = "CF";
forrest[42].substitute[8] = "RA";
forrest[42].deg = 1.0405445;
forrest[42].build_gen();


forrest[44].genotype[0] = "KJF[KF+Cj]";
for(int i=0; i<abc.length; i++) forrest[44].substitute[i]=str(abc[i]);
forrest[44].substitute[3] = "FF+J";
forrest[44].substitute[4] = "RbC";
forrest[44].substitute[5] = "[AA]";
forrest[44].substitute[6] = "rA-";
forrest[44].substitute[7] = "RFJF";
forrest[44].substitute[8] = "[J]";
forrest[44].deg = 1.724913;
forrest[44].build_gen();


  forrest[45].genotype[0] = "KF+R[AKARj]";
  for(int i=0; i<abc.length; i++) forrest[45].substitute[i]=str(abc[i]);
  forrest[45].substitute[4] = "F-";
  forrest[45].substitute[5] = "r-R";
  forrest[45].substitute[6] = "FC";
  forrest[45].substitute[7] = "-CJ";
  forrest[45].substitute[8] = "[C]";
  forrest[45].deg = 2.5565758;
  forrest[45].build_gen();

  forrest[46].genotype[0] = "KF+R+[AKARj]";
  for(int i=0; i<abc.length; i++) forrest[46].substitute[i]=str(abc[i]);
  forrest[46].substitute[3] = "[F]";
forrest[46].substitute[4] = "[A]";
forrest[46].substitute[5] = "C";
forrest[46].substitute[6] = "FC";
forrest[46].substitute[7] = "F";
forrest[46].substitute[8] = "[C]";
forrest[46].deg = 2.3350263;
forrest[46].build_gen();

forrest[48].genotype[0] = "CJr";
for(int i=0; i<abc.length; i++) forrest[48].substitute[i]=str(abc[i]);
forrest[48].substitute[4] = "R++";
forrest[48].substitute[5] = "+b";
forrest[48].substitute[6] = "F[RA]";
forrest[48].substitute[7] = "Cr-R";
forrest[48].substitute[8] = "+-";
forrest[48].deg = 1.4766102;
forrest[48].build_gen();



forrest[49].genotype[0] = "KrrA";
for(int i=0; i<abc.length; i++) forrest[49].substitute[i]=str(abc[i]);
forrest[49].substitute[3] = "-F+A";
forrest[49].substitute[4] = "R[-]";
forrest[49].substitute[5] = "[F]";
forrest[49].substitute[6] = "r";
forrest[49].substitute[7] = "[K]";
forrest[49].substitute[8] = "FA";
forrest[49].deg = 2.4121723;
forrest[49].build_gen();



  
  /*
  forrest[16].genotype[0] = "FF";
  for(int i=0; i<abc.length; i++) forrest[16].substitute[i]=str(abc[i]);
  forrest[16].substitute[3] = "F";
  forrest[16].substitute[4] = "A";
  forrest[16].substitute[5] = "J";
  forrest[16].substitute[6] = "C";
  forrest[16].substitute[7] = "R";
  forrest[16].substitute[8] = "K";
  forrest[16].deg=PI/6;
  forrest[16].build_gen();
  */

  
  //p=15; pmax=15;
  //for(int i=0; i<15; i++) make_parent(i);
 
}