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

% NIST certified values
target = zeros(dim,1);
target(1) = 1.6745063063E+00;
target(2) = -1.3927397867E-01;
target(3) = 2.5961181191E-03;
target(4) = -1.7241811870E-03;
target(5) = 2.1664802578E-05;

% start 1
ig1 = [2;-0.1;0.003;-0.001;0.00001]
% start 2
ig2 = [1.5;-0.15;0.0025;-0.0015;0.00002]

% points to fit
x = zeros(nb,1)
y = zeros(nb,1)
x(1)=9.65;	y(1)=0.0082
x(2)=10.74;	y(2)=0.0112
x(3)=11.81;	y(3)=0.0149
x(4)=12.88;	y(4)=0.0198
x(5)=14.06;	y(5)=0.0248
x(6)=15.28;	y(6)=0.0324
x(7)=16.63;	y(7)=0.042
x(8)=18.19;	y(8)=0.0549
x(9)=19.88;	y(9)=0.0719
x(10)=21.84;	y(10)=0.0963
x(11)=24;	y(11)=0.1291
x(12)=26.25;	y(12)=0.171
x(13)=28.86;	y(13)=0.2314
x(14)=31.85;	y(14)=0.3227
x(15)=35.79;	y(15)=0.4809
x(16)=40.18;	y(16)=0.7084
x(17)=44.74;	y(17)=1.022
x(18)=49.53;	y(18)=1.458
x(19)=53.94;	y(19)=1.952
x(20)=58.29;	y(20)=2.541
x(21)=62.63;	y(21)=3.223
x(22)=67.03;	y(22)=3.999
x(23)=71.25;	y(23)=4.852
x(24)=75.22;	y(24)=5.732
x(25)=79.33;	y(25)=6.727
x(26)=83.56;	y(26)=7.835
x(27)=87.75;	y(27)=9.025
x(28)=91.93;	y(28)=10.267
x(29)=96.1;	    y(29)=11.578
x(30)=100.28;	y(30)=12.944
x(31)=104.46;	y(31)=14.377
x(32)=108.66;	y(32)=15.856
x(33)=112.71;	y(33)=17.331
x(34)=116.88;	y(34)=18.885
x(35)=121.33;	y(35)=20.575
x(36)=125.79;	y(36)=22.32
x(37)=125.79;	y(37)=22.303
x(38)=128.74;	y(38)=23.46
x(39)=130.27;	y(39)=24.06
x(40)=133.33;	y(40)=25.272
x(41)=134.79;	y(41)=25.853
x(42)=137.93;	y(42)=27.11
x(43)=139.33;	y(43)=27.658
x(44)=142.46;	y(44)=28.924
x(45)=143.9;	y(45)=29.511
x(46)=146.91;	y(46)=30.71
x(47)=148.51;	y(47)=31.35
x(48)=151.41;	y(48)=32.52
x(49)=153.17;	y(49)=33.23
x(50)=155.97;	y(50)=34.33
x(51)=157.76;	y(51)=35.06
x(52)=160.56;	y(52)=36.17
x(53)=162.3;	y(53)=36.84
x(54)=165.21;	y(54)=38.01
x(55)=166.9;	y(55)=38.67
x(56)=169.92;	y(56)=39.87
x(57)=170.32;	y(57)=40.03
x(58)=171.54;	y(58)=40.5
x(59)=173.79;	y(59)=41.37
x(60)=174.57;	y(60)=41.67
x(61)=176.25;	y(61)=42.31
x(62)=177.34;	y(62)=42.73
x(63)=179.19;	y(63)=43.46
x(64)=181.02;	y(64)=44.14
x(65)=182.08;	y(65)=44.55
x(66)=183.88;	y(66)=45.22
x(67)=185.75;	y(67)=45.92
x(68)=186.8;	y(68)=46.3
x(69)=188.63;	y(69)=47
x(70)=190.45;	y(70)=47.68
x(71)=191.48;	y(71)=48.06
x(72)=193.35;	y(72)=48.74
x(73)=195.22;	y(73)=49.41
x(74)=196.23;	y(74)=49.76
x(75)=198.05;	y(75)=50.43
x(76)=199.97;	y(76)=51.11
x(77)=201.06;	y(77)=51.5
x(78)=202.83;	y(78)=52.12
x(79)=204.69;	y(79)=52.76
x(80)=205.86;	y(80)=53.18
x(81)=207.58;	y(81)=53.78
x(82)=209.5;	y(82)=54.46
x(83)=210.65;	y(83)=54.83
x(84)=212.33;	y(84)=55.4
x(85)=215.43;	y(85)=56.43
x(86)=217.16;	y(86)=57.03
x(87)=220.21;	y(87)=58
x(88)=221.98;	y(88)=58.61
x(89)=225.06;	y(89)=59.58
x(90)=226.79;	y(90)=60.11
x(91)=229.92;	y(91)=61.1
x(92)=231.69;	y(92)=61.65
x(93)=234.77;	y(93)=62.59
x(94)=236.6;	y(94)=63.12
x(95)=239.63;	y(95)=64.03
x(96)=241.5;	y(96)=64.62
x(97)=244.48;	y(97)=65.49
x(98)=246.4;	y(98)=66.03
x(99)=249.35;	y(99)=66.89
x(100)=251.32;	y(100)=67.42
x(101)=254.22;	y(101)=68.23
x(102)=256.24;	y(102)=68.77
x(103)=259.11;	y(103)=69.59
x(104)=261.18;	y(104)=70.11
x(105)=264.02;	y(105)=70.86
x(106)=266.13;	y(106)=71.43
x(107)=268.94;	y(107)=72.16
x(108)=271.09;	y(108)=72.7
x(109)=273.87;	y(109)=73.4
x(110)=276.08;	y(110)=73.93
x(111)=278.83;	y(111)=74.6
x(112)=281.08;	y(112)=75.16
x(113)=283.81;	y(113)=75.82
x(114)=286.11;	y(114)=76.34
x(115)=288.81;	y(115)=76.98
x(116)=291.08;	y(116)=77.48
x(117)=293.75;	y(117)=78.08
x(118)=295.99;	y(118)=78.6
x(119)=298.64;	y(119)=79.17
x(120)=300.84;	y(120)=79.62
x(121)=302.02;	y(121)=79.88
x(122)=303.48;	y(122)=80.19
x(123)=305.65;	y(123)=80.66
x(124)=308.27;	y(124)=81.22
x(125)=310.41;	y(125)=81.66
x(126)=313.01;	y(126)=82.16
x(127)=315.12;	y(127)=82.59
x(128)=317.71;	y(128)=83.14
x(129)=319.79;	y(129)=83.5
x(130)=322.36;	y(130)=84
x(131)=324.42;	y(131)=84.4
x(132)=326.98;	y(132)=84.89
x(133)=329.01;	y(133)=85.26
x(134)=331.56;	y(134)=85.74
x(135)=333.56;	y(135)=86.07
x(136)=336.1;	y(136)=86.54
x(137)=338.08;	y(137)=86.89
x(138)=340.6;	y(138)=87.32
x(139)=342.57;	y(139)=87.65
x(140)=345.08;	y(140)=88.1
x(141)=347.02;	y(141)=88.43
x(142)=349.52;	y(142)=88.83
x(143)=351.44;	y(143)=89.12
x(144)=353.93;	y(144)=89.54
x(145)=355.83;	y(145)=89.85
x(146)=358.32;	y(146)=90.25
x(147)=360.2;	y(147)=90.55
x(148)=362.67;	y(148)=90.93
x(149)=364.53;	y(149)=91.2
x(150)=367;	    y(150)=91.55
x(151)=371.3;	y(151)=92.2
