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