Dear fellows !
/first thanks to Neil from yesterday/
please consider the following as "Pseudocode"
clear
set obs 100
gen d =.
gen d1 =.
forvalues i = 1(1)100 {
local i = `i' * 0.1     <-  0.1 *1,2,3,4,..... *-multplikation-operator
range x 1(1)i          <-   range until i-value                       
gen double y=exp(-(`i'-x)^2/2)*0.1/(x + 0.1)    <- y as function of i-value
integ y x,gen(Sy)    <- integrate until i-value
*save new y- and Sy- value in two columns :
replace d[_n] = y   
replace d1[_n] = Sy
drop Sy y
}
please any advice that it works correctly.
      sincerely aschbacher andreas
-- 
+++ NEU bei GMX und erstmalig in Deutschland: T�V-gepr�fter Virenschutz +++
100% Virenerkennung nach Wildlist. Infos: http://www.gmx.net/virenschutz
*
*   For searches and help try:
*   http://www.stata.com/support/faqs/res/findit.html
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/