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

% NIST certified values
target = zeros(dim,1)
target(1) = 2.5906836021E+00
target(2) = 5.6177717026E-09
target(3) = -5.7701013174E-02

% start 1
ig1 = [2;0.0001;-0.01]
% start 2
ig2 = [2.5;0.000000005;-0.05]

% points to fit
x = zeros(nb,1)
y = zeros(nb,1)
x(1,1)=1;	x(1,2)=180;	y(1)=15
x(2,1)=1;	x(2,2)=180;	y(2)=17
x(3,1)=1;	x(3,2)=180;	y(3)=15.5
x(4,1)=1;	x(4,2)=180;	y(4)=16.5
x(5,1)=1;	x(5,2)=225;	y(5)=15.5
x(6,1)=1;	x(6,2)=225;	y(6)=15
x(7,1)=1;	x(7,2)=225;	y(7)=16
x(8,1)=1;	x(8,2)=225;	y(8)=14.5
x(9,1)=1;	x(9,2)=250;	y(9)=15
x(10,1)=1;	x(10,2)=250;	y(10)=14.5
x(11,1)=1;	x(11,2)=250;	y(11)=12.5
x(12,1)=1;	x(12,2)=250;	y(12)=11
x(13,1)=1;	x(13,2)=275;	y(13)=14
x(14,1)=1;	x(14,2)=275;	y(14)=13
x(15,1)=1;	x(15,2)=275;	y(15)=14
x(16,1)=1;	x(16,2)=275;	y(16)=11.5
x(17,1)=2;	x(17,2)=180;	y(17)=14
x(18,1)=2;	x(18,2)=180;	y(18)=16
x(19,1)=2;	x(19,2)=180;	y(19)=13
x(20,1)=2;	x(20,2)=180;	y(20)=13.5
x(21,1)=2;	x(21,2)=225;	y(21)=13
x(22,1)=2;	x(22,2)=225;	y(22)=13.5
x(23,1)=2;	x(23,2)=225;	y(23)=12.5
x(24,1)=2;	x(24,2)=225;	y(24)=12.5
x(25,1)=2;	x(25,2)=250;	y(25)=12.5
x(26,1)=2;	x(26,2)=250;	y(26)=12
x(27,1)=2;	x(27,2)=250;	y(27)=11.5
x(28,1)=2;	x(28,2)=250;	y(28)=12
x(29,1)=2;	x(29,2)=275;	y(29)=13
x(30,1)=2;	x(30,2)=275;	y(30)=11.5
x(31,1)=2;	x(31,2)=275;	y(31)=13
x(32,1)=2;	x(32,2)=275;	y(32)=12.5
x(33,1)=4;	x(33,2)=180;	y(33)=13.5
x(34,1)=4;	x(34,2)=180;	y(34)=17.5
x(35,1)=4;	x(35,2)=180;	y(35)=17.5
x(36,1)=4;	x(36,2)=180;	y(36)=13.5
x(37,1)=4;	x(37,2)=225;	y(37)=12.5
x(38,1)=4;	x(38,2)=225;	y(38)=12.5
x(39,1)=4;	x(39,2)=225;	y(39)=15
x(40,1)=4;	x(40,2)=225;	y(40)=13
x(41,1)=4;	x(41,2)=250;	y(41)=12
x(42,1)=4;	x(42,2)=250;	y(42)=13
x(43,1)=4;	x(43,2)=250;	y(43)=12
x(44,1)=4;	x(44,2)=250;	y(44)=13.5
x(45,1)=4;	x(45,2)=275;	y(45)=10
x(46,1)=4;	x(46,2)=275;	y(46)=11.5
x(47,1)=4;	x(47,2)=275;	y(47)=11
x(48,1)=4;	x(48,2)=275;	y(48)=9.5
x(49,1)=8;	x(49,2)=180;	y(49)=15
x(50,1)=8;	x(50,2)=180;	y(50)=15
x(51,1)=8;	x(51,2)=180;	y(51)=15.5
x(52,1)=8;	x(52,2)=180;	y(52)=16
x(53,1)=8;	x(53,2)=225;	y(53)=13
x(54,1)=8;	x(54,2)=225;	y(54)=10.5
x(55,1)=8;	x(55,2)=225;	y(55)=13.5
x(56,1)=8;	x(56,2)=225;	y(56)=14
x(57,1)=8;	x(57,2)=250;	y(57)=12.5
x(58,1)=8;	x(58,2)=250;	y(58)=12
x(59,1)=8;	x(59,2)=250;	y(59)=11.5
x(60,1)=8;	x(60,2)=250;	y(60)=11.5
x(61,1)=8;	x(61,2)=275;	y(61)=6.5
x(62,1)=8;	x(62,2)=275;	y(62)=5.5
x(63,1)=8;	x(63,2)=275;	y(63)=6
x(64,1)=8;	x(64,2)=275;	y(64)=6
x(65,1)=16;	x(65,2)=180;	y(65)=18.5
x(66,1)=16;	x(66,2)=180;	y(66)=17
x(67,1)=16;	x(67,2)=180;	y(67)=15.3
x(68,1)=16;	x(68,2)=180;	y(68)=16
x(69,1)=16;	x(69,2)=225;	y(69)=13
x(70,1)=16;	x(70,2)=225;	y(70)=14
x(71,1)=16;	x(71,2)=225;	y(71)=12.5
x(72,1)=16;	x(72,2)=225;	y(72)=11
x(73,1)=16;	x(73,2)=250;	y(73)=12
x(74,1)=16;	x(74,2)=250;	y(74)=12
x(75,1)=16;	x(75,2)=250;	y(75)=11.5
x(76,1)=16;	x(76,2)=250;	y(76)=12
x(77,1)=16;	x(77,2)=275;	y(77)=6
x(78,1)=16;	x(78,2)=275;	y(78)=6
x(79,1)=16;	x(79,2)=275;	y(79)=5
x(80,1)=16;	x(80,2)=275;	y(80)=5.5
x(81,1)=32;	x(81,2)=180;	y(81)=12.5
x(82,1)=32;	x(82,2)=180;	y(82)=13
x(83,1)=32;	x(83,2)=180;	y(83)=16
x(84,1)=32;	x(84,2)=180;	y(84)=12
x(85,1)=32;	x(85,2)=225;	y(85)=11
x(86,1)=32;	x(86,2)=225;	y(86)=9.5
x(87,1)=32;	x(87,2)=225;	y(87)=11
x(88,1)=32;	x(88,2)=225;	y(88)=11
x(89,1)=32;	x(89,2)=250;	y(89)=11
x(90,1)=32;	x(90,2)=250;	y(90)=10
x(91,1)=32;	x(91,2)=250;	y(91)=10.5
x(92,1)=32;	x(92,2)=250;	y(92)=10.5
x(93,1)=32;	x(93,2)=275;	y(93)=2.7
x(94,1)=32;	x(94,2)=275;	y(94)=2.7
x(95,1)=32;	x(95,2)=275;	y(95)=2.5
x(96,1)=32;	x(96,2)=275;	y(96)=2.4
x(97,1)=48;	x(97,2)=180;	y(97)=13
x(98,1)=48;	x(98,2)=180;	y(98)=13.5
x(99,1)=48;	x(99,2)=180;	y(99)=16.5
x(100,1)=48;	x(100,2)=180;	y(100)=13.6
x(101,1)=48;	x(101,2)=225;	y(101)=11.5
x(102,1)=48;	x(102,2)=225;	y(102)=10.5
x(103,1)=48;	x(103,2)=225;	y(103)=13.5
x(104,1)=48;	x(104,2)=225;	y(104)=12
x(105,1)=48;	x(105,2)=250;	y(105)=7
x(106,1)=48;	x(106,2)=250;	y(106)=6.9
x(107,1)=48;	x(107,2)=250;	y(107)=8.8
x(108,1)=48;	x(108,2)=250;	y(108)=7.9
x(109,1)=48;	x(109,2)=275;	y(109)=1.2
x(110,1)=48;	x(110,2)=275;	y(110)=1.5
x(111,1)=48;	x(111,2)=275;	y(111)=1
x(112,1)=48;	x(112,2)=275;	y(112)=1.5
x(113,1)=64;	x(113,2)=180;	y(113)=13
x(114,1)=64;	x(114,2)=180;	y(114)=12.5
x(115,1)=64;	x(115,2)=180;	y(115)=16.5
x(116,1)=64;	x(116,2)=180;	y(116)=16
x(117,1)=64;	x(117,2)=225;	y(117)=11
x(118,1)=64;	x(118,2)=225;	y(118)=11.5
x(119,1)=64;	x(119,2)=225;	y(119)=10.5
x(120,1)=64;	x(120,2)=225;	y(120)=10
x(121,1)=64;	x(121,2)=250;	y(121)=7.27
x(122,1)=64;	x(122,2)=250;	y(122)=7.5
x(123,1)=64;	x(123,2)=250;	y(123)=6.7
x(124,1)=64;	x(124,2)=250;	y(124)=7.6
x(125,1)=64;	x(125,2)=275;	y(125)=1.5
x(126,1)=64;	x(126,2)=275;	y(126)=1
x(127,1)=64;	x(127,2)=275;	y(127)=1.2
x(128,1)=64;	x(128,2)=275;	y(128)=1.2
