From my code (below) I have K1 and K2 set to zero in order to test this numerical method. However, K1 is on an interval from 0 to 32, and K2 is on an interval from 0 to 8. I would like to know if there is a way to set these intervals to the associated variables, so that it will run the program for the interval itself. That way I don't have to run a combination of 32 or 8 separate calculations, etc.
Here is my code pasted directly from the nb.
Subscript[x, 0] := 0; Subscript[y, 0] := 1;
Subscript[u, 3] := 0; h := 0.1;
Subscript[K, 1] := 0; Subscript[K, 2] := 0
x[n_] := (Subscript[x, n] = Subscript[x, 0] + n h);
y[n_] := (Subscript[y, n ] = Subscript[y, 0] + n h);
Table[x[n], {n, 0, 5}];
Table[y[n], {n, 0, 5}]
f[n_] := (Subscript[f,
n] = -Subscript[K, 1] Subscript[x,
n] + -Subscript[K, 2] Subscript[y, n]);
Table[f[n], {n, 0, 5}]
System = {Subscript[u,
0] == {1/
80 (-27 h Subscript[f, 0] - 102 h Subscript[f, 1] -
72 h Subscript[f, 2] - 42 h Subscript[f, 3] +
3 h Subscript[f, 4] + 80 Subscript[u, 3])},
Subscript[u,
1] == {1/
90 (h Subscript[f, 0] - 34 h Subscript[f, 1] -
114 h Subscript[f, 2] - 34 h Subscript[f, 3] +
h Subscript[f, 4] + 90 Subscript[u, 3])},
Subscript[u,
2] == {1/
720 (-11 h Subscript[f, 0] + 74 h Subscript[f, 1] -
456 h Subscript[f, 2] - 346 h Subscript[f, 3] +
19 h Subscript[f, 4] + 720 Subscript[u, 3])},
Subscript[u,
4] == {1/
720 (-19 h Subscript[f, 0] + 106 h Subscript[f, 1] -
264 h Subscript[f, 2] + 646 h Subscript[f, 3] +
251 h Subscript[f, 4] + 720 Subscript[u, 3])},
Subscript[u,
5] == {1/
90 (29 h Subscript[f, 0] - 146 h Subscript[f, 1] +
294 h Subscript[f, 2] - 266 h Subscript[f, 3] +
269 h Subscript[f, 4] + 90 Subscript[u, 3])}};
variables = {Subscript[u, 0], Subscript[u, 1], Subscript[u, 2],
Subscript[u, 4], Subscript[u, 5]};
Solve[System, variables]
Subscript
. It will however copy and paste correctly. $\endgroup$