How do I find the best parameter estimates with fmincon without knowledge of an initial point?

I am in the final stages of my research and the deadline is closing in. All I need is to get the best estimation for the 15 parameters I have in my objective function. The problem is, I get very large differences in estimates from small changes to the initial point x0. Is there any way for me to find the best estimate or the best initial point? Here is the code that I have developed so far with lots of help from Walter Roberson in another thread. The perms part of the code is used for the constraint that x1 to x5 all have to be distinct. The other constraints can be seen in the function defined at the bottom.

format long g
t = [1;2;3;4;5;6;7;8;9;10;11;12;13;14;15;16;17;18;19;20;21];
Potcap = [0.986702168; 0.97272091; 0.957930961; 0.942424001; 0.926404439; 0.909375829; 0.891635315; 0.873329497; 0.854253377; 0.834442537; 0.813592045; 0.791234153; 0.768357994; 0.744306568; 0.719094568; 0.692488019; 0.66411648; 0.634292184; 0.602498593; 0.569711546; 0.534931653];
x0 = [0.25;0.24;0.23;0.22;0.21;0.24;0.23;0.22;0.21;0.1;0.1;0.1;0.1;0.1;0.3];
lb = zeros(15,1); lb(1:9) = eps(realmin)^(1/5);
ub = [Inf Inf Inf Inf Inf Inf Inf Inf Inf 1 1 1 1 1 Inf];
N = length(x0);
sets = perms(1:5);
Nsets = size(sets,1);
sets(1:5,:)
bestx = nan(1,N);
bestfval = inf;
bestperm = nan(1,5);
fnc = @(x) Fungsi(x,t,Potcap);

options = optimoptions(@fmincon);
options.MaxFunctionEvaluations = 1e4;
options.Algorithm = 'sqp';
disp(options)

