serie_II_2.mw

Serie 1 (Suter)

Funktion mass zur berechnung der Atommasse in MeV/c^2

> mass := proc(A,Z)

> # Konstanten

> # Umrechnung in MeV/c^2

> faktor := 5.6096E29; #/2.99792458E8; #???

> m_P:=1.67262158E-27 * faktor;

> m_N:=1.67492716E-27 * faktor;

> m_e:=9.10938188E-31 * faktor;

> m_H:=m_P+m_e; # nicht ganz, aber beinahe

> # a_x

> a_V:=15.85;

> a_S:=18.34;

> a_C:=0.71;

> a_A:=92.86;

> a_P:=11.46;

> # gg=+1, gu=0, uu=-1?

> rho:=0;

> if( type(Z,even) and type(A-Z,even) ) then

> rho:=+1;

> end if;

> if( type(Z,odd) and type(A-Z,odd) ) then

> rho:=-1;

> end if;

> # Energie berechnen (in MeV/c^2)

> Z*m_H + (A-Z) * m_N - a_V * A + a_S * A^(2/3) + a_C * Z^2 * A^(-1/3) + a_A * (Z-A/2)^2/A + rho * a_P * A^(-1/2);

> end proc;

Warning, `faktor` is implicitly declared local to procedure `mass`

Warning, `m_P` is implicitly declared local to procedure `mass`

Warning, `m_N` is implicitly declared local to procedure `mass`

Warning, `m_e` is implicitly declared local to procedure `mass`

Warning, `m_H` is implicitly declared local to procedure `mass`

Warning, `a_V` is implicitly declared local to procedure `mass`

Warning, `a_S` is implicitly declared local to procedure `mass`

Warning, `a_C` is implicitly declared local to procedure `mass`

Warning, `a_A` is implicitly declared local to procedure `mass`

Warning, `a_P` is implicitly declared local to procedure `mass`

Warning, `rho` is implicitly declared local to procedure `mass`

