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

% NIST certified values
target = zeros(dim,1);
target(1) = 1.6657666537E-01;
target(2) = 5.1653291286E-03;
target(3) = 1.2150007096E-02;

% start 1
ig1 = [0.1;0.01;0.02]

% start 2
ig2 = [0.15;0.008;0.010]

% points to fit
x = zeros(nb,1)
y = zeros(nb,1)
x(1)=0.5;	y(1)=92.9;
x(2)=1;	y(2)=57.1;
x(3)=1.75;	y(3)=31.05;
x(4)=3.75;	y(4)=11.5875;
x(5)=5.75;	y(5)=8.025;
x(6)=0.875;	y(6)=63.6;
x(7)=2.25;	y(7)=21.4;
x(8)=3.25;	y(8)=14.25;
x(9)=5.25;	y(9)=8.475;
x(10)=0.75;	y(10)=63.8;
x(11)=1.75;	y(11)=26.8;
x(12)=2.75;	y(12)=16.4625;
x(13)=4.75;	y(13)=7.125;
x(14)=0.625;	y(14)=67.3;
x(15)=1.25;	y(15)=41;
x(16)=2.25;	y(16)=21.15;
x(17)=4.25;	y(17)=8.175;
x(18)=0.5;	y(18)=81.5;
x(19)=3;	y(19)=13.12;
x(20)=0.75;	y(20)=59.9;
x(21)=3;	y(21)=14.62;
x(22)=1.5;	y(22)=32.9;
x(23)=6;	y(23)=5.44;
x(24)=3;	y(24)=12.56;
x(25)=6;	y(25)=5.44;
x(26)=1.5;	y(26)=32;
x(27)=3;	y(27)=13.95;
x(28)=0.5;	y(28)=75.8;
x(29)=2;	y(29)=20;
x(30)=4;	y(30)=10.42;
x(31)=0.75;	y(31)=59.5;
x(32)=2;	y(32)=21.67;
x(33)=5;	y(33)=8.55;
x(34)=0.75;	y(34)=62;
x(35)=2.25;	y(35)=20.2;
x(36)=3.75;	y(36)=7.76;
x(37)=5.75;	y(37)=3.75;
x(38)=3;	y(38)=11.81;
x(39)=0.75;	y(39)=54.7;
x(40)=2.5;	y(40)=23.7;
x(41)=4;	y(41)=11.55;
x(42)=0.75;	y(42)=61.3;
x(43)=2.5;	y(43)=17.7;
x(44)=4;	y(44)=8.74;
x(45)=0.75;	y(45)=59.2;
x(46)=2.5;	y(46)=16.3;
x(47)=4;	y(47)=8.62;
x(48)=0.5;	y(48)=81;
x(49)=6;	y(49)=4.87;
x(50)=3;	y(50)=14.62;
x(51)=0.5;	y(51)=81.7;
x(52)=2.75;	y(52)=17.17;
x(53)=0.5;	y(53)=81.3;
x(54)=1.75;	y(54)=28.9;