par DamDam » 01 Juin 2007, 14:38
Voici le programme:
nd=NormalDistribution[60632.8136,24496.0932]
ld1=LogNormalDistribution[10.9575, 0.3742]
a=.
k=CDF[nd,x]
l=CDF[ld1,x- a]
"Calcul des aires"
a =.;
i=Integrate[k/. x->b, {b,-1000000, 1000000}]
j=Integrate[l/. x-> c, {c,-1000000, 1000000}]
réponses:
i = 939367.
j = If[Re[a]>=1.10*10^6 Re[a]<= -1*10^6 Im[a]=!0,1.*10^6+30772.2Erf[-20.9704+1.88965 Log[-1.*10^6-1. a]] + (500000. +0.5a) Erf[-20.7058 + 1.88965 Log[-1.*10^6-1.a]]-30772.2 Erf[-20.9704+1.88965 Log[1.*10^6-1.a]], Integrate[0.5+0.5Erf[-20.7058+1.88965 Log[-1.a+c]], {c, -10000000, 10000000}, Assumptions -> Im[a]==0&&-1.*10^6
Le but est alors que la différence entre i et j soit nulle.