Хаотическая динамика. Система диф. уравнений Лоренца

public class ChaosCryption {
 
    static float k1, k2, k3, k4, m1, m2, m3, m4, r1, r2, r3, r4;
    static float t, U, S, B;
    static float x = (float)0.01, y = 0, z = 0;
    
    private static float Encode(float point) {
        U = y + point; 
        
        k1=fX(x,y,z, U); 
        m1=fY(x,y,z);
        r1=fZ(x,y,z);
 
        k2=fX(x+k1/2,y+m1/2,z+r1/2, U);
        m2=fY(x+k1/2,y+m1/2,z+r1/2);
        r2=fZ(x+k1/2,y+m1/2,z+r1/2);
 
        k3=fX(x+k2/2,y+m2/2,z+r2/2, U);
        m3=fY(x+k2/2,y+m2/2,z+r2/2);
        r3=fZ(x+k2/2,y+m2/2,z+r2/2);
 
        k4=fX(x+k3,y+m3,z+r3, U);
        m4=fY(x+k3,y+m3,z+r3);
        r4=fZ(x+k3,y+m3,z+r3);
 
        x=x+(k1+k2*2+k3*2+k4)/6;
        y=y+(m1+m2*2+m3*2+m4)/6;
        z=z+(r1+r2*2+r3*2+r4)/6;        
        
        return U + y;
    }
    
    
    private static float Decode(float point) {
        B = point;
        S = B - y;
        
        k1=fX(x,y,z, B);
        m1=fY(x,y,z);
        r1=fZ(x,y,z);
 
        k2=fX(x+k1/2,y+m1/2,z+r1/2, B);
        m2=fY(x+k1/2,y+m1/2,z+r1/2);
        r2=fZ(x+k1/2,y+m1/2,z+r1/2);
 
        k3=fX(x+k2/2,y+m2/2,z+r2/2, B);
        m3=fY(x+k2/2,y+m2/2,z+r2/2);
        r3=fZ(x+k2/2,y+m2/2,z+r2/2);
 
        k4=fX(x+k3,y+m3,z+r3, B);
        m4=fY(x+k3,y+m3,z+r3);
        r4=fZ(x+k3,y+m3,z+r3);
 
        x=x+(k1+k2*2+k3*2+k4)/6 ;
        y=y+(m1+m2*2+m3*2+m4)/6;
        z=z+(r1+r2*2+r3*2+r4)/6;
            
        return S;
    }
    
    private static void Reset() {
        x = (float)0.01;
        y = 0;
        z = 0;
    }
    
    private static float fX(float x, float y, float z, float U) {
        return 10*(U-x);
    }
 
    private static float fY(float x, float y, float z) {
        return 28*x-y-x*z;
    }
 
    private static float fZ(float x, float y, float z) {
        return -(8/3)*z+x*y;
    }
 
    
    
    // интерфейс
    public static float[] Encode(float[] encoding) {
        float[] encoded = new float[encoding.length];
        for(int i=0; i<encoding.length; i++) {
            encoded[i] = Encode(encoding[i]);
        }
        Reset();
        return encoded;
    }
    
    // интерфейс
    public static float[] Decode(float[] decoding) {
        float[] decoded = new float[decoding.length];
        for(int i=0; i<decoding.length; i++) {
            decoded[i] = Decode(decoding[i]);
        }
        Reset();
        return decoded;
    }
}