function solve_nonleq
clear all;clc
global y
yi=[20 40 60 80 100];
x0=1;
for i=1:4
    y=yi(i);
    [x1(i) fval1(i)]=fsolve(@nonleq1,x0);
    f1(i)=nonleq1(x1(i));
end
x1
fval1


for i=1:4
    y=yi(i);
    [x2(i) fval2(i)]=fsolve(@nonleq2,x0);
    f2(i)=nonleq2(x2(i));
end
x2
fval2

for i=1:4
    y=yi(i);
    [x3(i) fval3(i)]=fsolve(@nonleq3,x0);
    f3(i)=nonleq3(x3(i));
end
x3
fval3




function f=nonleq1(x)
global y
f=143.611*x^5.9285+19.0650*x^0.1139-y;


function f=nonleq2(x)
global y
f=135.3485*x^5.9870+19.3083*x^0.2480-y;

function f=nonleq3(x)
global y
f=140.932*x^5.7776+13.5877*x^0.1920-y;