Исправить программу решения системы уравнений на Python

import math
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
 
def dx_dt(x1,y,z, t, x2):
    return -y-z + (Eps*(x1-x2))
def dy_dt(x,y,z, t):
    return x+0.2*y
def dz_dt(x,y,z, t):
    return 0.2+z*(x-5.7)
Eps = 10
def Resslers(n = 200, h = 0.05):
    fig = plt.figure()
    ax = plt.axes(projection='3d')
    ax.set_title('Система Рёсслера 3D')
    def Ressler1(x1 = 0, y1 = 0.8, z1 = 0):
        out = open("out_1.txt", "w")
        a = 0.2
        b = 0.2
        c = 5.7
        g = 0
        o = 0
        t = 0
        x1_list = []
        y1_list = []
        z1_list = []
        x2 = 0
        while t <= n:
            x1 = x1 + dx_dt(x1, y1, z1, t, x2) * h
            y1 = y1 + dy_dt(x1, y1, z1, t) * h
            z1 = z1 + dz_dt(x1, y1, z1, t) * h
            t += h
            x1_list.append(x1)
            y1_list.append(y1)
            z1_list.append(z1)
            out.write("t = %.2f, x = %.19f, y = %.19f,z = %.19f\n" % (t, x1, y1, z1))
        ax.plot3D(x1_list, y1_list, z1_list, 'red')
    Ressler1()
 
 
    #x1_list, y1_list, z1_list = Ressler1()
    def Ressler2(x2 = 0, y2 = 1, z2 = 0):
        out = open("out_2.txt", "w")
        a = 0.2
        b = 0.2
        c = 5.7
        g = 0
        o = 0
        t = 0
        x2_list = []
        y2_list = []
        z2_list = []
        x1 = 0
        while t <= n:
            x2 = x2 + dx_dt(x2, y2, z2, t, x1) * h
            y2 = y2 + dy_dt(x2, y2, z2, t) * h
            z2 = z2 + dz_dt(x2, y2, z2, t) * h
            t += h
            x2_list.append(x2)
            y2_list.append(y2)
            z2_list.append(z2)
            out.write("t = %.2f, x = %.19f, y = %.19f,z = %.19f\n" % (t, x2, y2, z2))
        ax.plot3D(x2_list, y2_list, z2_list)
        plt.xlabel("Ось Х")
        plt.ylabel("Ось Y")
        plt.show()
    Ressler2()
 
Resslers()