% Thurber
% http://www.itl.nist.gov/div898/strd/nls/data/LINKS/DATA/Thurber.dat
dim = 7 % dimension of the minimization problem
nb = 37 % number of points

% NIST certified values
target = zeros(dim,1);
target(1) = 1.2881396800E+03;
target(2) = 1.4910792535E+03;
target(3) = 5.8323836877E+02;
target(4) = 7.5416644291E+01;
target(5) = 9.6629502864E-01;
target(6) = 3.9797285797E-01;
target(7) = 4.9727297349E-02;

% start 1
ig1 = [1000;1000;400;40;0.7;0.3;0.03]

% start 2
ig2 = [1300;1500;500;75;1;0.4;0.05]

% points to fit
x = zeros(nb,1)
y = zeros(nb,1)
x(1)=-3.067;	y(1)=80.574
x(2)=-2.981;	y(2)=84.248
x(3)=-2.921;	y(3)=87.264
x(4)=-2.912;	y(4)=87.195
x(5)=-2.84;	y(5)=89.076
x(6)=-2.797;	y(6)=89.608
x(7)=-2.702;	y(7)=89.868
x(8)=-2.699;	y(8)=90.101
x(9)=-2.633;	y(9)=92.405
x(10)=-2.481;	y(10)=95.854
x(11)=-2.363;	y(11)=100.696
x(12)=-2.322;	y(12)=101.06
x(13)=-1.501;	y(13)=401.672
x(14)=-1.46;	y(14)=390.724
x(15)=-1.274;	y(15)=567.534
x(16)=-1.212;	y(16)=635.316
x(17)=-1.1;	y(17)=733.054
x(18)=-1.046;	y(18)=759.087
x(19)=-0.915;	y(19)=894.206
x(20)=-0.714;	y(20)=990.785
x(21)=-0.566;	y(21)=1090.109
x(22)=-0.545;	y(22)=1080.914
x(23)=-0.4;	y(23)=1122.643
x(24)=-0.309;	y(24)=1178.351
x(25)=-0.109;	y(25)=1260.531
x(26)=-0.103;	y(26)=1273.514
x(27)=0.01;	y(27)=1288.339
x(28)=0.119;	y(28)=1327.543
x(29)=0.377;	y(29)=1353.863
x(30)=0.79;	y(30)=1414.509
x(31)=0.963;	y(31)=1425.208
x(32)=1.006;	y(32)=1421.384
x(33)=1.115;	y(33)=1442.962
x(34)=1.572;	y(34)=1464.35
x(35)=1.841;	y(35)=1468.705
x(36)=2.047;	y(36)=1447.894
x(37)=2.2;	y(37)=1457.628


