      d((26),j)=( ((( ( ((-(0.))/((R*T))-( (0.))*(((eta(12))-delta_psi_0
     *))/((R*T))**2)*(Faraday))*(bbb))*EXP(bbb*Faraday*((eta(12))-delta_
     *psi_0)/(R*T)))*(exp(alog((ff(1))/(eta(11))))/((1.0+KNa/(ff(2)))**n
     *nn*(1.0+KCa/(eta(11)))))+((((F_delta(1,j)/((eta(11)))-(0.)*2)/((ff
     *(1))/(eta(11))))*EXP(alog((ff(1))/(eta(11)))))/(((1.0+KNa/(ff(2)))
     ***nnn*(1.0+KCa/(eta(11)))))-((( (-F_delta(2,j)*(KNa)/((ff(2)))**2)
     *)*(nnn)*((1.0+KNa/(ff(2))))**(nnn-1))*((1.0+KCa/(eta(11))))+( (-(0
     *.)*2))*(((1.0+KNa/(ff(2)))**nnn)))*(exp(alog((ff(1))/(eta(11)))))/
     *(((1.0+KNa/(ff(2)))**nnn*(1.0+KCa/(eta(11)))))**2)*(exp(bbb*Farada
     *y*((eta(12))-delta_psi_0)/(R*T))))*(Vmax_NaCa))




      d((26),j)=((((F_delta(1,j)/((eta(11))))/((ff(1))/(eta(11))))*EXP(alog((ff(1))/(eta(11)))))
/((1.0+KNa/(ff(2)))**nnn*(1.0+KCa/(eta(11))))
               -((( (-F_delta(2,j)*(KNa)/((ff(2)))**2)
     )*(nnn)*((1.0+KNa/(ff(2))))**(nnn-1))*((1.0+KCa/(eta(11)))))*(exp(alog((ff(1))/(eta(11)))))/
     (((1.0+KNa/(ff(2)))**nnn*(1.0+KCa/(eta(11)))))**2)

*exp(bbb*Faraday*((eta(12))-delta_psi_0)/(R*T))*(Vmax_NaCa)

