Mathematica 9 is now available
Student Support Forum
-----
Student Support Forum: 'mathematica ndsolve - no output' topicStudent Support Forum > General > Archives > "mathematica ndsolve - no output"

< Previous Comment | Next Comment >Help | Reply To Comment | Reply To Topic
Author Comment/Response
D
02/21/13 11:52am

I don't believe the problem is in the code because the program runs fully through on another computer. I believe it has something to do with how my computer is analyzing the notebook when it comes to NDsolve.

xM = 15

x0=-8

a = 0.5;
barrierHeight = 5.5 10^-7;

barrierV[x_] := Which[x < 0, 0, x > a, 0, True, barrierHeight]

c = 3000;
hbarc = 0.00197;
hbar = hbarc/ c

mass = 1000;
hbarSqOver2mass = (hbarc^2)/ (2 mass)

E1 = 4.9 10^-7
kzero = Sqrt[E1/hbarSqOver2mass ]
lambda = 2 \[Pi]/ kzero

v = E1/ (hbar)
\[Sigma] = 1
\[Omega]1 = c hbarc kzero/ mass
\[Omega]2 = c hbarc/ mass

\[Psi]pkt[v_, \[Sigma]_, \[Omega]1_, \[Omega]2_][x_,
t_] = (\[Sigma]^2/\[Pi])^(1/4)
Exp[I (kzero x - v t)]/((\[Sigma]^2 + I \[Omega]2 t )^(1/2))
Exp[-((x - \[Omega]1 t)^2/(2 (\[Sigma]^2 + I \[Omega]2 t) ))];

\[Chi][x_] := \[Psi]pkt[v, \[Sigma], \[Omega]1, \[Omega]2][x - x0, 0]

enforceBCs[\[Psi]_, x_,
xM_] := \[Psi] - (((\[Psi] /. x -> -xM) - (\[Psi] /. x -> xM))/
2*(Cos[(x + xM)/(2 xM) Pi] + 1) + (\[Psi] /. x -> xM))

schroedingerEq =
I hbar D[\[CapitalPsi][x, t], {t, 1}] == -
hbarSqOver2mass D[\[CapitalPsi][x, t], {x, 2}] +
barrierV[x] \[CapitalPsi][x, t]

(nsol = NDSolve[{schroedingerEq, \[CapitalPsi][x, 0] ==
enforceBCs[\[Chi][x], x, xM], \[CapitalPsi][xM, t] ==
0, \[CapitalPsi][-xM, t] == 0}, \[CapitalPsi][x, t], {x, -xM,
xM}, {t, 0, 200}, AccuracyGoal -> 3,
PrecisionGoal -> 3]) // Timing

URL: ,

Subject (listing for 'mathematica ndsolve - no output')
Author Date Posted
mathematica ndsolve - no output D. 02/20/13 9:31pm
Re: mathematica ndsolve - no output Forum Modera... 02/21/13 11:02am
Re: mathematica ndsolve - no output D 02/21/13 11:52am
Re: mathematica ndsolve - no output yehuda 02/21/13 7:02pm
Re: Re: mathematica ndsolve - no output D 02/23/13 9:58pm
Re: Re: Re: mathematica ndsolve - no output yehuda 02/24/13 6:31pm
Re: Re: Re: Re: mathematica ndsolve - no output Forum Modera... 02/24/13 7:33pm
Re: Re: Re: mathematica ndsolve - no output Bill Simpson 02/24/13 7:13pm
Re: Re: Re: Re: mathematica ndsolve - no output yehuda 02/25/13 01:48am
Re: Re: Re: Re: Re: mathematica ndsolve - no ou... D 02/27/13 9:54pm
Re: Re: Re: Re: Re: Re: mathematica ndsolve - n... jf 02/28/13 10:53am
Re: mathematica ndsolve - no output yehuda 02/28/13 5:07pm
< Previous Comment | Next Comment >Help | Reply To Comment | Reply To Topic