时间有限，写了个简单的代码，可以初步算出大概的结果：

clear
set obs 1000000
gen d = .
gen x = .
local j = 1
forvalues x = 2(0.000001)3{
  qui replace x = `x' in `j'
  local d = `x' +`x'^2 +`x'^3 +`x'^4 +`x'^5 - 100
  qui replace d = abs(`d') in `j++'
  *dis in y "." _c
}
sort d
br
sort x
line d x   //有点动画的效果

非常感谢连老师的回答，一直在线等着的。可惜sata对这样的问题处理太复杂了，不适合进行课堂教学。但仍然感谢

连老师将stata这一工具发挥到了那是相当高的水平。
学习中。