# 3. přednáška - předvedení přímého výpočtu Lorenzovy křivky rm(list=ls()) # vyčistí pracovní prostor # adresář nazvaný data obsahující datový soubor tolary.txt # je umístěn v pracovním adresáři (./data.tolary.txt) data(tolary) # načte datový soubor ls() # seznam objektů v prac. prostoru summary(tolary) # informace o datové tabulce tolary attach(tolary) # zpřístupnění datové tabulky y # výpis hodnot y1,...,yn sort(y) # variační řada (Y = cumsum(sort(y))) # kumulativní součty variační řady (n = length(y)) # počet hodnot plot((0:n)/n,c(0,Y)/Y[n],type="l") # lomená čára = Lorenzův oblouk points((0:n)/n,c(0,Y)/Y[n],pch=3) # vyznačení spojovaných bodů abline(a=0,b=1) # úhlopříčka čtverce abline(h=0:1,v=0:1,lty=3) # jednotkový čtverec