É concedida permissão para copiar, distribuir e/ou modificar este documento nos termos da Licença de Documentação Livre GNU, versão 1.2 ou qualquer versão posterior publicada pela Free Software Foundation; sem Secções Invariantes, sem textos de Capa e sem textos de Contra-Capa. É incluída uma cópia da licença na secção intitulada GNU Free Documentation License.http://www.gnu.org/copyleft/fdl.htmlGFDLGNU Free Documentation Licensetruetrue

	A utilização deste ficheiro é regulada nos termos da licença Creative Commons - Atribuição-CompartilhaIgual 3.0 Não Adaptada.
	Pode: partilhar – copiar, distribuir e transmitir a obra recombinar – criar obras derivadas De acordo com as seguintes condições: atribuição – Tem de fazer a devida atribuição da autoria, fornecer uma hiperligação para a licença e indicar se foram feitas alterações. Pode fazê-lo de qualquer forma razoável, mas não de forma a sugerir que o licenciador o apoia ou subscreve o seu uso da obra. partilha nos termos da mesma licença – Se remisturar, transformar ou ampliar o conteúdo, tem de distribuir as suas contribuições com a mesma licença ou uma licença compatível com a original.
	Esta marca de licenciamento foi adicionada a este ficheiro durante a atualização da licença GFDL.http://creativecommons.org/licenses/by-sa/3.0/CC BY-SA 3.0Creative Commons Attribution-Share Alike 3.0truetrue

R-Quelltext:

#Orenstein-Uhlenbeck Prozess mit normalverteiltem Anfangswert
ou<-function(theta=1, mu=1.2, sigma=0.3, init=c(mu, sigma^2/(2*theta)), from=0, to=2, steps=500*(to-from)) {
  t<-seq(from,to,length=steps)
  oup<-rnorm(n=1,mean=init[1],sd=sqrt(init[2]))
  for (i in (1:(steps-1))) {
    oup[i+1]=oup[i] + theta*(mu-oup[i])*(to-from)/steps+sigma*rnorm(n=1,mean=0,sd=sqrt(sigma^2*(to-from)/steps))
  }
  return(oup)
}

cols=c("navy", "darkolivegreen4", "red")
t<-seq(0,2,length=1000)
set.seed(311051)

png(filename="OrnsteinUhlenbeck3.png", width=1200, height=900, pointsize=12)
par(bg="lightgrey")
matplot(t,ou(init=c(0,0)),type="l", xlab="", ylab="",col=cols[1], lwd=2, ylim=c(0,2))
lines(t,ou(init=c(2,0)),col=cols[2], lwd=2)
lines(t,ou(),col=cols[3], lwd=2)
title(main="Ornstein - Uhlenbeck", cex.main=2.5)
text(x=1.85,y=0,expression(list(theta==1, mu==1.2, sigma == 0.3)) )
text(x=1.55,y=1.75,expression(paste("d",X[t]==theta(mu-X[t]),"dt+",sigma,"d", W[t])),cex=2.5)
text(x=0,y=0.07, expression(X[0]==0),col=cols[1])
text(x=0,y=1.92, expression(X[0]==2),col=cols[2])
text(x=0.03,y=1.29, expression(paste(X[0], "~N(", mu,",", frac(sigma^2,2*theta), ")") ),col=cols[3])
abline(h=1.2, lty=3, col="grey25")
dev.off()

 Este(a) PNG gráfico foi criado com o R.