for K = 1 : Nsets
   ts = sets(K,:);
   R = zeros(4,N);
   R(sub2ind(size(R), 1:4, ts(1:4))) = 1;
   R(sub2ind(size(R), 1:4, ts(2:5))) = -1;
   b(1:4,1) = -eps(realmin);
   [p, fval, exitflag, output] = fmincon(fnc, x0, R, b, [], [], lb, ub, @mycon, options);
   if exitflag < 0; continue; end
   if fval == bestfval
      bestx = [bestx; p(:).'];
      bestperm = [bestperm; ts];
      fprintf('permutation [%d,%d,%d,%d,%d] equaled fval %g\n', ts, fval);
   elseif fval < bestfval
       bestx = p(:).';
       bestperm = ts;
       bestfval = fval;
       fprintf('permutation [%d,%d,%d,%d,%d] improved fval to %g\n', ts, fval);
   end
end
function f = Fungsi(x, t, Potcap)
   f = sum(((x(10).*exp(-(t.*log(x(1).*x(15) + 1))./x(15)) + x(11).*exp(-(t.*log(x(2).*x(15) + 1))./x(15)) + x(12).*exp(-(t.*log(x(3).*x(15) + 1))./x(15)) + x(13).*exp(-(t.*log(x(4).*x(15) + 1))./x(15)) + x(14).*exp(-(t.*log(x(5).*x(15) + 1))./x(15)) + (x(13).*exp(-(t.*log(x(4).*x(15) + 1))./x(15)).*exp(-(t.*log(x(5).*x(15) + 1))./x(15)).*(x(9).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(9).*exp((t.*log(x(5).*x(15) + 1))./x(15))))./(x(4) - x(5)) + (x(6).*x(10).*exp(-(t.*log(x(1).*x(15) + 1))./x(15)).*exp(-(t.*log(x(2).*x(15) + 1))./x(15)).*(exp((t.*log(x(1).*x(15) + 1))./x(15)) - exp((t.*log(x(2).*x(15) + 1))./x(15))))./(x(1) - x(2)) + (x(7).*x(11).*exp(-(t.*log(x(2).*x(15) + 1))./x(15)).*exp(-(t.*log(x(3).*x(15) + 1))./x(15)).*(exp((t.*log(x(2).*x(15) + 1))./x(15)) - exp((t.*log(x(3).*x(15) + 1))./x(15))))./(x(2) - x(3)) + (x(8).*x(12).*exp(-(t.*log(x(3).*x(15) + 1))./x(15)).*exp(-(t.*log(x(4).*x(15) + 1))./x(15)).*(exp((t.*log(x(3).*x(15) + 1))./x(15)) - exp((t.*log(x(4).*x(15) + 1))./x(15))))./(x(3) - x(4)) + (x(12).*exp(-(t.*log(x(3).*x(15) + 1))./x(15)).*exp(-(t.*log(x(4).*x(15) + 1))./x(15)).*exp(-(t.*log(x(5).*x(15) + 1))./x(15)).*(x(3).*x(8).*x(9).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(3).*x(8).*x(9).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(4).*x(8).*x(9).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) + x(4).*x(8).*x(9).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(5).*x(8).*x(9).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(5).*x(8).*x(9).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15))))./((x(3) - x(4)).*(x(3) - x(5)).*(x(4) - x(5))) + (x(6).*x(7).*x(10).*exp(-(t.*log(x(1).*x(15) + 1))./x(15)).*exp(-(t.*log(x(2).*x(15) + 1))./x(15)).*exp(-(t.*log(x(3).*x(15) + 1))./x(15)).*(x(1).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)) - x(1).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)) - x(2).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)) + x(2).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)) + x(3).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)) - x(3).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15))))./((x(1) - x(2)).*(x(1) - x(3)).*(x(2) - x(3))) + (x(7).*x(8).*x(11).*exp(-(t.*log(x(2).*x(15) + 1))./x(15)).*exp(-(t.*log(x(3).*x(15) + 1))./x(15)).*exp(-(t.*log(x(4).*x(15) + 1))./x(15)).*(x(2).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)) - x(2).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(3).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)) + x(3).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) + x(4).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(4).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15))))./((x(2) - x(3)).*(x(2) - x(4)).*(x(3) - x(4))) - (x(11).*exp(-(t.*log(x(2).*x(15) + 1))./x(15)).*exp(-(t.*log(x(3).*x(15) + 1))./x(15)).*exp(-(t.*log(x(4).*x(15) + 1))./x(15)).*exp(-(t.*log(x(5).*x(15) + 1))./x(15)).*(x(2).*x(3)^2.*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(2)^2.*x(3).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(2).*x(3)^2.*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(2).*x(4)^2.*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) + x(2)^2.*x(3).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(2)^2.*x(4).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) + x(3).*x(4)^2.*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(3)^2.*x(4).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) + x(2).*x(4)^2.*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(2).*x(5)^2.*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(2)^2.*x(4).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(2)^2.*x(5).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(2).*x(5)^2.*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(3).*x(5)^2.*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(2)^2.*x(5).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(3)^2.*x(5).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(3).*x(4)^2.*x(7).*x(8).*x(9).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(3)^2.*x(4).*x(7).*x(8).*x(9).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(3).*x(5)^2.*x(7).*x(8).*x(9).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(4).*x(5)^2.*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(3)^2.*x(5).*x(7).*x(8).*x(9).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(4)^2.*x(5).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(4).*x(5)^2.*x(7).*x(8).*x(9).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(4)^2.*x(5).*x(7).*x(8).*x(9).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15))))./((x(2) - x(3)).*(x(2) - x(4)).*(x(2) - x(5)).*(x(3) - x(4)).*(x(3) - x(5)).*(x(4) - x(5))) - (x(6).*x(7).*x(8).*x(10).*exp(-(t.*log(x(1).*x(15) + 1))./x(15)).*exp(-(t.*log(x(2).*x(15) + 1))./x(15)).*exp(-(t.*log(x(3).*x(15) + 1))./x(15)).*exp(-(t.*log(x(4).*x(15) + 1))./x(15)).*(x(1).*x(2)^2.*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)) - x(1)^2.*x(2).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)) - x(1).*x(2)^2.*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(1).*x(3)^2.*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)) + x(1)^2.*x(2).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) + x(1)^2.*x(3).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)) + x(2).*x(3)^2.*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)) - x(2)^2.*x(3).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)) + x(1).*x(3)^2.*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) + x(1).*x(4)^2.*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(1)^2.*x(3).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(1)^2.*x(4).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(1).*x(4)^2.*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(2).*x(4)^2.*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) + x(1)^2.*x(4).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) + x(2)^2.*x(4).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(2).*x(3)^2.*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) + x(2)^2.*x(3).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) + x(2).*x(4)^2.*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) + x(3).*x(4)^2.*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(2)^2.*x(4).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(3)^2.*x(4).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(3).*x(4)^2.*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) + x(3)^2.*x(4).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15))))./((x(1) - x(2)).*(x(1) - x(3)).*(x(1) - x(4)).*(x(2) - x(3)).*(x(2) - x(4)).*(x(3) - x(4))) - (x(10).*exp(-(t.*log(x(1).*x(15) + 1))./x(15)).*exp(-(t.*log(x(2).*x(15) + 1))./x(15)).*exp(-(t.*log(x(3).*x(15) + 1))./x(15)).*exp(-(t.*log(x(4).*x(15) + 1))./x(15)).*exp(-(t.*log(x(5).*x(15) + 1))./x(15)).*(x(1).*x(2)^2.*x(3)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(1).*x(2)^3.*x(3)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(1)^2.*x(2).*x(3)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) + x(1)^2.*x(2)^3.*x(3).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) + x(1)^3.*x(2).*x(3)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(1)^3.*x(2)^2.*x(3).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(1).*x(2)^2.*x(3)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(1).*x(2)^2.*x(4)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) + x(1).*x(2)^3.*x(3)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1).*x(2)^3.*x(4)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) + x(1)^2.*x(2).*x(3)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1)^2.*x(2).*x(4)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(1)^2.*x(2)^3.*x(3).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(1)^2.*x(2)^3.*x(4).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(1)^3.*x(2).*x(3)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(1)^3.*x(2).*x(4)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) + x(1)^3.*x(2)^2.*x(3).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1)^3.*x(2)^2.*x(4).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) + x(1).*x(3)^2.*x(4)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(1).*x(3)^3.*x(4)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(1)^2.*x(3).*x(4)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) + x(1)^2.*x(3)^3.*x(4).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) + x(1)^3.*x(3).*x(4)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(1)^3.*x(3)^2.*x(4).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) + x(1).*x(2)^2.*x(4)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1).*x(2)^2.*x(5)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(1).*x(2)^3.*x(4)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(1).*x(2)^3.*x(5)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(2).*x(3)^2.*x(4)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) + x(2).*x(3)^3.*x(4)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(1)^2.*x(2).*x(4)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(1)^2.*x(2).*x(5)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1)^2.*x(2)^3.*x(4).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1)^2.*x(2)^3.*x(5).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1)^3.*x(2).*x(4)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1)^3.*x(2).*x(5)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(1)^3.*x(2)^2.*x(4).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(1)^3.*x(2)^2.*x(5).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(2)^2.*x(3).*x(4)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(2)^2.*x(3)^3.*x(4).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(2)^3.*x(3).*x(4)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) + x(2)^3.*x(3)^2.*x(4).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)) - x(1).*x(2)^2.*x(5)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1).*x(2)^3.*x(5)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(1).*x(3)^2.*x(5)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1).*x(3)^3.*x(5)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1)^2.*x(2).*x(5)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1)^2.*x(3).*x(5)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(1)^2.*x(2)^3.*x(5).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(1)^2.*x(3)^3.*x(5).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(1)^3.*x(2).*x(5)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(1)^3.*x(3).*x(5)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1)^3.*x(2)^2.*x(5).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1)^3.*x(3)^2.*x(5).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(1).*x(3)^2.*x(4)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1).*x(3)^3.*x(4)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(2).*x(3)^2.*x(5)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(2).*x(3)^3.*x(5)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1)^2.*x(3).*x(4)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(1)^2.*x(3)^3.*x(4).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(1)^3.*x(3).*x(4)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1)^3.*x(3)^2.*x(4).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(2)^2.*x(3).*x(5)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(2)^2.*x(3)^3.*x(5).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(2)^3.*x(3).*x(5)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(2)^3.*x(3)^2.*x(5).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1).*x(3)^2.*x(5)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(1).*x(3)^3.*x(5)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1).*x(4)^2.*x(5)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(1).*x(4)^3.*x(5)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(1)^2.*x(3).*x(5)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(1)^2.*x(4).*x(5)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1)^2.*x(3)^3.*x(5).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1)^2.*x(4)^3.*x(5).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1)^3.*x(3).*x(5)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1)^3.*x(4).*x(5)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(1)^3.*x(3)^2.*x(5).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(1)^3.*x(4)^2.*x(5).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(1).*x(4)^2.*x(5)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1).*x(4)^3.*x(5)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(2).*x(3)^2.*x(4)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(2).*x(3)^3.*x(4)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(2).*x(4)^2.*x(5)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(2).*x(4)^3.*x(5)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1)^2.*x(4).*x(5)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(1)^2.*x(4)^3.*x(5).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(1)^3.*x(4).*x(5)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(1)^3.*x(4)^2.*x(5).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(2)^2.*x(3).*x(4)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(2)^2.*x(4).*x(5)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(2)^2.*x(3)^3.*x(4).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(2)^2.*x(4)^3.*x(5).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(2)^3.*x(3).*x(4)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(2)^3.*x(4).*x(5)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(2)^3.*x(3)^2.*x(4).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(2)^3.*x(4)^2.*x(5).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(2).*x(3)^2.*x(5)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(2).*x(3)^3.*x(5)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(2)^2.*x(3).*x(5)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(2)^2.*x(3)^3.*x(5).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(2)^3.*x(3).*x(5)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(2)^3.*x(3)^2.*x(5).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(2).*x(4)^2.*x(5)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(2).*x(4)^3.*x(5)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(3).*x(4)^2.*x(5)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(3).*x(4)^3.*x(5)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(2)^2.*x(4).*x(5)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(2)^2.*x(4)^3.*x(5).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(2)^3.*x(4).*x(5)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(2)^3.*x(4)^2.*x(5).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(3)^2.*x(4).*x(5)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(3)^2.*x(4)^3.*x(5).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(3)^3.*x(4).*x(5)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(3)^3.*x(4)^2.*x(5).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(1).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(3).*x(4)^2.*x(5)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(3).*x(4)^3.*x(5)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(3)^2.*x(4).*x(5)^3.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(3)^2.*x(4)^3.*x(5).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) - x(3)^3.*x(4).*x(5)^2.*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15)) + x(3)^3.*x(4)^2.*x(5).*x(6).*x(7).*x(8).*x(9).*exp((t.*log(x(2).*x(15) + 1))./x(15)).*exp((t.*log(x(3).*x(15) + 1))./x(15)).*exp((t.*log(x(4).*x(15) + 1))./x(15)).*exp((t.*log(x(5).*x(15) + 1))./x(15))))./((x(1) - x(2)).*(x(1) - x(3)).*(x(1) - x(4)).*(x(2) - x(3)).*(x(1) - x(5)).*(x(2) - x(4)).*(x(2) - x(5)).*(x(3) - x(4)).*(x(3) - x(5)).*(x(4) - x(5))))-Potcap).^2);
