;
; create data set ( 3 gaussian + noise)
;
sigmas = [10.0, 50, 10.0]/( 2.0*sqrt(2*alog(2)) )
x=makeArray1(100,0,100) 
y1= 5.*exp( -0.5*( (x-10.)/sigmas[0] )^2 ) 
y2= 10.*exp( -0.5*( (x-50.)/sigmas[1] )^2 )
y3= 30.*exp( -0.5*( (x-70.)/sigmas[2] )^2 )
ynoise =  0.5*randomu(seed,100)
y=y1+y2+y3+ynoise
;
; write data to a file (not needed)
;
OpenW,unit,'fakedata.dat',/Get_Lun
FOR i=0,99 DO printf,unit,x[i],y[i]
Free_Lun,unit
itmp = Dialog_Message('File fakedata.dat written to disk.')