Modérateur : Politburo
Exact, j'ai oublié une ligne :Marcus von Cube a écrit :Mémoire R n'est jamais modifié. Il manque une instruction comme...
Code : Tout sélectionner
R+List 2[J]->Z
Int(Z/M)->R
Z-MR->List 2[J]:NextCode : Tout sélectionner
StoPic "Pics\Fib"Euh c'est quoi l'intérêt de cette formule ???Gilles59 a écrit :Gege, juste pour l'anecdote. En remplaçant 1 par 'n' en (2,1) dans la matrice, FIBO(400) donne :
200.n^200 +1333300.n^199+ 2666333 etc...
Aucun si ce n'est voir comment la HP50 calcule et que c'est tres lourd lol ...gege a écrit :Exact, j'ai oublié une ligne :Marcus von Cube a écrit :Mémoire R n'est jamais modifié. Il manque une instruction comme...Mea culpaCode : Tout sélectionner
R+List 2[J]->Z Int(Z/M)->R Z-MR->List 2[J]:Next
Sympa ça :Code : Tout sélectionner
StoPic "Pics\Fib"Euh c'est quoi l'intérêt de cette formule ???Gilles59 a écrit :Gege, juste pour l'anecdote. En remplaçant 1 par 'n' en (2,1) dans la matrice, FIBO(400) donne :
200.n^200 +1333300.n^199+ 2666333 etc...
G.E. perplexe
Mais le fichier qu'on trouve quand la calculatrice est branché sur le PC ne sert à rien. Le format est de type g3p. J'ai pas essayé de le analyser.gege a écrit : Sympa ça :Code : Tout sélectionner
StoPic "Pics\Fib"
Code : Tout sélectionner
« "" "" -> Ogr Zero @ Ogr : Objet graphique - Zero : Ligne remplie de 0
« TIME @ Heure de lancement du prog
1. 70. START "0" 'Zero' STO+ NEXT @ Génére ligne de pixels vides. 1 quartet = 4 pixels. 70 quartets = 280 pixels
1 1 @ Les 2 premiers nombres de Fibonacci
1. 398. START DUP2 + NEXT @ Génére les 400 premiers nbres de Fibonacci sur la pile en ~ 1 sec
1. 400. START @ la boucle dépile tout et génère le code hexa pour affichage binaire.
"" SWAP @ Initialise graphique de la ligne
WHILE 16 IDIV2 DUP2 OR REPEAT @ Conversion hexadécaimale au format graph GROB
"084C2A6E195D3B7F" SWAP 1 + @ nécessaire car le premier pixel d'un octet est le bit de
DUP SUB ROT + SWAP @ poids fort de l'octet : 1000000b pour afficher le premier pixel
END @ et pas 00000001b comme on pourait le croire
DROP2
Zero SWAP DUP SIZE 70. @ Calcul la longueur de la ligne à afficher
SWAP - 1. + SWAP REPL @ et rempli le début avec des 0
'Ogr' STO+ @ Concaténe dans l'objet graphique
NEXT
"GROB 280 400 " Ogr + OBJ-> @ Génére l'objet graphique sous forme de chaine de caractères
TIME SWAP @ Heure de fin et place le graph au niveau 1 de la pile
»
»
Code : Tout sélectionner
'ProgramMode:RUN
32->P
ClrGraph:AxesOff:LabelOff
2^P->M:187 Int/ P+1->Q
Q->Dim List 1
Q->Dim List 2
1->List 2[1]
"DEBUT"?->N
For 1->I To N+378
List 1+List 2->List 3
List 2->List 1
List 3->List 2
0->R:For 1->J To Q
List 2[J]+R->Z
Z Int/ M->R
Z-MR->List 2[J]:Next
If N<=I:Then
For 1->J To Q
List 2[J]->Z:I+1-N->C
PJ->K:K->L:K>187=>187->L
For K-P+1->R To L
Z/2->Z
If 0<Frac Z
Then Black PxlOn R,C
'_Else __Yellow __PxlOn _R,C
IfEnd
Int Z->Z:
Z=0=>Break
Next:Next
IfEnd
Next
StoPict "Pict\Fib"
Code : Tout sélectionner
'ProgramMode:RUN
32->P
ClrGraph:AxesOff:LabelOff
2^P->M:64/P->Q
Q->Dim List 1
Q->Dim List 2
1->List 2[1]
"DEBUT"?->N
For 1->I To N+126
List 1+List 2->List 3
List 2->List 1
List 3->List 2
0->R:For 1->J To Q
List 2[J]+R->Z
Int (Z/M)->R
Z-MR->List 2[J]:Next
If N<=I
Then For 1->J To Q
List 2[J]->Z:I+1-N->C
PJ->K:K->L:K>63=>63->L
For K-P+1->R To L
Z/2->Z
0<Frac Z=>PxlOn R,C
Int Z->Z
Z=0=>Break
Next:Next
IfEnd
Next
StoPict 1
Code : Tout sélectionner
'ProgramMode:RUN
24->P
ClrGraph
ViewWindow 1,79,0,47,1,0
2^P->M:2->Q
{0,0}->List 1
{1,0}->List 2
"DEBUT"?->N
For 1->I To N+78
List 1+List 2->List 3
List 2->List 1
List 3->List 2
0->R:For 1->J To Q
List 2[J]+R->Z
Int (Z/M)->R
Z-MR->List 2[J]:Next
If N<=I
Then For 1->J To Q
List 2[J]->Z:I+1-N->C
PJ->K:K->L:K>47=>47->L
For K-P+1->R To L
Z/2->Z
0<Frac Z=>Plot C,R
Int Z->Z
Z=0=>Break
Next:Next
IfEnd
Next
Intéressant! mais j'ai pas trop compris l'"astuce des pixels" :Ogege a écrit :Bonjour,
Avec l'astuce des pixels j'obtiens le graphe de 1 à 127 en 2 minute 10 secondes sur la Graph 85.
J'ai piqué l'idée de Marcus de faire le graphe horizontalement, on voit bien mieux.
Le programme ici ce soir, 144 octets tout mouillé et 5 variables (pas de liste) !!
Cool - à plusieurs on trouve des trucs.
G.E.
Code : Tout sélectionner
«
1 1
1. 131. START DUP2 + NEXT
#0d #130d FOR i
#0d 'j' STO
WHILE 2 IDIV2 DUP2 OR REPEAT
IF THEN i j 2 ->LIST PIXON END
#1d 'j' STO+
END
DROP2
NEXT
»
Donc on ne stocke rien, on allume des pixels et pour savoir lesquels doivent l'être, on teste ceux des deux lignes au-dessus.C.ret a écrit :En effet, la relation de récurence F(n)=F(n-1) + F(n-2) signifie que chaque ligne du graphique peut être générée à partir des deux lignes précèdentes. Et oui.
Code : Tout sélectionner
ClrGraph:AxesOff
PxlOn 1,1:PxlOn 1,2
For 3->I To 127
Min({63,.7I})->Q
0->R:For 1->J To Q
R+PxlTest(J,I-2)+PxlTest(J,I-1)->P:P/2->P
0<Frac P=>PxlOn J,I
Int P->R:Next
Next:StoPic 14En fait, les temps explosent pour les "grandes tailles" de dessin à cause de la conversion hexadécimale qui nécessite de diviser des entiers long et visiblement c'est (relativement) très lent sur HP49/50. Je suppose que c'est du à la façon dont les entiers 'précision infini' sont codé dans la machine (par contre les sommes sont très très rapide).gege a écrit :Hello,
2 minutes c'est bien ! J'avais l'impression que c'était beaucoup plus d'après tes dires.
G.E.
Gege, Je converti tel quel. C'est quoi le .7I dans la fonction Min ? 0.7 * i ?gege a écrit :Hello,
Programme :144 octets.Code : Tout sélectionner
ClrGraph:AxesOff PxlOn 1,1:PxlOn 1,2 For 3->I To 127 Min({63,.7I})->Q 0->R:For 1->J To Q R+PxlTest(J,I-2)+PxlTest(J,I-1)->P:P/2->P 0<Frac P=>PxlOn J,I Int P->R:Next Next:StoPic 14
Voilà,
G.E.
Min() c'est minimum d'une liste. Le facteur .7 est une "chiffre magique". C'est l'inclinaison de la droite qui limite le graphe.Gilles59 a écrit : C'est quoi le .7I dans la fonction Min ? 0.7 * i ?
Min comme minimum je suppose...
Code : Tout sélectionner
'ProgramMode:RUN
ClrGraph:AxesOff
Black PxlOn 1,1:Black PxlOn 1,2
For 3->I To 379
Min({187,.7I})->Q
0->R:For 1->J To Q
(R+PxlTest(J,I-2)+PxlTest(J,I-1))/2->P
0<Frac P=>Black PxlOn J,I
Int P->R:Next
Next
StoPict "Pict\Fib2"Code : Tout sélectionner
(R+(0<PxlTest(J,I-2))+(0<PxlTest(J,I-1)))/2->P
J'attends la prochaine version du OS. Code : Tout sélectionner
«
TIME
ERASE {#0d,#0d} PVIEW
{#1d,#1d} PIXON
{#2d,#1d} PIXON
#3d #127d FOR i
#0d
i .7 *
#1d SWAP FOR j
i #2d - j 2 ->LIST PIX? i #1d - j 2 ->LIST PIX? + + DUP
IF #1d AND #0d <> THEN i j 2 ->LIST PIXON END
ASR
NEXT
DROP
NEXT
TIME
»Code : Tout sélectionner
IF #1d AND #0d ‹ THENCode : Tout sélectionner
IF #1d AND #1d == THENC'est le signe différent (un égal barré) en fait mais il n'existe pas sur PC et j'ai oublié de le corriger. J'ai modifié au dessusMarcus von Cube a écrit :Cela ne marche pas. J'ai le changé comme ça :Code : Tout sélectionner
IF #1d AND #0d ‹ THENCode : Tout sélectionner
IF #1d AND #1d == THEN
Comme je transfère avec l'émulateur et fait un copier-coller,je dois faire les modifs en collant ici Non, parce-que j'ai aussi changé le chiffre : #0d =/= contre #1d ==Gilles59 a écrit : Avec == tu dois l'avoir en inversion vidéo par rapport à moi
C'est interessant.Gilles59 a écrit :J'utilise HPUserEdit v5 qui est vraiment excellent et a un lien direct avec l'émulateur HP50G.
Code : Tout sélectionner
%%HP: T(3)A(R)F(.);
\<<
TIME
PICT SIZE
#1d - SWAP #1d -
\-> t h w
\<<
ERASE
{ #0d #0d } PVIEW
{ #0d #0d } PIXON
{ #1d #0d } PIXON
#2d w
FOR i
#0d
#0d
i B\->R .7 *
h B\->R
MIN R\->B
FOR j
i #2d - j 2. \->LIST PIX? +
i #1d - j 2. \->LIST PIX? +
IF DUP #1d AND #1d == THEN
i j 2. \->LIST PIXON
END
ASR
NEXT
DROP
NEXT
PICT RCL 'FIBPICT' STO
TIME t -
\->HMS 100 * 2 RND
PICTURE
\>>
\>>