Solving an array with NDSolve
- To: mathgroup at smc.vnet.net
- Subject: [mg89245] Solving an array with NDSolve
- From: svend <domdey at tphys.uni-heidelberg.de>
- Date: Sun, 1 Jun 2008 03:36:15 -0400 (EDT)
hello all,
i want to solve a PDE by discretizing in one variable and solving for the other one with NDSolve. unfortunately this second variable also appears in my discretization points and mathematica cannot handle this in a straightforward way.
let me post a simple example:
n = 10; x0 = 1; x1 = 10;
eq1 := {D[ff[j][x], x] == ff[problem][x], ff[j][x0] == j^2}
problem := IntegerPart[x]
eqsystem := Flatten[Table[eq1, {j, 1, n}]]
funs := Flatten[Table[ff[j], {j, 1, n}]];
test := NDSolve[eqsystem, funs, {x, x0, x1}][[1]];
i want to solve for ff[j][x] but mathematica cannot evaluate problem during the execution of NDSolve. Does anyone know how to fix this?
Many thanks
Svend