*
* asymmetry.kumac 
* Estimateur pour l'asymmetrie angulaire
*
*                               Markus 04/01

a=0.03
*                               Parametre alpha
ne=10001
v/cr x([ne])
si n=array([ne],1#[ne])

fun1 10 fasym.f 50 -1. 1. ' '

do i=1,[ne]
 v/in x([i]) $CALL('HRNDM1(10)')
enddo

*
* Estimation de alpha a partir de la moyenne
*
si sx=4.*sumv(x)/n+.3
*
* Estimation avec l'observable "asymmetrie"
*
si p=x gt 0.
si sp=8./3.*(2.*sumv(p)/n -1.)+.1
*
* Optimal observable g(x)
*
si g=x/(1+x*x)
*
* Esperance de g(x)
si sg=16./(3.*(4-3.14159))*sumv(g)/n-.1
si gg=g**3
*
* Esperance de g(x)^3
*
si sgg=128./(3.*(3.*3.14159-8))*sumv(gg)/n-.3

*
* Presentation graphique
*

for/fil 66 asym_003.ps
meta 66 -111
title 'Estimateurs pour asymmetrie angulaire'
zone 1 2
opt nstat
*
1d 11 'Echantillon' 50 -1. 1.
v/hfill x 11; norm 11 25.; h/pl 11 ep1
fun1 12 fasym.f 50 -1. 1. s
itx 0.6 0.15 '[a]=0.03'
*
opt logx; null 40. [ne] [a]-.6 [a]+.6
set hcol 1; v/dr sx ! s
set plci 1;dline 40. [ne] [a]+.3 [a]+.3
set hcol 2; v/dr sp ! s
set plci 2;dline 40. [ne] [a]+.1 [a]+.1
set hcol 3; v/dr sg ! s
set plci 3;dline 40. [ne] [a]-.1 [a]-.1
set hcol 4; v/dr sgg ! s
set plci 4;dline 40. [ne] [a]-.3 [a]-.3
set hcol; set plci; opt linx
itx [ne]*1.1 [a]+.3 'moyenne'
itx [ne]*1.1 [a]+.1 'asymmetrie'
itx [ne]*1.1 [a]-.1 'E(g(x))'
itx [ne]*1.1 [a]-.3 'E(g^3!(x))'
atit 'Taille echantillon' 'Estimateur pour [a]'
close 66