Hi everyone,
I have a system of PDE's.
PDE1 = D[u0[x, t], t] - Exp[-x] == 0
PDE2 = D[u1[x, t], t] + D[u0[x, t], x] -
2*D[u0[x, t], {x, 2}, {t, 1}] + Exp[-x] == 0
PDE3 = D[u2[x, t], t] + D[u1[x, t], x] -
2*D[u1[x, t], {x, 2}, {t, 1}] == 0
PDE4 = D[u3[x, t], t] + D[u2[x, t], x] -
2*D[u2[x, t], {x, 2}, {t, 1}] == 0
PDE5 = D[u4[x, t], t] + D[u3[x, t], x] -
2*D[u3[x, t], {x, 2}, {t, 1}] == 0
PDE6 = D[u5[x, t], t] + D[u4[x, t], x] -
2*D[u4[x, t], {x, 2}, {t, 1}] == 0
ics = {u5[x, 0] == 0, u4[x, 0] == 0, u3[x, 0] == 0, u2[x, 0] == 0,