Annäherung an Pi mit Monte carlo, paintbox

[Borgschüler]

Guten Tag!

Habe morgen matura, habe ein Problem mit unserem Monte carlo Beispiel!

Und zwar schaut der Quelltext für mich ganz richtig aus, aber es wird nur die Hälfte von der Paintbox belegt und ich weis absolut nicht warum!

aber seht selbst:

Code:

procedure TForm1.Button1Click(Sender: TObject);
var i, count,x,y : integer;
begin
randomize;
count := 0;
for i := 1 to 100000 do
 begin
  x:= random(200);
  y:= random(200);
   if (sqrt(x*x+y*y) < 200) then
        paintbox1.canvas.pixels[x,y] := clred
     else
        paintbox1.canvas.pixels[x,y] := clgreen;
 end;
end;

Also die paintbox hat width/length jeweils 200. Das Problem ist nur, das Programm belegt nur exakt linke Hälfte, die eig ganz richtig ausschaut. Wenn ich duie paintbox 100 mal 100 mache (length width), dann scheint es zu funktionieren, problem ist nur dass es dann plötzlich kein richtiges quadrat mehr ist sondern ein rechteck.
Angefangen hab ich mit length/width 400/400, dann hate ich links einen streifen wo es funktioniert und rechts.

Ich dreh langsam durch. Kann es sein das width/length nicht die Pixel angibt?

Bitte dringend um Lösungsvorschläge und bitte nicht kleine Hinweise, damit ich dann selber draufkommen muss, weil ich probier schon verdammt lang um und komm einfach nicht auf die lösung... und morgen ist ja matura...


Pi ausgerechnet hab ich noch nicht, aber das dürfte dann kein problem sein.

[einbeliebigername]

Hallo,

nimm mal lieber ein TImage statt dem TPaintBox. Bei einer PaintBox verschwindet das was man nicht im OnPaint zeichnet wieder, wenn es neugezeichnet werden soll. Das könnte schon dein Problem sein, weil irgend was die PaintBox zum teilweisen Neuzeichnen bringt.

Einbeliebigername.

[Klaus01]

Hallo,

bitte randomize nur einmal (beim Programmstart) aufrufen.

Was bezweckst du mit?

Delphi-Quellcode:

if sqrt(x*x + y*y) < 200 then

Grüße
Klaus

[Borgschüler]

hmm in der schule haben wir das immer mit der paintbox gemacht.

randomize wird im button eh nur einmal aufgerufen. wenn ich es bei formcreate reinschreibe, ändert isch leider garnichts. er macht irgendwie immer nur die hälfte von der paintbox. es sei denn du machst die paintbox richtig klein, aber das kanns auch nicht sein!

[Klaus01]

Delphi-Quellcode:

 if (sqrt(x*x+y*y) < 200)

->

Delphi-Quellcode:

if sqrt(sqr(x)+sqr(y)) <= sqrt(sqr(paintbox.width)+(sqr(paintbox.heigth)) then
  //im radius
else
 //asserhalb des Radius

Grüße
Klaus

[einbeliebigername]

Hallo,

Zitat:

Zitat von Klaus01 (Beitrag 1146639)

Delphi-Quellcode:

if sqrt(sqr(x)+sqr(y)) <= sqrt(sqr(paintbox.width)+(sqr(paintbox.heigth)) then

Das ist aber nicht im Sinne des Erfinders. Denn damit ist der Radius des Kreis so groß wie die Diagonale der PaintBox. Das hat dann die Folge, dass alle Punkte, die man sieht, innerhalb des Kreises liegen. Besser ist das:

Delphi-Quellcode:

  R:= Min(paintbox.width, paintbox.heigth)
  ...
    x:= random(R);
    y:= random(R);
    if (sqrt(x*x+y*y) < R) then
  ...

Einbeliebigername.

[TiGü]

Wie soll es denn richtig aussehen?
Wenn ich deinen Code ausführe, dann kommt das bei rum:

[guinnes]

Zitat:

Zitat von TiGü (Beitrag 1146642)

Wie soll es denn richtig aussehen?

Genau so ! : http://de.wikipedia.org/wiki/Monte-Carlo-Simulation

[Borgschüler]

tigü ja genau, so soll es aussehen!
tigü, hast du meinen quelltext imn ersten beitrag eins zu eins übernommen?
hmm, dann probier ich mal, den radius mit paintbox width auszudürcken!

EDIT: Ich werd noch wahnsinnig, jetzt hab ich < paintbox1.width und es geht immer noch nicht!


übrigens:

als zum berechnen von pi dachte ich mir

K/Q = 200²Pi/200²

erklärung: Verhältnis der punkte von Kreis und Quadrat(gesamtes) = verhältnis von Viertelkreisgleichung mit festgelegtem radius und fläche von quadrat

und daraus drück ich mir aus:

K/Q mal 4 = Pi

[Borgschüler]

so so schaut das aus bei mir, siehe anhang