to resolve the following equations under maxima 5.21.1 with wxmaxima as gui:
eq1:(y[i]-y[1])/(y[1]-y[2])=y[1]*(1+K[L]*rho*y[2])/y[2]/(1+K[L]*rho*y[1]);
eq2:(y[i]-y[1])/(y[2]-y[f])=y[1]*(1+K[L]*rho*y[f])/y[f]/(1+K[L]*rho*y[1]);
ev(solve([eq1,eq2],[y[1],y[2]]),float,y[i]=0.158,y[f]=0.047,rho=1000,K[L]=0.001294);
wxmaxima just shows `maxima is calculating` after almost ten hours.
while maple 13.02 gives the result almost immediately:
eq1:(y[i]-y[1])/(y[1]-y[2])=y[1]*(1+K[L]*rho*y[2])/y[2]/(1+K[L]*rho*y[1])
eq2:(y[i]-y[1])/(y[2]-y[f])=y[1]*(1+K[L]*rho*y[f])/y[f]/(1+K[L]*rho*y[1])
evalf(eval(solve({eq1, eq2}, [y[1], y[2]]), [y[i] = 0.158, y[f] = 0.47, rho
= 1000, K[L] = 0.129]))
[[y[1] = 0.1067788429, y[2] = 0.07118961648]]
So, is it possible adjust maxima parameters to faster the calculate, or it's
just due to lack of certain optimization?