function c=nagurney2_costs(x)  
c(1)=0.01*x(1)^2+x(1)+x(14)+100;
c(2)=0.01*x(2)^2+x(2)+0.5*x(11)+0.25*x(12)+0.25*x(14)+50;
c(3)=0.01*x(3)^2+x(3)+0.25*x(10)+x(21)+0.25*x(22)+150;
c(4)=0.01*x(4)^2+x(4)+x(10)+50;
c(5)=0.01*x(5)^2+x(5)+0.25*x(32)+0.25*x(33)+x(34)+100;
c(6)=0.01*x(6)^2+x(6)+x(32)+75;
c(7)=0.01*x(7)^2+x(7)+x(36)+75;
c(8)=0.01*x(8)^2+x(8)+0.5*x(27)+0.25*x(28)+0.25*x(36)+100;
c(9)=0.01*x(9)^2+x(9)+x(22)+100;
c(10)=0.01*x(10)^2+x(10)+0.5*x(3)+0.25*x(4)+0.5*x(22)+50;
c(11)=0.01*x(11)^2+x(11)+0.25*x(2)+x(13)+0.25*x(14)+50;
c(12)=0.01*x(12)^2+x(12)+x(2)+50;
c(13)=0.01*x(13)^2+x(13)+x(11)+50;
c(14)=0.01*x(14)^2+x(14)+0.25*x(2)+x(1)+0.25*x(11)+50;
c(15)=0.01*x(15)^2+x(15)+x(26)+50;
c(16)=0.01*x(16)^2+x(16)+0.25*x(19)+0.5*x(20)+0.25*x(26)+50;
c(17)=0.01*x(17)^2+x(17)+x(30)+50;
c(18)=0.01*x(18)^2+x(18)+0.5*x(30)+0.25*x(23)+x(24)+50;
c(19)=0.01*x(19)^2+x(19)+0.5*x(16)+x(25)+0.25*x(26)+50;
c(20)=0.01*x(20)^2+x(20)+x(16)+50;
c(21)=0.01*x(21)^2+x(21)+x(3)+50;
c(22)=0.01*x(22)^2+x(22)+0.25*x(3)+0.25*x(10)+50;
c(23)=0.01*x(23)^2+x(23)+x(17)+0.25*x(18)+0.25*x(30)+x(29)+100;
c(24)=0.01*x(24)^2+x(24)+x(18)+100;
c(25)=0.01*x(25)^2+x(25)+x(19)+100;
c(26)=0.01*x(26)^2+x(26)+x(20)+0.25*x(19)+0.25*x(16)+x(15)+75;
c(27)=0.01*x(27)^2+x(27)+0.25*x(28)+0.25*x(36)+x(35)+75;
c(28)=0.01*x(28)^2+x(28)+x(8)+50;
c(29)=0.01*x(29)^2+x(29)+x(23)+50;
c(30)=0.01*x(30)^2+x(30)+0.25*x(18)+0.5*x(17)+0.25*x(23)+50;
c(31)=0.01*x(31)^2+x(31)+x(33)+75;
c(32)=0.01*x(32)^2+x(32)+0.25*x(33)+x(6)+0.25*x(5)+75;
c(33)=0.01*x(33)^2+x(33)+0.5*x(32)+0.25*x(31)+0.25*x(5)+50;
c(34)=0.01*x(34)^2+x(34)+x(5)+50;
c(35)=0.01*x(35)^2+x(35)+x(27)+40;
c(36)=0.01*x(36)^2+x(36)+x(7)+0.25*x(8)+0.25*x(27)+75;
c=c';