end
function [c,ceq] = mycon(x)
   c(1) = x(6)-x(1);
   c(2) = x(7)-x(2);
   c(3) = x(8)-x(3);
   c(4) = x(9)-x(4);
   ceq = (x(10)+x(11)+x(12)+x(13)+x(14))-1;
end

Any ideas on determining x0 for the problem above? Thanks in advance.

[John D'Errico]

You have apparently only 21 pieces of information. You want to estimate 15 unknown parameters from that data, given a model that has no way to find an intelligent set of starting points, and in fact, a model that surely is going to produce serious numerical problems, just to evaluate it for some if not most sets of parameters.

Is there a valid reason why you think all of those terms, all those parameters are necessary?

Sadly, there is no magic way to find good initial estimates for the parameters of a general nonlinear model. The only thing you might think to do is to try a linear Taylor series approximation, but that has absolutely no value, since a tool like fmincon computes the local gradient of your objective sums of squares. Effectively, it is already using a Taylor series approximation, at every iteration. So to hope that a linear approximation would help is meaningless, since fmincon already does exactly that.

I would suggest fervent prayer, as that may be your best chance. In the end, this is not something I would want to hang my research on, as if I was on the thesis committee, you would get some very strong objections to your approach.

[Rex Mantiri]

Should I try other solvers such as GlobalSearch or MultiStart, in hopes of finding the global minimum maybe?

[John D'Errico]

Multi-start in some form is a good idea, probably your only chance in a million light years of finding a "valid" solution to this problem. But personally, I would strongly, seriously recommend you have a long, heartfelt discussion with your advisor on this modeling choice. If your advisor does not have experience on this sort of problem, then I would strongly suggest you invite an expert in from the statistics or mathematics department to sit in on such a discussion.

Just my opinion.

[Matt J]

It might be worth pointing out that all the expressions exp(-(t.*log(x(j).*x(15) + 1))./x(15)) throughout your objective function can be replaced, after a change of variables, with z(j).^t where z is the new unknown variable. Moreover, since t ranges from 1 to 21, the hgiher order (in t) terms will either be insanely sensitive to changes in z(j) when z(j)>1 or super-insensitive when z(j)<1.

This does not bode well for the optimization, so I think you're going to have to either reduce the range of t or pre-estimate z somehow, converting it from an unknown variable to a known one.

[Matt J]

I would also point out that your c,ceq constraints are all linear expressions, so it would be better to implement them as linear constraints.

[Rex Mantiri]

Oh right, I think I can change the constraints to Aeqx=beq. For the inequality constraints I guess I just add 4 more rows to R right? What would defining z(j) do for the model in this case? Thanks for all the comments btw.

[Matt J]

What would defining z(j) do for the model in this case?

Well, you would have fewer unknowns.

[Rex Mantiri]

Oh yeah that would reduce the total parameters from 15 to 14. The problem is that x(15) is actually one of the main points in the model and it's value is something I have to compute.

[Matt J]

No, if the z(j) are pre-estimated somehow, you would reduce the number of unknowns by half.