-------------
import numpy as np
import matplotlib.pyplot as plt
import scipy.optimize as opt
import scipy.integrate as quade
boltzmann_const = 1.38e-23
planck_const = 6.62e-34
hbar = planck_const / ( 2 * np.pi )
transition_temp = 9.2
gap_energy_at_zero_kelvin = 3.528 / ( 2 * transition_temp * boltzmann_const )
debye_freq = ( 296 * boltzmann_const ) / hbar
# For subtracting from root_of_integral
a_const = np.log( ( 1.13 * hbar * debye_freq ) / ( boltzmann_const * transition_temp) )
# For simplifying function f.
b_const = ( hbar * debye_freq ) / ( 2 * boltzmann_const)

x_values = np.arange(0.01,0.001,0.0001) 
delta = []
for x in x_values:
    def fun(CE):
        Radius = np.sqrt(CE * CE + x * x) 
        return np.tanh(1477.92 * Radius) / Radius
    integral = quade.quad(fun, 0, 1)
    delta_val = opt.fsolve(lambda CE:integral, 1e-23) - a_const
    delta.append(delta_val)
    delta=np.array(delta)
    plt.plot(delta, x_values)

--------
내가 만들려고하는 코드도 올리겠음

내가 만들고 싶은코드는 저 2번째 그래프 형태임.


그걸 만들기위해서 첫번째 매스매티카 자료를 참조했고


그를 바탕으로 저렇게 만듬


근데 파이썬을 실행하면


저 메세지가 뜨고 plot창이 안뜸


아는 사람 없음?