function c=nagurney1_costs(x)    
c(1)=0.00005*x(1)^4+5*x(1)+2*x(2)+500;
c(2)=0.00003*x(2)^4+4*x(2)+x(1)+200;
c(3)=0.00005*x(3)^4+3*x(3)+x(4)+350;
c(4)=0.00003*x(4)^4+6*x(4)+3*x(5)+400;
c(5)=0.00006*x(5)^4+6*x(5)+4*x(6)+600;
c(6)=7*x(6)+3*x(7)+500;
c(7)=0.00008*x(7)^4+8*x(7)+2*x(8)+400;
c(8)=0.00004*x(8)^4+5*x(8)+2*x(9)+650;
c(9)=0.00001*x(9)^4+6*x(9)+2*x(10)+700;
c(10)=4*x(10)+x(12)+800;
c(11)=0.00007*x(11)^4+7*x(11)+4*x(12)+650;
c(12)=8*x(12)+2*x(13)+700;
c(13)=0.00001*x(13)^4+7*x(13)+3*x(18)+600;
c(14)=8*x(14)+3*x(15)+500;
c(15)=0.00003*x(15)^4+9*x(15)+2*x(14)+200;
c(16)=8*x(16)+5*x(12)+300;
c(17)=0.00003*x(17)^4+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)=0.00003*x(20)^4+6*x(20)+x(21)+300;
c(21)=0.00004*x(21)^4+4*x(21)+x(22)+400;
c(22)=0.00002*x(22)^4+6*x(22)+x(23)+500;
c(23)=0.00003*x(23)^4+9*x(23)+2*x(24)+350;
c(24)=0.00002*x(24)^4+8*x(24)+x(25)+400;
c(25)=0.00003*x(25)^4+9*x(25)+3*x(26)+450;
c(26)=0.00006*x(26)^4+7*x(26)+8*x(27)+300;
c(27)=0.00003*x(27)^4+8*x(27)+3*x(28)+500;
c(28)=0.00003*x(28)^4+7*x(28)+3*x(27)+650;
c=c';
