Exemplos de nuclei em Python

Exemplo n.º 1

import calc as ca

# euler;rk;emm;exact

test = ca.nuclei(10000, 10, 0.001)
test.exact(100)

rs = []
h = []

for i in range(1000):
    step = float((i + 1)) / 1000.
    test.step = step
    test.exact(100)
    test.rk(100)
    tot = 0
    for j in range(len(test.nu)):
        tot += (test.base[j] - test.nu[j]) * (test.base[j] - test.nu[j])
    #print tot,",",len(test.nu)
    tot = tot / float(len(test.nu))
    rs.append(tot)
    h.append(step)

ca.plt.plot(h, rs)
ca.plt.show()

Exemplo n.º 2

Arquivo: few_steps.py Projeto: aragornranger/computationalphysics_N2013301020051

import calc as ca

test = ca.nuclei(1000, 10, 1)
test.exact(10)
test.euler(10)
ca.plt.plot(test.t, test.nu, label="euler", color="blue")
ca.plt.plot(test.t, test.base, label="exact", color="black")
test.rk(10)
ca.plt.plot(test.t, test.nu, label="r-k", color="red")
test.emm(10)
ca.plt.plot(test.t, test.nu, label="mid-point", color="green")

ca.plt.legend(loc="upper right")
ca.plt.show()

Exemplo n.º 3

Arquivo: plot_rk.py Projeto: aragornranger/computationalphysics_N2013301020051

import calc as ca

# euler;rk;emm;exact

test = ca.nuclei(10000,10,0.001)
test.exact(100)
    
rs = []
h = []

for i in range(1000):
    step = float((i + 1)) / 1000.
    test.step = step
    test.exact(100)
    test.rk(100)
    tot = 0
    for j in range(len(test.nu)):
        tot += (test.base[j] - test.nu[j]) * (test.base[j] - test.nu[j])
    #print tot,",",len(test.nu)
    tot = tot / float(len(test.nu))
    rs.append(tot)
    h.append(step)
    
        

ca.plt.plot(h,rs)
ca.plt.show()    

Exemplo n.º 4

Arquivo: plot_emm.py Projeto: aragornranger/computationalphysics_N2013301020051

import calc as ca

# euler;rk;emm;exact

test = ca.nuclei(10000,10,0.1)
test.exact(100)
test.euler(100)

rs = []
h = []

for i in range(1000):
    step = float((i + 1)) / 1000.
    test.step = step
    test.exact(100)
    test.emm(100)
    tot = 0
    for j in range(len(test.nu)):
        tot += (test.base[j] - test.nu[j]) * (test.base[j] - test.nu[j])
    #print tot,",",len(test.nu)
    tot = tot / float(len(test.nu))
    rs.append(tot)
    h.append(step)
    
        

ca.plt.plot(h,rs)
ca.plt.show()