Dear all, I'm learning how to use implicit TS solvers to calculate the flux. The example code is http://www.mcs.anl.gov/petsc/petsc-current/src/ts/examples/tutorials/ex10.c.html I'm dying in understanding Lines 232-255. Is there anyone who knows how to work this out? Thank you in advance and sorry for bothering others. Best regards, Feng-