解が小さすぎてすべて０になります．

Question

import sympy as sp
import math
import cv2
import numpy as np
import matplotlib.pyplot as plt
import csv
import os
import glob
import pandas as pd

options = {}
options['strings_to_formulas'] = False
options['strings_to_urls'] = False
with pd.ExcelWriter('eq.xlsx', mode='w') as writer:
    results = []
    for i in range(28) :
        print(i)
        U = i*0.25 + 0.25
        print(U)
        π = math.pi
        m = 0.0156
        δ = 0.0356167639226128
        D = 0.016
        H = 0.04
        L = 0.3
        fc = 9.0
        ρ = 1.18
        M = (m/(ρ*H*D*L))*(D/H)
        Vr = U/(fc*H)
        A = H*L
        
        lists = []
        sp.init_printing()
        sp.var('t')
        y = sp.Function('y')(t) #xはtの関数(Function)
        α = (4*π*δ)/Vr
        ω = (2*π)/Vr
        dCL = 0.3 #dCL/dα
        T = Vr #=2π/ωn
        eq2 = sp.Eq( sp.diff(y, t ,2)+(α+((ρ*U*A)/2*M)*dCL)*sp.diff(y, t, 1)+(ω*ω)*y,  0 )
        yre = sp.dsolve(eq2, ics={y.subs(t,0):0, y.diff(t,2).subs(t,0):0})
        
        for j in range(0, 4001):
            print(j)
            lists.append(yre.subs(t, j*0.0005))
        results.append(lists)

    df = pd.DataFrame(results).T
    
    df.to_excel(writer)

上記の運動方程式を解きたいのですが，解が小さすぎるのか表示がすべて0になってしまいます．
どうしたらよいでしょうか？

Answer 1

print(pd.options.display.precision)
6
pd.options.display.precision = 10

などで指定してみるのはどうでしょうか？

https://note.nkmk.me/python-pandas-option-display/

型はfloatになっていますか？
.info()で各要素確認してみてください。
.astype(float)で型変換してみるのもいいかもしれません。