Pour illustrer le calcul d'une intégrale avec la méthode de Monte-Carlo

Version sans commentaires, utilisable avec profit seulement avec les commentaires faits au cours MAT1681 et après avoir assisté à sa construction pas à pas.

dessintir[f_,base_,hauteur_,d_,n_]:= Module[{x,b1,b2,h1,h2,tir,graphef},
    b1 = base[[1]]; b2 = base[[2]]; h1 = hauteur[[1]]; h2 = hauteur[[2]];
    tir = Graphics[{PointSize[d], Table[Point[{Random[Real,{b1,b2}],Random[Real,{h1,h2}]}],{n}]}];
    graphef=Plot[f[x],{x,b1,b2},DisplayFunction\[Rule]Identity];
    res=Show[graphef,tir,DisplayFunction\[Rule]$DisplayFunction]
    ]

Exemple : Dessin du graphe de la fonction qui calcule le carré de son argument ( fonction écrite comme la fonction pure #^2 &) et d'un tir de 1000 points (de "diamètre relatif" d = 0.01), le tout dans la cible [0,2] x [0,4] :
dessintir[#^2 &,{0,2},{0,4},.01,1000]

Retour à ma page personnelle



©jmt : création le 031015