mass := proc (A, Z) local faktor, m_P, m_N, m_e, m_H, a_V, a_S, a_C, a_A, a_P, rho; faktor := 0.56096e30; m_P := 0.167262158e-26*faktor; m_N := 0.167492716e-26*faktor; m_e := 0.910938188e-30*faktor; m...mass := proc (A, Z) local faktor, m_P, m_N, m_e, m_H, a_V, a_S, a_C, a_A, a_P, rho; faktor := 0.56096e30; m_P := 0.167262158e-26*faktor; m_N := 0.167492716e-26*faktor; m_e := 0.910938188e-30*faktor; m...mass := proc (A, Z) local faktor, m_P, m_N, m_e, m_H, a_V, a_S, a_C, a_A, a_P, rho; faktor := 0.56096e30; m_P := 0.167262158e-26*faktor; m_N := 0.167492716e-26*faktor; m_e := 0.910938188e-30*faktor; m...mass := proc (A, Z) local faktor, m_P, m_N, m_e, m_H, a_V, a_S, a_C, a_A, a_P, rho; faktor := 0.56096e30; m_P := 0.167262158e-26*faktor; m_N := 0.167492716e-26*faktor; m_e := 0.910938188e-30*faktor; m...mass := proc (A, Z) local faktor, m_P, m_N, m_e, m_H, a_V, a_S, a_C, a_A, a_P, rho; faktor := 0.56096e30; m_P := 0.167262158e-26*faktor; m_N := 0.167492716e-26*faktor; m_e := 0.910938188e-30*faktor; m...mass := proc (A, Z) local faktor, m_P, m_N, m_e, m_H, a_V, a_S, a_C, a_A, a_P, rho; faktor := 0.56096e30; m_P := 0.167262158e-26*faktor; m_N := 0.167492716e-26*faktor; m_e := 0.910938188e-30*faktor; m...mass := proc (A, Z) local faktor, m_P, m_N, m_e, m_H, a_V, a_S, a_C, a_A, a_P, rho; faktor := 0.56096e30; m_P := 0.167262158e-26*faktor; m_N := 0.167492716e-26*faktor; m_e := 0.910938188e-30*faktor; m...mass := proc (A, Z) local faktor, m_P, m_N, m_e, m_H, a_V, a_S, a_C, a_A, a_P, rho; faktor := 0.56096e30; m_P := 0.167262158e-26*faktor; m_N := 0.167492716e-26*faktor; m_e := 0.910938188e-30*faktor; m...mass := proc (A, Z) local faktor, m_P, m_N, m_e, m_H, a_V, a_S, a_C, a_A, a_P, rho; faktor := 0.56096e30; m_P := 0.167262158e-26*faktor; m_N := 0.167492716e-26*faktor; m_e := 0.910938188e-30*faktor; m...mass := proc (A, Z) local faktor, m_P, m_N, m_e, m_H, a_V, a_S, a_C, a_A, a_P, rho; faktor := 0.56096e30; m_P := 0.167262158e-26*faktor; m_N := 0.167492716e-26*faktor; m_e := 0.910938188e-30*faktor; m...mass := proc (A, Z) local faktor, m_P, m_N, m_e, m_H, a_V, a_S, a_C, a_A, a_P, rho; faktor := 0.56096e30; m_P := 0.167262158e-26*faktor; m_N := 0.167492716e-26*faktor; m_e := 0.910938188e-30*faktor; m...mass := proc (A, Z) local faktor, m_P, m_N, m_e, m_H, a_V, a_S, a_C, a_A, a_P, rho; faktor := 0.56096e30; m_P := 0.167262158e-26*faktor; m_N := 0.167492716e-26*faktor; m_e := 0.910938188e-30*faktor; m...mass := proc (A, Z) local faktor, m_P, m_N, m_e, m_H, a_V, a_S, a_C, a_A, a_P, rho; faktor := 0.56096e30; m_P := 0.167262158e-26*faktor; m_N := 0.167492716e-26*faktor; m_e := 0.910938188e-30*faktor; m...mass := proc (A, Z) local faktor, m_P, m_N, m_e, m_H, a_V, a_S, a_C, a_A, a_P, rho; faktor := 0.56096e30; m_P := 0.167262158e-26*faktor; m_N := 0.167492716e-26*faktor; m_e := 0.910938188e-30*faktor; m...mass := proc (A, Z) local faktor, m_P, m_N, m_e, m_H, a_V, a_S, a_C, a_A, a_P, rho; faktor := 0.56096e30; m_P := 0.167262158e-26*faktor; m_N := 0.167492716e-26*faktor; m_e := 0.910938188e-30*faktor; m...mass := proc (A, Z) local faktor, m_P, m_N, m_e, m_H, a_V, a_S, a_C, a_A, a_P, rho; faktor := 0.56096e30; m_P := 0.167262158e-26*faktor; m_N := 0.167492716e-26*faktor; m_e := 0.910938188e-30*faktor; m...mass := proc (A, Z) local faktor, m_P, m_N, m_e, m_H, a_V, a_S, a_C, a_A, a_P, rho; faktor := 0.56096e30; m_P := 0.167262158e-26*faktor; m_N := 0.167492716e-26*faktor; m_e := 0.910938188e-30*faktor; m...

> # Testen (Werte in MeV)

> evalf(mass(194,78));

180674.5187

> evalf(mass(1,0));

965.2721397

> # Die Ergebnisse stimmen mit der Nuklid-Tabelle überein jedoch sind die Werte in MeV???

> evalf( mass(235,92) - ( mass(95,38) + mass(140,54) ) );

170.4535305

> # Vergleich mit Tabelle (Werte in MeV)

> ( 235.0439231 - ( 94.9193582 + 139.9216357 ) ) * 9.3149E2;

189.0265205

> # Auch nicht schlecht

> # Weiter Zerfall (Werte in MeV)

> # Sr-95->(Beta-) Y-95 ->(Beta-) Zr-95 ->(Beta-) Nb-95 ->(Beta-) Mo-95 (stable)  94.9058415 amu

> # Xe-140 ->(Beta-) Cs-140 ->(Beta-) Ba-140 ->(Beta-) La-140 ->(Beta-) Ce-140 (stable) 139.9054340 amu

> U_Zerfall := ( 235.0439231 - ( 94.9058415 + 139.9054340 ) ) * 9.3149E2;

U_Zerfall := 216.7089129

> # AKW

> # in einem Gramm angereichertes Uran sind wieviele U-235 Atome?

> N_A:=6.02214199E23;

N_A := 0.602214199e24

> U_Energie_MeV := 1 * 0.03 * N_A * U_Zerfall;

U_Energie_MeV := 0.3915155532e25

> # in Giga Joule

> U_Energie_GJJ := U_Energie_MeV * 1.6022E-13 * 1E-9;

U_Energie_GJJ := 627.2862193

> # Uran verbrauch in kg

> evalf((3 * Pi*1E7) / U_Energie_GJ)/1000;

150.2468518

>