/*
 * @(#)Graph.java	1.9 99/08/04
 * MODIFS PYL fev 2007!!!!!!!!!!
 *
 * This software is designed or intended for use in on-line control of
 * aircraft, air traffic, aircraft navigation or aircraft communications; or in
 * the design, construction, operation or maintenance of any nuclear
 * facility. 
 */
 
//javac guiBrachistochrone.java Kul.java   ;java guiBrachistochrone  
//javac guiBrachistochrone.java Kul.java   ; open index.html  

import java.util.*;
import java.awt.*;
import java.applet.Applet;
import java.awt.event.*;
import javax.swing.*;
import java.io.*;

 
 
 
class Kul extends Thread    
 {
 public double t,Tfinal,Trecord=9999;;
 private double dt=0;
 public double dx,xfinal=4,Ldom=4;
 public double T[],To[],B[]; 
 public double xp,yp,v,vxp,vyp,s;
 public int n,ifinal,flag,visi=0;
  
 boolean pleaseWait = true;
 
private int property;
public int getProperty  (         ){ return property; }
public void setProperty (int value){ property = value; }
public void setdt (double value){ dt = value;System.out.println("------------>>>dt"+dt); }
public double getdt  (         ){/*System.out.println("-------------<<<dt"+dt);*/return dt; }

public void Kul( )
 {    
 }

public void raz0()
{
 xp=0; 
 yp=T[0];
 t=0;
 vxp=0;
 vyp=0; 
 v=0;
 s=0;
 flag=0;
 }

   
 public void avance()
 { int i;
   double g=1;
   double sina=0,cosa,tana=0,ds;
   
    
    i=(int)(xp/dx);
	//System.out.println(i);
	if((i>1)&&(i<n)){
	tana=(T[i+1]-T[i-1])/2/dx; 
	}
	
	if(i<=1){tana=(T[i+1]-T[i])/dx;}//System.out.println(" " +tana  );}

 	
	cosa=Math.sqrt(1/(1+tana*tana));
	sina=tana*cosa;
	
	v = v - dt * g * sina;
	ds = dt * v;
	s = s + ds;
	
	xp = xp + ds * cosa;
	yp = yp + ds * sina;
	
	t=t+dt;
    i=(int)(xp/dx);
    if(yp<T[i])yp=T[i];
    yp=T[i];
    if(i>=ifinal){if(flag==0){Tfinal=t;flag=1; if(Tfinal<Trecord){Trecord=Tfinal;} }}
	
	if(xp>Ldom){raz0(); }
	
	
	
// System.out.println("mouvement  "+xp+"  " +yp +" " +t + "  j " +cosa + " " +sina  );
  } 
  
 public void prems()
 { int i;double x=0;
    //pleaseWait = true;
    t=0;
	n=250;
	dx=Ldom/n;
	ifinal=0;
   System.out.println(" initialize "+ xfinal); 
   T  = new double[n+1];
   B  = new double[n+1];
   To = new double[n+1];
   for(i=0;i<=n;i++)
		{ x=i*dx;
			if(x<xfinal){ifinal++;
			T[i]=1-(x)/(xfinal);}
			 else{T[i]=0;} 
			 To[i]=T[i];
			 }
	 double R=1,the=0;
	 i=0;
	 
	  int month ;
	  month=(int)(2*xfinal);
        switch (month) {
            case 1:  R=1.2028;   break;
            case 2:  R=0.572917; break;
            case 3:  R=0.500663; break;
            case 4:  R=0.5172;   break;
            case 5:  R=0.562182; break;
            case 6:  R=0.619687; break;
            case 7:  R=0.683811; break;
            case 8:  R=0.751873; break;
            default: System.out.println("Invalid month.");break;
        }
	  
	 
	  while(i<n)
	  	{the=the+.01;
	  	 x=R*(the-Math.sin(the));
	     i=(int)(x/dx);
	     B[i]=R*(Math.cos(the)-1)+1;
	     if(x>xfinal){B[i]=0;}
	     // System.out.println("params lisse" + i + " " + x + " " + B[i]  );
	  	  }		  
    raz0();
    }
 
   public void run() {
   int i,it=0;
   
   System.out.println("params lisse" + dt  );
	 		
	while (!interrupted()) {
	    if(dt>0){
		avance();
		try {
		Thread.sleep(1);
	    } catch (InterruptedException e) {}
                synchronized (this) {
                    while (pleaseWait) {
                        try {
                            wait();
                        } catch (Exception e) {}
                    }
                }
    }}
             
}
		
 }