1

最初のモジュールは次の 3 つです。

import numpy as np
import matplotlib as plt
import math
import beginvalues as bv


def func(lijst,t):
    x = lijst[0]
    y = lijst[1]
    z = lijst[2]

    G = 6,67e-14
    MBH = 9,94e38
    N = 1 + x**2 + y**2/bv.a**2 + z**2/bv.b**2
    NT = (1e-3 + x**2 + y**2 + z**2)**(3/2)

    a_x = -(220**2)*x/N - (G * MBH * x)/NT
    a_y = -(220**2)*y/(N*bv.a**2) - (G * MBH * y )/NT
    a_z = -(220**2)*z/(N*bv.b**2) - (G * MBH * z )/NT
    return [lijst[3],lijst[4],lijst[5],a_x,a_y,a_z]

私の2番目のもの:

import numpy as np


# Ingeven van de parameterwaarden a en b. Indien deze niet relevant zijn, wordt een melding gegeven. 
tekst = "Geef waarden in voor a en b, waarbij 1>a>b:"
print tekst
a= float(input("Geef een waarde in voor a: "))
b= float(input("Geef een waarde in voor b: "))


if b<(a/(1+a**2)**(0.5)) or b>a:
    print " Dit zijn geen goede waarden voor a en b. \n Het model is niet fysisch relevant. \n Gelieve andere waarden te kiezen."



# Ingeven van de beginvoorwaarden.
# We maken een array aan waarin de beginwaarden geplaatst worden.

begin_x = 1
begin_y = 1
begin_z = 1
v0_x=1
v0_y=1
v0_z=1
init_val = [begin_x,begin_y,begin_z,v0_x,v0_y,v0_z]

そして私の3番目のもの:

import bvgzwartgat as bvg
from scipy.integrate import odeint
import numpy as np
import matplotlib.pyplot as plt
import beginvalues as bv


# We voeren een tijdsintegratie door van 0 tot 13*10^9 in 10000 stappen.
# De beginvoorwaarden worden geinitialiseerd in de module "beginvalues".
# TIJDSEENHEID?????????????????? kpc/km/s = 978 10^6 jaar
t=np.linspace(0,15,10000)
I=odeint(bvg.func,bv.init_val,t)




# 2d plotjes 
X=[]
Y=[]
Z=[]
for i in range(len(I)):
    kolom = I[i]
    X.append(kolom[0])
    Y.append(kolom[1])
    Z.append(kolom[2])



plt.figure("x-y projectie")
plt.plot(X,Y)

plt.figure("x-z projectie")
plt.plot(X,Z)

plt.figure("y-z projectie")
plt.plot(Y,Z)
plt.show( )

3番目のメインモジュールからプログラムを実行しようとすると、次のエラーが発生します: in func a_x = -(220**2)*x/N - (G * MBH * x)/NT TypeError: can'tシーケンスに 'tuple' 型の非 int を掛ける

だから私は、プログラムが a_x がタプルであると考えていると思います。

このコードは、以前は a_x の (G*MBH*x)/NT 部分がなくても完全に機能していました。したがって、プログラムが a_x をシーケンスと見なす理由を理解する必要があると思います。

4

1 に答える 1