# Fisica para Ciencias Biologicas _ marco 2009
# NC
# este programa corre em octave, linux debian.
# para correr em Windows: va para o prompt de octave e digite
# edit logistico2009.m
# para que voce possa ler o que tem dentro (este arquivo)
# volte ao prompt do octave e digite
# logistico2009


# Programa para fazer graficos gerados por iteracao do mapa logistico
# O que estiver aa direita do simbolo do jogo da velha eh ignorado 
# pelo computador, mas serve como comentario para o programador.

# estas linhas em branco foram ignoradas

# antes de comecar limpe tudo

clear

#INPUT

# defina o numero de iteracoes maximo
t=100
# as linhas que nao terminam em ponto e virgula sao ecoadas na tela
# por isso aparece o valor de t.

#declare a variavel x: abra um espaco na memoria.
x = ones (1, t);

# defina a variavel mu, que aqui tera o nome a

a=3.9

# defina o valor inicial x(1)

x(1)=.7;

# AQUI CORRE O PROGRAMA
# comeca a iteracao 
for i = 2:t;
  x(i)=a*x(i-1)*(1.-x(i-1));
endfor;

# acabou a iteracao

# precisamos escrever algo para ver o que o computador fez
#OUTPUT
#defina a grade de pontos para fazer o grafico
k = linspace (1, t, t);
# ha muitas formas de plotar os resultados

#retire os comentarios duplos  (jogo  da velha duas vezes)
# das linhas abaixo para ver o resultado dos diferentes tipos
# de graficos 

# forma mais simples
##figure(1)
##plot(k,x)

# nao queremos perder os dados, entao abrimos a figura dois
##figure(2)
#onde plotamos usando simbolos nao unidos por uma linha
##plot(k,x,"@12")

#mas queremos simbolos e linhas
##figure(3)
##plot(k,x,"-@12")

# e queremos indicar o que significam os eixos
#figure(4)
xlabel("tempo, valor de mu = 3.7")
ylabel("x")
title("Mapa Logistico")
plot(k,x,"-@12")


# mas voce quer gerar um arquivo que guarde para sempre
# esta informacao


# no formato encapsulated postcript
##print -deps serietemporalLogistico_mu_002_x37.eps

# ou no formato png
print -dpng serietemporalLogistico_mu_002_x37.png

# mude o valor de a, ou o valor inicial, isto eh, x(1)
# e gere outros graficos. Nao esqueca de mudar
#o  nome do grafico gerado de forma que o anterior nao seja apagado
# use nomes informativos.
