在python中求解多个值的非线性方程

Question

问题：无法为非线性方程中的多个已知值求解未知

def iv(p):

    Io=2.02*10**-6
    n=1.6
    Isc_ref=4.37
    Il=Isc_ref
    voc_ref=42.93
    Ns=72
    a=(Ns*n*k*Tref)/q
    Rsh=(voc_ref)/(Il-(Io*(np.exp(voc_ref/a)-1)))
    RS=0.01  
    k=1.38*10**-23
    q=1.602*10**-19
    Tref=25+273


    I=p
    v=33.5

    diode=Il-(Io*(np.exp((v+(I*RS))/a)-1))-((v+(I*RS))/Rsh)-I

    return diode

I=opt.fsolve(iv,3)  

（I）是未知的，（v）是已知值，其余是常数。我可以用v的一个值来解决这个问题。
需要解决方案：（1）如何解决（I）多个v值的问题？                      （2）如果常数（Rs，Io，Il，Rsh，a）在8760行和91列的不同数据帧中，我如何解决（I）？