I need to plot
"Ik" and "rs" are known
with some plotting softwares
My suggestion is that if you can't separate the y variable against the others, then you should use implicit differentiation to get dy/dx and then do a numerical integration based on these values.
Are you familiar with numerical integration?