Hi Bhavesh,
The likely cause of the error is that the “exp” function expects dimension-less input. From the error message you provided, it seems that your unit of the gas constant (R) does not align with unit of “Ea” which is making the result of “-Ea / (R * Tc)” to have some dimension. Make sure that the units of Ea and R align to make it work. You may consider expressing Ea in “J/mol” and R in “m^2*kg/s^2/K/mol”. This will make sure that result of “-Ea / (R * Tc)” is dimensionless and then you can pass it to the “exp” function as argument.
I hope this solves the problem!
