My OperaEllips
Saturday, November 15, 2008 3:40:55 AM
> Elips1:=(x,y,z)->(x^2/25)+(y^2/16)+(z^2/9)=1;
pertaam - tama kita definisikan dulu fungsinya. kita definisikan satu persatu dan gambar dulu satu persatu baru setelah itu digabung. kita menggbarnya dengan menggunakan perintah implicitplot3d.
> implicitplot3d((x^2/25)+(y^2/16)+(z^2/9)=1,x=-5..5,y=-4..4,z=-5..5);
Selanjutnya kita definisikan elips yang kedua dan kita gambarkan juga dengan menggunakan perintah implicitplot3d.
> Elips2:=(x,y,z)->(x^2/9)+(y^2/16)+(z^2/25)=1;
> implicitplot3d((x^2/9)+(y^2/16)+(z^2/25)=1,x=-5..5,y=-4..4,z=-5..5);
> _seed:=rand();
> randomX:=()->rand()/10.0^12;
> randomY:=()->rand()*.25/10.0^12;
> PasUr:={seq([randomX(),randomY()],k=1..1000)}:
> nDalamS:=0:
> for i from 1 to 1000 do
if PasUr [2]<=f(PasUr[1]) then nDalamS:=nDalamS+1 fi od;
Error, cannot determine if this expression is true or false: .1357769704e-1-f(.1081886648e-1) <= 0
> nDalamS;
>
pertaam - tama kita definisikan dulu fungsinya. kita definisikan satu persatu dan gambar dulu satu persatu baru setelah itu digabung. kita menggbarnya dengan menggunakan perintah implicitplot3d.
> implicitplot3d((x^2/25)+(y^2/16)+(z^2/9)=1,x=-5..5,y=-4..4,z=-5..5);
Selanjutnya kita definisikan elips yang kedua dan kita gambarkan juga dengan menggunakan perintah implicitplot3d.
> Elips2:=(x,y,z)->(x^2/9)+(y^2/16)+(z^2/25)=1;
> implicitplot3d((x^2/9)+(y^2/16)+(z^2/25)=1,x=-5..5,y=-4..4,z=-5..5);
> _seed:=rand();
> randomX:=()->rand()/10.0^12;
> randomY:=()->rand()*.25/10.0^12;
> PasUr:={seq([randomX(),randomY()],k=1..1000)}:
> nDalamS:=0:
> for i from 1 to 1000 do
if PasUr [2]<=f(PasUr[1]) then nDalamS:=nDalamS+1 fi od;
Error, cannot determine if this expression is true or false: .1357769704e-1-f(.1081886648e-1) <= 0
> nDalamS;
>
Unregistered user # Monday, September 28, 2009 3:16:56 AM
Putra Pratama # Saturday, November 14, 2009 5:23:21 AM
Originally posted by anonymous:
hemm.. gimana yak? soalnya aku juga udah jarang praktikum pake softwarenya. Jadi rada lupa gitu dah.. hehe:(