       difference equations

• To: mathgroup at smc.vnet.net
• Subject: [mg20567] difference equations
• From: psalzman at landau.ucdavis.edu
• Date: Sat, 30 Oct 1999 00:14:12 -0400
• Sender: owner-wri-mathgroup at wolfram.com

```dear all,

i wrote a difference equation at low order, and looked at the same thing at
higher order.  it's absurd.   so it's Mathematica or nothing since i have no grad
students to pawn off the work to (i am a grad student).

here's the basic starting point:

(* a = dt/(2 hbar)  and  b = -hbar^2/(2 m)   *)
(1 + I a * H) Psi[j + 1, i] == (1 - I a * H) Psi[j, i]

i would like
H -> b(1/r[i]^2) deriv_r(r[i]^2 deriv_r) + V[j + 1, i]
on the left, but
H -> b(1/r[i]^2) deriv_r(r[i]^2 deriv_r) + V[j, i]
on the right.   how can i do this?

also, my 'stubs' for the derivative, deriv_r shouldn't commute with
anything, and Mathematica likes to commute stuff.  is there a way to tell it that
i want deriv_r to remain on the left side of whatever it starts out to the
left of?

lastly, i'm using Psi[j, i], etc  to mean Psi at time j, grid i.   is there
a better way of doing this?

these are some basic questions, i know.  but up until now all i've ever used
Mathematica for is to solve integrals, plot graphs and spacetimes.

thanks!
pete

```

• Prev by Date: notebooks