@Benny
Thanks for your help! I did correct the typo of b. I followed a different approach where I copied example from git repo and changed one parameter at a time. It works for my case except 1 operation. I cannot seem to add coefficients to objective function.
#Defining objective functions
fs1 = fill(x -> (x)^b, nvar) #Need to fill with cost coefficients
dfs1 = fill(x -> b/x^(1-b), nvar) #Need to fill with cost coefficients
fs = [x -> fs1[i]*coeffs[i] for i=1:(states*actions)]
dfs = [x -> dfs1[i]*coeffs[i] for i=1:(states*actions)]
If I use fs1 and dfs1 in the sigmoidal programming LinearSP object, it works and gives me solution. But If I use fs and dfs instead it gives me error that var and float64 cannot be multiplied. There is no matching method.
I’m attaching the .jl file. Kindly let me know how to go about with this.
Thank you for all the help! I seriously appreciate it.
Steady_state_EoS.jl (2.6 KB)