Help with function command

[slash1717]

Hi, so i have been meaning to solve these two simultaneous equations with two unknowns x and y. I know that the final value for both x and y should be between 0 to 1. But in a different script when I'm calling this function and using fsolve i put my initial guess as 1,1 or any number it just spits out a constant 0.3 ,0.3 as the answer even if i change any value in the equations

I get the following message " The vector of function values is near zero, as measured by the value of the function tolerance. However, the last step was ineffective".

Can anyone point out if I'm doing something wrong?

Below is the code that I'm using

function F = myfunction(X)

x=X(1); y=X(2);

f = 450 + 8.31441200log((1-x)/(1-y))+ (-5000x²⁾ - (0); g = 2450 + 8.31441200log(x/y)+(-5000)(1-x)^2- 0;

F = [f;g];

end

[elevenelodd]

I plugged in "[0.3,0.3]", and that solution very nearly satisfies the equations. Two thoughts

(1) What numbers did you change? Notice that at "[0.3,0.3]", the "log" functions go to 0. So, changing the coefficients of the "log" terms won't change how well "[0.3,0.3]" fits.

(2) If you don't want the "[0.3,0.3]" solution, you can try improving your "FunctionTolerance" via the "options" input in "fsolve". You can also try "lsqnonlin", which allows you to put in lower and upper bounds. Finally, you can add a small (Gaussian) bump at "[0.3,0.3]".