Esempi in Python per nuclei

Esempio 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()

Esempio n. 2

File: few_steps.py Progetto: 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()

Esempio n. 3

File: plot_rk.py Progetto: 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()    

Esempio n. 4

File: plot_emm.py Progetto: 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()