However after the solution (1xn array) is achieved, i perform some other calculation and then i want to feed the end result back into the PDE (in the next iteration in a for loop) as initial condition.
Define
icfun = @(x)interp1(X,U,x)
where X is the mesh vector and U is the vector of initial values corresponding to X you want to feed back into the solver for the next iteration.
Additinally how can additional paramters be passed into the pde function?
