[maple]collecter des valeurs...

[Romain-des-Bois]

Bonjour à tous !

je suis vraiment pas très fort en maple et là, j'ai un petit problème (qui doit se résoudre très simplement je suis sûr !)

J'ai une fonction f (définie avant)

Je construis une fonction indicatrice :

Je génère deux nombres aléatoires x et y.
Je calcule f(x)
Si (Max(f(x),y) = f(x)) Alors Ind=1 Sinon Ind=0
Je lui demande de m'afficher Ind.

Chaque fois que j'exécute cette commande, maple m'affiche 1 ou 0 (c'est nickel)

La commande :
>x:=h():
y:=h():
If (max(f(x),y) = f(x)) Then Ind:=1 Else Ind:=0:
Ind;

(H est la procédure qui me génère mes nombres aléatoires)

Le problème, c'est que je voudrais répéter ça 100 ou 1000 fois ! Et je voudrais que maple me donne simplement le nombre de 1 qui sont sortis

Si quelqu'un a une piste, qu'il n'hésite pas !

Merci

Romain

PS : l'objectif de ce petit programme, c'est de calculer une intégrale par la méthode de Monte-Carlo

PPS : quelqu'un connait-il la commande maple pour donner le nombre de chiffres d'un nombre ?

Merci !
-----

[invitec053041c]

re bonjour !

Une petite boucle for ferait l'affaire je pense.

c:=0

for i from 1 to 100 do

x:=h():
y:=h():
If (max(f(x),y) = f(x))
then c:=c+1;

end for;
return c;
end proc;

EDIT: j'imagine que c'est un proc(f) non ?

[Romain-des-Bois]

Salut !

je vois ce que tu veux faire, tu ruses un peu (je veux dire par là qu'avec cette boucle, j'ai aucune idée de la position des 1 bon, je commence à faire mon difficile !)

je vais essayer pour voir si ça me convient !


Merci beaucoup !

Romain

[invitec053041c]

Je ne comprend pas, ton programme me sort que des 0 .(au passage j'ai oublié un end if )

Sinon, si tu veux la liste des 0 et des 1 sortis, tu peux créer une liste (et garder ton compteur en même temps aussi).

L:=NULL;

for i from 1 to 100 do

x:=h():
y:=h():
If (max(f(x),y) = f(x))
then L:=L,1
else L:=L,0

end if;
end for;

blabla
return [L]

[A voir en vidéo sur Futura]

[invitec053041c]

Mdrrr
Je pensais que h() était une fonction Maple! Je pouvais attendre longtemps des 1 .

[Romain-des-Bois]

Re-

j'ai apporté les deux-trois toutes petites modifications qui s'imposaient à ta boucle, et ça marche nickel ! En plus le temps de calcul est pas long, même jusqu'à 10000.

Finalement, je m'en fous un peu de la position des 0 (excuse moi de t'avoir fait réfléchir à ça pour rien !)
Etant donné que le temps de calcul est pas long, je vais appliquer la boucle jusqu'à 10000, peut-être 100000, et pas question d'avoir une liste de 100000 0 et 1 !!!

Déjà, J'ai des valeurs vraiment pas mal pour 10000 !


Merci pour ton précieux coup de main


Romain

EDIT : si tu veux éditer des nombres au hasard, tu utilises rand(). Ma procédure h() (comme hasard), c'était pour les contraindre entre 0 et 1.

[invitec053041c]

D'accord, c'est un plaisir .

Il ne faut pas avoir peur de le faire travailler le maple , le tout c'est d'y aller crescendo sinon on peut avoir de belles surprises...!
J'avais fait un petit programme qui vérifiait si jusqu'au nombre n la conjecture de Syrracuse était vérifiée .
Je l'ai fait aller loin, très loin! Puis je l'ai fait tourner pendant 1/2h...
C'est pas possible qu'elle ne soit pas vérifiée bon sang !

EDIT: d'accord, merci!

[Romain-des-Bois]

Tiens ! J'avais oublié cette conjecture (je ne me rappelais que du nom en fait ) Bon elle est vérifié alors ?

---

Si quelqu'un a une intégrale compliquée à calculer (avec des valeurs numériques, hein, comme en physique), qu'il n'hésite pas ! Ca me permettra de tester mon programme ! (erreur de moins de 0,005% pour cinq secondes de calcul )


Romain

[invitec053041c]

Tiens ! J'avais oublié cette conjecture (je ne me rappelais que du nom en fait ) Bon elle est vérifié alors ?

Jusqu'à quelques millions oui !

Si quelqu'un a une intégrale compliquée à calculer (avec des valeurs numériques, hein, comme en physique), qu'il n'hésite pas ! Ca me permettra de tester mon programme ! (erreur de moins de 0,005% pour cinq secondes de calcul )

Euh...

[FonKy-]

Jusqu'à quelques millions oui !





Euh...

loool

c'est bien jolie mais maple te donne directement le resultat

0.2760636949034663665639292115 074525604962180968216935737871 560562064308826076931536656382 330179699785

Edit: ah oui par contre j'ai essayer de la calculer pour 1337 chiffres et ca met un peu plus de temps

[FonKy-]

waouh il m'a pondu la reponse mais qu'au bout de un quart d'heure oO, ca donne quoi chez toi ?

0.2760636949034663665639292115 074525604962180968216935737871 560562064308826076931536656382 330179699784916254992008708778 052370772945726685780031995177 815235475152438139399562122060 397178064289288152720751369234 346980755818171283573800017094 240375949364211377039790307088 275603732957195663676262885238 595570341064151620687187083776 556270743274113996903373413813 924091968028741900258367396480 072828738517087333812506400813 873111720543518394488100286310 310018605296408442213447735786 766402361820569381978915415307 397575383391814907864759468992 551122834132451630888450067086 447578385428945166524827327179 495750619771083697615944741660 117124690996027873552850485650 907166409528205179697803887202 637659595633956766907526016422 662016769356366391644115979022 799117631048844127380001221436 329603049629987009979087128957 501241784134447958188503457411 408069665888074741659359313315 500835772504061164246044931149 816390942132084430179811992460 143459749270464313745335615198 796666595091243615841845313256 838995371178845064332935068971 494295068835760812358590157045 441317995924643818820487354597 215059451839480764281947232116 052331913161147000684471206098 380470138916751223449331539178 660244582132912075145550254809 619279349777344362165198599967 004363637837196502510581017295 508030364344436760603934287134 980904606658279953523078719875 3867338086705930339

Serait-il possible que tu partage ce programme ?

Cordialement, FonKy-

edit: j'ai pas mon meilleur pc sous la main aussi ^^

[invitec053041c]

, merci Fonky !
Je n'aurais pas pu dormir ce soir sans toutes ces décimales!

[FonKy-]

, merci Fonky !
Je n'aurais pas pu dormir ce soir sans toutes ces décimales!

C'est bizarre je viens de refaire le calcul a la main et je trouve que maple a fait une erreur sur la 1213e décimale.

ok je sors

[Romain-des-Bois]

Merci Ledescat pour ton exemple !

Je la programme demain, et je vous donnerai le résultat ! (il sera moins précis que maple tout cours)

Le programme sera disponible bientôt sur mon site (avec des explications et tout et tout...)


Romain