versioninfo()

Pkg.init()

Pkg.add("Unitful")

#import Pkg; 
Pkg.status()

Pkg.clone("https://github.com/laguer/PhysicalConstant.jl.git")

Pkg.add("PhysicalConstant")

using Unitful
using Measurements
using MathPhysicalConstants.SI

big(SI.PlanckConstantHbar)/(big(SI.MassElectron)*big(SI.SpeedOfLight))


big(SI.SanchezElectrConstant)

measurement(SI.PlanckConstantHbar)

big(measurement(SI.PlanckConstantHbar))

using Unitful
using Measurements
using MathPhysicalConstants.SI

a1 = big(e^e^e) 

a1^e

a = SI.SanchezElectrConstant

using Unitful
using Measurements
using Cosmology
using PyCall
@pyimport numpy as np

np.log(137)/ np.log(a/137)

using Unitful
using Measurements
using Cosmology
using PyCall
@pyimport numpy as np
@pyimport sympy as sy
@pyimport decimal

#np.exp(2^e^2.5)
np.exp(2^e)

exp(2^(e^2.5))

using Unitful
using Measurements
using Cosmology
using PyCall
@pyimport numpy as np
@pyimport sympy as sy
@pyimport decimal

# decimal.Decimal(np.exp(2^e^2.5))
decimal.Decimal(np.exp(e^(2*e)+e^2))

using Unitful
using Measurements
using Cosmology
using PyCall
@pyimport numpy as np
@pyimport sympy as sy
@pyimport scipy as sci
@pyimport matplotlib.pyplot as plt
@pyimport decimal
println(np.e^np.e^np.e)


(SI.PlanckConstantH*SI.SpeedOfLight)/(SI.Boltzmann * 1 * 2.8977729e-3)



using Unitful
using Measurements
using Cosmology
using PyCall
@pyimport math
println(math.pi)


i = e
exp(i*pi)

Sum=0.0
for n = 1:2:10000
    Sum += 1/n - 1/(n+1)
end
Sum # an expression on the last line (if it doesn't end with ";") is printed as "Out"

Out[26] - log(2)

println(math.e^128-1)


using PhysConsts
α.quantity


big((α.value)^-1)


big(log((α.value)^-1))


big(math.exp(699))


import Pkg; Pkg.status()

using PhysConsts.DATA

c.quantity

a0.value

G.value

g.quantity

R = rand(50,50)
R' * R

println("Hello world!\n")
println(STDERR, "Testing, some unicode output on stderr...\n")

ccall(:printf, Cint, (Ptr{UInt8},), "Hello from C!!\n");

using Unitful
using Measurements
using Cosmology
using PyCall

@pyimport math
@pyimport numpy as np
f(x) = np.power(np.log(2), (x / 4 ))

using Unitful
using Measurements
using Cosmology
using PyCall
@pyimport math
@pyimport numpy as np
println(f(2))
f([1.1,1.2,2,3,5,8])

using RDatasets
sleep = dataset("lme4","sleepstudy")

using Gadfly
plot(sleep, x = "Days", y = "Reaction", Geom.point)

plot(sleep, x = "Reaction", Geom.density, color = "Subject", Scale.x_continuous(minvalue= 0, maxvalue= 500))

plot(x = sleep[:Reaction], Geom.histogram(bincount = 30), Scale.x_continuous(minvalue = 200), color = sleep[:Days])

display("text/html", """Hello <b>world</b> in <font color="red">HTML</font>!""")

type HTML
   s::String
end
import Base.show
show(io::IO, MIME, x::HTML) = print(io, x.s)

x = HTML("<ul> <li> Hello from a bulleted list! </ul>")

display(x)
#println(x)