Решение нелинейных уравнений для нескольких значений в 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) Как и тогда, как я могу решить для (I), если константы (Rs,Io,Il,Rsh,a) находятся в разных кадрах данных из 8760 строк и 91 столбца?