Gaussian


In [11]:
from __future__ import division
import numpy as np, scipy as sp, seaborn as sns, pandas as pd
from matplotlib import pyplot as plt
%matplotlib inline
sns.set_style("darkgrid")

In [12]:
def gaussian(x,A,x0,sigma):
    return A*np.exp(-(x-x0)**2/(2*sigma**2))

In [13]:
x = np.linspace(-3,3,100)

In [14]:
plt.plot(x,gaussian(x,A=1,x0=1,sigma=1))
plt.xlabel('$x$')
plt.ylabel('$P(x)$')


Out[14]:
<matplotlib.text.Text at 0x10a8fdc90>

In [15]:
from IPython.html.widgets import interact

In [10]:
def plotgaussian(A,x0,sigma):
    x = np.linspace(-3,3,1000)
    plt.plot(x,gaussian(x,A=A,x0=x0,sigma=sigma),label='${a}*\exp(-(x-{x0})^2/(2*{sigma}^2))$'.format(a=A,x0=x0,sigma=sigma))
    plt.legend()
    plt.xlabel('$x$')
    plt.ylabel('$P(x)$')
    plt.title('gaussian, $A$={a},$x_0$={x0},$\sigma$={sigma}'.format(a=A,x0=x0,sigma=sigma))
interact(plotgaussian,A=(.1,10),
         x0=(-2,2),
         sigma=(0.1,1),
         )



In [32]:
%%script perl
print $_ . "$_" for (0 ..1000);


0011223344556677889910101111121213131414151516161717181819192020212122222323242425252626272728282929303031313232333334343535363637373838393940404141424243434444454546464747484849495050515152525353545455555656575758585959606061616262636364646565666667676868696970707171727273737474757576767777787879798080818182828383848485858686878788888989909091919292939394949595969697979898999910010010110110210210310310410410510510610610710710810810910911011011111111211211311311411411511511611611711711811811911912012012112112212212312312412412512512612612712712812812912913013013113113213213313313413413513513613613713713813813913914014014114114214214314314414414514514614614714714814814914915015015115115215215315315415415515515615615715715815815915916016016116116216216316316416416516516616616716716816816916917017017117117217217317317417417517517617617717717817817917918018018118118218218318318418418518518618618718718818818918919019019119119219219319319419419519519619619719719819819919920020020120120220220320320420420520520620620720720820820920921021021121121221221321321421421521521621621721721821821921922022022122122222222322322422422522522622622722722822822922923023023123123223223323323423423523523623623723723823823923924024024124124224224324324424424524524624624724724824824924925025025125125225225325325425425525525625625725725825825925926026026126126226226326326426426526526626626726726826826926927027027127127227227327327427427527527627627727727827827927928028028128128228228328328428428528528628628728728828828928929029029129129229229329329429429529529629629729729829829929930030030130130230230330330430430530530630630730730830830930931031031131131231231331331431431531531631631731731831831931932032032132132232232332332432432532532632632732732832832932933033033133133233233333333433433533533633633733733833833933934034034134134234234334334434434534534634634734734834834934935035035135135235235335335435435535535635635735735835835935936036036136136236236336336436436536536636636736736836836936937037037137137237237337337437437537537637637737737837837937938038038138138238238338338438438538538638638738738838838938939039039139139239239339339439439539539639639739739839839939940040040140140240240340340440440540540640640740740840840940941041041141141241241341341441441541541641641741741841841941942042042142142242242342342442442542542642642742742842842942943043043143143243243343343443443543543643643743743843843943944044044144144244244344344444444544544644644744744844844944945045045145145245245345345445445545545645645745745845845945946046046146146246246346346446446546546646646746746846846946947047047147147247247347347447447547547647647747747847847947948048048148148248248348348448448548548648648748748848848948949049049149149249249349349449449549549649649749749849849949950050050150150250250350350450450550550650650750750850850950951051051151151251251351351451451551551651651751751851851951952052052152152252252352352452452552552652652752752852852952953053053153153253253353353453453553553653653753753853853953954054054154154254254354354454454554554654654754754854854954955055055155155255255355355455455555555655655755755855855955956056056156156256256356356456456556556656656756756856856956957057057157157257257357357457457557557657657757757857857957958058058158158258258358358458458558558658658758758858858958959059059159159259259359359459459559559659659759759859859959960060060160160260260360360460460560560660660760760860860960961061061161161261261361361461461561561661661761761861861961962062062162162262262362362462462562562662662762762862862962963063063163163263263363363463463563563663663763763863863963964064064164164264264364364464464564564664664764764864864964965065065165165265265365365465465565565665665765765865865965966066066166166266266366366466466566566666666766766866866966967067067167167267267367367467467567567667667767767867867967968068068168168268268368368468468568568668668768768868868968969069069169169269269369369469469569569669669769769869869969970070070170170270270370370470470570570670670770770870870970971071071171171271271371371471471571571671671771771871871971972072072172172272272372372472472572572672672772772872872972973073073173173273273373373473473573573673673773773873873973974074074174174274274374374474474574574674674774774874874974975075075175175275275375375475475575575675675775775875875975976076076176176276276376376476476576576676676776776876876976977077077177177277277377377477477577577677677777777877877977978078078178178278278378378478478578578678678778778878878978979079079179179279279379379479479579579679679779779879879979980080080180180280280380380480480580580680680780780880880980981081081181181281281381381481481581581681681781781881881981982082082182182282282382382482482582582682682782782882882982983083083183183283283383383483483583583683683783783883883983984084084184184284284384384484484584584684684784784884884984985085085185185285285385385485485585585685685785785885885985986086086186186286286386386486486586586686686786786886886986987087087187187287287387387487487587587687687787787887887987988088088188188288288388388488488588588688688788788888888988989089089189189289289389389489489589589689689789789889889989990090090190190290290390390490490590590690690790790890890990991091091191191291291391391491491591591691691791791891891991992092092192192292292392392492492592592692692792792892892992993093093193193293293393393493493593593693693793793893893993994094094194194294294394394494494594594694694794794894894994995095095195195295295395395495495595595695695795795895895995996096096196196296296396396496496596596696696796796896896996997097097197197297297397397497497597597697697797797897897997998098098198198298298398398498498598598698698798798898898998999099099199199299299399399499499599599699699799799899899999910001000

In [ ]: