package rezonans;


import java.util.TimerTask;

public class Zadanie extends TimerTask{
    //#### ZMIENNE  do Bariery
    double A=0.3;
    double PI=3.14;
    int t=0;
    double omega=4*PI;
    double sigma=0.6;


    //####### ZMIENNE DO CZASTECZKI
     int N=1;
     double Xn=0.0;
     double h=0.007;
     //####   FUNKCJE
     public double polozenieCzasteczki(){
       // Metoda iteracyjna obliczająca kolejne położenie cząsteczki w osi OX
       // na podstawie poprzedniego położenie oraz szumu
       double Xnpow3=Xn*Xn*Xn;    //X^3
       double Xnpow2=Xn*Xn;               //X^2
       double temp=Xnpow3-Xnpow2+A*Math.sin(omega*t);// 1 człon równania
       double szum=sigma*Math.sqrt(h)*gauss();  //2 człon równania,szum
       double XnPlus1=Xn-h*(temp) + szum; // Przypisane wartości
       N=N+1;
       Xn=XnPlus1;
       return XnPlus1;
     }

     public double rownanieBariery(double x){
         double temp=A*x*Math.sin(omega*t);
         return 0.25*x*x*x*x  - 0.5*x*x + temp;
     }
    public double gauss(){
        double x1,x2;
        x1=Math.random();
        x2=Math.random();
        double y1,y2;
        double ln =  Math.sqrt(-2*Math.log(x1) );
         y1=ln*Math.cos(2*Math.PI*x2);
        return   y1;
    }


    public void run() {
        t+=5;
        NewClass.aplet.repaint();
    }


}