package termico;
import java.util.*;

public class Termico {
    private int N1;
    private int N2;
    private double T1;
    private double T2;
    private double[] energia1;
    private double[] energia2;
    private int t;  //tiempo
//números aleatorios
    private Random rnd=new Random();

    public Termico(int N1, int N2, double T1, double T2) {
        this.N1=N1;
        this.N2=N2;
        this.T1=T1;
        this.T2=T2;
        energia1=new double[N1];
        energia2=new double[N2];
//energía de las partículas de cada subsistema
        for(int i=0; i<N1; i++){
            energia1[i]=T1;
        }
        for(int i=0; i<N2; i++){
            energia2[i]=T2;
        }
        t=0;
    }
    public double getTemperatura1(){return T1;};
    public double getTemperatura2(){return T2;};

    public void evolucion(int dt){
        double Et, Ea, Eb;
        int i, j;
        for(int k=0; k<dt; k++){
            i=(int)(rnd.nextDouble()*N1);
            j=(int)(rnd.nextDouble()*N2);
            Et=energia1[i]+energia2[j];
            Ea=rnd.nextDouble()*Et;
            Eb=Et-Ea;
            T1+=(Ea-energia1[i])/N1;
            T2+=(Eb-energia2[j])/N2;
            energia1[i]=Ea;
            energia2[j]=Eb;
        }
        t+=dt;
    }
}