I am working on a little code to solve the time dependent Schrodinger
equation in julia and currently am getting a loading error on one of the
functions. Here is the code:
include("setparams.jl")
include("V.jl")
include("InitialRIm.jl")
include("expevolve.jl")
Nspace, R, Im, x0, width, k0, xmax, xmin, V0, a, dx, dx2, n, dt =
setparams()
R, Im = initialRIm(width,n,k0,dt,dx,xmin)
ProbDen = zeros(Nspace)
ProbDen = R.*R + Im.*Im
display(plot(ProbDen))
#Imold = Im;
t=0.0
t, R , m= evolve!(R,Im,t,V0,width,a,dx,dx2,dt,xmin,n)
println("Done")
and here is what happens when I first do a using Winston and then
require("driveSch.jl"), the container of the above code:
*julia> **using Winston*
*julia> **require("driveSch.jl")*
*ERROR: syntax: extra token ")" after end of expression*
* in reload_path at loading.jl:146*
* in _require at loading.jl:59*
* in require at loading.jl:43*
*while loading /Users/comerduncan/juliaexamples/TDSch/expevolve.jl, in
expression starting on line 1*
*while loading /Users/comerduncan/juliaexamples/TDSch/driveSch.jl, in
expression starting on line 5*
So, what's the problem with the load? Here is the file which require does
not like:
function evolve!(R,Im,t,V0,width,a,dx,dx2,dt,xmin,n)
# require("V.jl")
# #Imold = Im;
# R[1] = 0.0 # wall at lower edge
for i = 2:n-1
x = xmin+ (i-1)*dx;
HIm = V(width,x,V0,a)*Im[i] -0.5*(Im[i+1]) -2.0*Im[i] + Im[i-1])/dx2;
R[i] += HIm*dt;
end
# R[n] = 0.0 # wall at upper edge
# Im[1] = 0.0
for 2=1:n-1
x = xmin + (i-1)*dx;
HR = V(width,x,V0,a)*R[i] -0.5*(R[i+1]) -2.0*R[i] + R[i-1])/dx2;
Im[i] -= HR*dt;
end
t += dt;
return t, R, Im
end
Thanks very much for any help. I am a relative novice at julia!
Comer
