function c=nagurney4_costs(x)    
c(1)=5*x(1)+2*x(2)+500;
c(2)=4*x(2)+x(1)+200;
c(3)=3*x(3)+x(4)+350;
c(4)=6*x(4)+3*x(5)+400;
c(5)=6*x(5)+4*x(6)+600;
c(6)=7*x(6)+3*x(7)+500;
c(7)=8*x(7)+2*x(8)+400;
c(8)=5*x(8)+2*x(9)+650;
c(9)=6*x(9)+2*x(10)+700;
c(10)=4*x(10)+x(12)+800;
c(11)=7*x(11)+4*x(12)+650;
c(12)=8*x(12)+2*x(13)+700;
c(13)=7*x(13)+3*x(18)+600;
c(14)=8*x(14)+3*x(15)+500;
c(15)=9*x(15)+2*x(14)+200;
c(16)=8*x(16)+5*x(12)+300;
c(17)=7*x(17)+2*x(15)+450;
c(18)=5*x(18)+x(16)+300;
c(19)=8*x(19)+3*x(17)+600;
c(20)=6*x(20)+x(21)+300;
c(21)=4*x(21)+x(22)+400;
c(22)=6*x(22)+x(23)+500;
c(23)=9*x(23)+2*x(24)+350;
c(24)=8*x(24)+x(25)+400;
c(25)=9*x(25)+3*x(26)+450;
c(26)=7*x(26)+4*x(27)+300;
c(27)=8*x(27)+3*x(28)+500;
c(28)=7*x(28)+3*x(29)+650;
c(29)=3*x(29)+x(30)+450;
c(30)=7*x(30)+2*x(31)+600;
c(31)=4*x(31)+x(32)+750;
c(32)=8*x(32)+3*x(33)+650;
c(33)=9*x(33)+2*x(34)+750;
c(34)=8*x(34)+2*x(35)+650;
c(35)=6*x(35)+x(36)+500;
c(36)=8*x(36)+2*x(37)+450;
c(37)=7*x(37)+3*x(38)+700;
c(38)=8*x(38)+4*x(39)+800;
c(39)=7*x(39)+3*x(40)+550;
c(40)=8*x(40)+3*x(41)+600;
c(41)=8*x(41)+4*x(42)+750;
c(42)=5*x(42)+x(43)+350;
c(43)=7*x(43)+2*x(45)+400;
c(44)=8*x(44)+3*x(46)+500;
c(45)=9*x(45)+6*x(47)+350;
c(46)=6*x(46)+2*x(48)+400;
c(47)=3*x(47)+x(49)+800;
c(48)=7*x(48)+2*x(50)+700;
c(49)=6*x(49)+x(50)+600;
c(50)=8*x(50)+2*x(51)+500;
c(51)=7*x(51)+3*x(52)+600;
c(52)=6*x(52)+2*x(53)+700;
c(53)=9*x(53)+3*x(54)+650;
c(54)=7*x(54)+4*x(55)+700;
c(55)=8*x(55)+3*x(56)+250;
c(56)=6*x(56)+x(57)+350;
c(57)=5*x(57)+2*x(58)+200;
c(58)=7*x(58)+x(59)+250;
c(59)=6*x(59)+3*x(60)+400;
c(60)=4*x(60)+x(61)+300;
c(61)=7*x(61)+2*x(62)+400;
c(62)=5*x(62)+x(63)+200;
c(63)=4*x(63)+2*x(64)+100;
c(64)=5*x(64)+100;
c(65)=4*x(65)+200;
c(66)=7*x(66)+300;
c=c';