%-- 02/07/13 07:03 --%
mkdir aula3
cd aula3
edit triangulo.m
[ x, y ] = triangulo( 1e5 );
plot( x, y, '.' )
plot( x, y, '.', 'markersize', 1 )
tic, [ x, y ] = triangulo( 1e5 ); toc
tic, [ x, y ] = triangulo( 1e6 ); toc
plot( x, y, '.', 'markersize', 1 )
tic, [ x, y, r ] = triangulo( 1e6 ); toc
figure
hold on
plot( x(r==1), y(r==1), '.b', 'markersize', 1 )
plot( x(r==2), y(r==2), '.r', 'markersize', 1 )
plot( x(r==3), y(r==3), '.g', 'markersize', 1 )
edit plota_regra
plota_regra( x, r, r )
plota_regra( x, y, r )
plota_regra( x, r, r )
whos
r(1:15)
r(1:10)
r_a = [ 0; r(2:end-1) ];
whos
r(1:10)
r_a(1:10)
[ r(1:10 r_a(1:10) ]
[ r(1:10) r_a(1:10) ]
r_a = [ 0; r(1:end-1) ];
[ r(1:10) r_a(1:10) ]
plot( x((r==3).*(r_a==3)), y((r==3).*(r_a==3)), '.k', 'markersize', 1 )
plot( x((r==3)&(r_a==3)), y((r==3)&(r_a==3)), '.k', 'markersize', 1 )
plot( x((r==3)&(r_a==2)), y((r==3)&(r_a==2)), '.m', 'markersize', 1 )
edit triangulo.m
help cumsum
[ x, y, r ] = triangulo( 1e5 );
plota_regra( x, r, r )
plota_regra( x, y, r )
[ x, y, r ] = triangulo( 1e5 );
help nargin
[ x, y, r ] = triangulo( 1e5 );
plota_regra( x, y, r )
[ x, y, r ] = triangulo( 1e5, [ .5 .25 .25] );
plota_regra( x, y, r )
[ x, y, r ] = triangulo( 1e5, [ .8 .1 .1] );
plota_regra( x, y, r )
edit aleatorio_especial
aleatorio_especial( [.5 .25 .25] )
aleatorio_especial( [.5 .25 .2 .05] )
aleatorio_especial( [.5 .25 .2 .04 .01] )
aleatorio_especial( [.5 .25 .2 .04 .01], 10 )
r = aleatorio_especial( [.5 .25 .2 .04 .01], 1e4 );
hist( r )
hist( r, 1:5 )
help hist
[h]=hist( r, 1:5 )
r = aleatorio_especial( [.5 .25 .2 .04 .01], 1e4 );
[h]=hist( r, 1:5 )
r = aleatorio_especial( [.5 .25 .2 .04 .01], 1e4 );
[h]=hist( r, 1:5 )
r = aleatorio_especial( [.5 .25 .2 .04 .01], 1e4 );
[h]=hist( r, 1:5 )
r = aleatorio_especial( [.5 .25 .2 .04 .01], 1e4 );
[h]=hist( r, 1:5 )
r = aleatorio_especial( [.5 .25 .2 .04 .01], 1e4 );
[h]=hist( r, 1:5 )
r = aleatorio_especial( [.5 .25 .2 .04 .01], 1e4 );
[h]=hist( r, 1:5 )
r = aleatorio_especial( [.5 .25 .2 .04 .01], 1e4 );
[h]=hist( r, 1:5 )
r = aleatorio_especial( [.5 .25 .2 .04 .01], 1e4 );
[h]=hist( r, 1:5 )
r = aleatorio_especial( [.5 .25 .2 .04 .01], 1e4 );
[h]=hist( r, 1:5 )
r = aleatorio_especial( [.5 .25 .2 .04 .01], 1e4 );
[h]=hist( r, 1:5 )
sqrt(1e4*.01*(1-.01))
1/sqrt(2e3)
1/sqrt(1e3)
help switch
edit samambaia
tic, [ x, y, r ] = samambaia( 1e4 ); toc
tic, [ x, y, r ] = samambaia( 1e5 ); toc
plota_regra( x, y, r )
edit plota_regra
plota_regra( x, y, r, 'brgk' )
plota_regra( x, y, r, 'brkg' )
tic, [ x, y, r ] = samambaia( 1e5 ); toc
plota_regra( x, y, r, 'brkg' )
tic, [ x, y, r ] = samambaia( 1e5, [.01 .17 .17 .65] ); toc
plota_regra( x, y, r, 'brkg' )
tic, [ x, y, r ] = samambaia( 1e5 ); toc
plota_regra( x, y, r, 'brkg' )
