微分方程式を解き，計算結果を出力したいです．

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

for Vr in np.arange(0.5,0.5,24) :
    a = print(Vr)
    
    def my_eq():

        π = math.pi
        m = 0.0156
        Cn = 6.83879
        D = 0.016
        H = 0.04
        L = 0.3
        ρ = 1.18
        M = (m/(ρ*H*D*L))*(D/H)

        sp.init_printing()
        sp.var('t, C, CL')
        y = sp.Function('y')(t) #xはtの関数(Function)
        C = (4*π*Cn)/Vr
        F = CL/(2*M)
        eq2 = sp.Eq( sp.diff(y, t ,2)+C*sp.diff(y, t, 1),  F )
        sp.dsolve(eq2)
        print(sp.dsolve(eq2))

with pd.ExcelWriter('output.xlsx', mode='w') as writer:
    eq_df = sp.dsolve(eq2)

    eq_df.to_excel(writer, sheet_name=result)

上記のように微分方程式を解き，エクセルに出力したいです．
しかし，どのようにすれば　 sp.dsolve(eq2)　をうまくエクセルに出力できるかわかりません．
まず，sp.dsolve(eq）が毎回書き換わっているような気もします．

また、下記のようにエラーが出てしまします．
解決方法を教えていただきたいです．

IndexError: At least one sheet must be visible