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/