Hello, I am trying to solve the system of differential equations with the "method of lines" method. I have an array in the shape of
QU[t_] = Table[qu[i][t], {i, 0, n}];
I want certain points to have a value of 50 at a certain time, but the remaining points must be zero in this time interval. For example, for i = 1,2,3,4, let p [i] [t] = 50 for 0 <t <10 and p [i] [t] = 0 for the remaining points. I also want all the points to be 0 for t> 10. How to do that. Thank you