Lyapunov exponent of the logistic map (Mathematica Code)

In a previous post I'd shown a way to get the Lyapunov exponent from the time series data of any map. In this quick tutorial, I'll show you a cleaner way to get the Lyapunov exponent for the specific case of the logistic map, and then using a really short script in Mathematica, plot it …

Continue reading Lyapunov exponent of the logistic map (Mathematica Code)

Lyapunov Exponent of Logistic Map

Calculating the Lyapunov Exponent of a Time Series (with python code)

(In a later post I discuss a cleaner way to calculate the Lyapunov exponent for maps and particularly the logistic map, along with Mathematica code.) I found this method during my Masters while recreating the results of an interesting paper on how some standard tests for chaos fail to distinguish chaos from stochasticity (Stochastic neural network …

Continue reading Calculating the Lyapunov Exponent of a Time Series (with python code)

Target Distribution in Gnuplot

R code for multivariate random-walk Metropolis sampling

I couldn't find a simple R code for random-walk Metropolis sampling (the symmetric proposal version of Metropolis Hastings sampling) from a multivariate target distribution in arbitrary dimensions, so I wrote one. This is also my first R code. It requires the package MASS to sample from the multivariate normal proposal distribution using the mvrnorm function. …

Continue reading R code for multivariate random-walk Metropolis sampling

Javascript Slideshow Code

I prefer writing all my web-design code from scratch. So I use only text editors to write HTML directly, and raw JavaScript (no jQuery etc). While working on a javascript slideshow for Artarium, I came across a lot of problems in transitioning the images. My javascript was changing the src for the image, but I needed to …

Continue reading Javascript Slideshow Code