ceptreee

function lorenz(h,N,σ,ρ,β,x0,y0,z0,a,b,c)
    # stage
    s = length(b)
    # Number of Equations
    M = 3
    #--------------------
    # Equations
    #--------------------
    f₁(x,y,z) = -σ * x + σ * y
    f₂(x,y,z) = -x * z + ρ * x - y

ceptreee

function lorenz(h,n,a,b,c,x0,y0,z0)
    #--------------------
    # Equation
    #--------------------
    f₁(x,y,z) = -a * x + a * y
    f₂(x,y,z) = -x * z + b * x - y
    f₃(x,y,z) =  x * y - c * z
    #--------------------
    # Initialize
    #--------------------

ceptreee

import numpy as np 
import matplotlib.pyplot as plt 
import sympy as sb
from scipy import integrate

plt.rcParams["font.size"] = 10

# 複素変数
z = sb.Symbol('z')

ceptreee

import numpy as np 
import matplotlib.pyplot as plt 
import sympy as sb

# カラーバー用
from mpl_toolkits.axes_grid1 import make_axes_locatable

# フォントサイズ
plt.rcParams["font.size"] = 10

ceptreee

import numpy as np 
import matplotlib.pyplot as plt 

from matplotlib.colors import LinearSegmentedColormap

clist=[(1,1,1),(1,0,0)]
cm = LinearSegmentedColormap.from_list('wr',clist)
plt.register_cmap(cmap=cm)

clist=[(0,0,1),(1,1,1)]

ceptreee

using Plots; gr(); 
ENV["PLOTS_TEST"] = "true"
using IntervalArithmetic

N    = 27

x0 = @interval(0.11)
x0_sol = 0.11

f(x) = 4*x*(1-x)

ceptreee

import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
from scipy.integrate import odeint

def Lorenz(X,t):
    x, y, z = X    
    dx = -σ * x + σ * y
    dy = ρ * x - y - x * z
    dz = -β * z + x * y

ceptreee

ceptreee

import numpy as np 
import matplotlib.pyplot as plt 

def motion(event):  
    x = event.xdata
    y = event.ydata
    if (x is  None) or (y is  None):
        return
    ln.set_data(x,y)
    print('(x,y) = (%.2f,%.2f)' % (x,y))

ceptreee

import numpy as np 
import matplotlib.pyplot as plt 

def motion(event):  
    x = event.xdata
    y = event.ydata
    if (x is  None) or (y is  None):
        return
    ln.set_data(x,y)
    print('(x,y) = (%.2f,%.2f)' % (x,y))