Diffusion pNbody (master) Browse

pNbody/examples/ic/box/plot1.py
a2d7a927d00cmaster

plot1.py
View Options

			#!/usr/bin/env python

			import Ptools as pt
			from numpy import *

			y = arange(0,1,0.001)

			rho1 = 1.
			rho2 = 2.
			sigma = 0.025
			y1 = 0.50

			rho = rho1 + (rho2-rho1)/( (1+exp(-2.*(y-y1)/sigma)) )


			pt.plot(y,rho)
			pt.show()