STATA入门10 随机模拟(3)

}

end

simulate " pf" z,reps(10000)

sum *============================end====================================

结果约为： 14102.65

clear

mat A=J(11,11,.)

STATA入门

forvalue i=1/11 {

forvalue j=1/11 {

if `i'>=`j'{

mat A[`i',`j']=1000*(9+`i')

}

else if `i'<`j' {

mat A[`i',`j']=500*(18+`i'+`j')

}

}

}

mat list A

mat a=J(1,11,1)

mat C=a*A*a'/121 mat list C

* 15909.091

*当进货量和销货量取1位小数时

clear

mat A=J(101,101,.)

forvalue i=1/101 {

forvalue j=1/101 {

if `i'>=`j'{

mat A[`i',`j']=1000*(10+(`i'-1)/10)

}

else if `i'<`j' {

mat A[`i',`j']=500*(20+(`i'-1)/10+(`j'-1)/10)

}

}

}

mat a=J(1,101,1)

mat C=a*A*a'/10201 mat list C

15841.584

*============================end==================================== 请思考，上述程序有什么问题？为什么答案有如此大的差距？