Failure in initial objective function evaluation. FSOLVE cannot continue.

I am having this problem in MATLAB.

B=[0,0.0574612898073843,0.110961034961411,0.157844338275811,0.198326084545702,0.233505068208522,0.264466281579592,0.291951607426117,0.316302322044398,0.337552164906992,0.355641291820062]; E=[0.00821682460973364,0.0230168805223958,0.0353571800836649,0.0452377232935416,0.0526585101520253,0.0576195406591166,0.0601208148148149,0.0601623326191205,0.0577440940720333,0.0528660991735539,0.0454999999999998]; F=[0.0789623507805327,0.0592582389552085,0.0395296398326703,0.0197765534129168,-1.02030405058073e-06,-0.0198030813182328,-0.0396296296296297,-0.0594806652382411,-0.0793561881440668,-0.0992561983471074,-0.119200000000000]; Q0=[0.248000000000000,0.256216824609734,0.279233705132129,0.314590885215794,0.359828608509336,0.412487118661361,0.470106659320478,0.530227474135293,0.590389806754413,0.648133900826446,0.701000000000000]; Q1=[0.448000000000000,0.535179175390266,0.617454294867871,0.692341114784206,0.757355391490664,0.810012881338639,0.847829340679523,0.868320525864708,0.869002193245587,0.847390099173554,0.801000000000000];

X0=ones(1,11); A=fsolve(@(X)Func(X,B,E,F,Q0,Q1),X0,optimoptions('fsolve','Display','iter','MaxIterations',10000,'MaxFunctionEvaluations',50000,'FunctionTolerance',10^-8,'Diagnostic','on'));

tmp = -((Q1-Q0*exp(A))./(E+((EA+F)(0.5+(A./6)+((A.^2)/24))))); differenza=B+tmp

function [Y] = Func(X,B,E,F,Q0,Q1) Y= B-((Q1-Q0*exp(X))./(E+((EX+F)(0.5+(X./6)+((X.^2)/24))))); end

When I run it , it says:

Failure in initial objective function evaluation. FSOLVE cannot continue.

CAN ANYONE HELP ME, PLEASE?

Answer 1

[Star Strider]

The ‘Y’ evaluation is missing a few multiplication(?) operators, and needs to do element-wise operations elsewhere:

Y = B-((Q1-Q0*exp(X))./(E+((EX+F)(0.5+(X./6)+((X.^2)/24)))));
            ↑              ↑    ↑

Rewriting ‘Func’ as an anonymous function and correcting those (and similar errors in ‘tmp’):

Func = @(X,B,E,F,Q0,Q1) B-((Q1-Q0.*exp(X))./(E+((E.*X+F).*(0.5+(X./6)+((X.^2)/24)))));

as well as making your code actually readable:

Func = @(X,B,E,F,Q0,Q1) B-((Q1-Q0.*exp(X))./(E+((E.*X+F).*(0.5+(X./6)+((X.^2)/24)))));
B=[0,0.0574612898073843,0.110961034961411,0.157844338275811,0.198326084545702,0.233505068208522,0.264466281579592,0.291951607426117,0.316302322044398,0.337552164906992,0.355641291820062]; 
E=[0.00821682460973364,0.0230168805223958,0.0353571800836649,0.0452377232935416,0.0526585101520253,0.0576195406591166,0.0601208148148149,0.0601623326191205,0.0577440940720333,0.0528660991735539,0.0454999999999998]; 
F=[0.0789623507805327,0.0592582389552085,0.0395296398326703,0.0197765534129168,-1.02030405058073e-06,-0.0198030813182328,-0.0396296296296297,-0.0594806652382411,-0.0793561881440668,-0.0992561983471074,-0.119200000000000]; 
Q0=[0.248000000000000,0.256216824609734,0.279233705132129,0.314590885215794,0.359828608509336,0.412487118661361,0.470106659320478,0.530227474135293,0.590389806754413,0.648133900826446,0.701000000000000]; 
Q1=[0.448000000000000,0.535179175390266,0.617454294867871,0.692341114784206,0.757355391490664,0.810012881338639,0.847829340679523,0.868320525864708,0.869002193245587,0.847390099173554,0.801000000000000];
X0=ones(1,11);
A=fsolve(@(X)Func(X,B,E,F,Q0,Q1),X0,optimoptions('fsolve','Display','iter','MaxIterations',10000,'MaxFunctionEvaluations',50000,'FunctionTolerance',10^-8,'Diagnostic','on'));
tmp = -((Q1-Q0.*exp(A))./(E+((E.*A+F).*(0.5+(A./6)+((A.^2)/24)))));
differenza=B+tmp

produces:

A =
  0.59136      0.72882       0.7792      0.77019       0.7238       0.6548       0.5718        0.479      0.37742      0.26558      0.13939