Hi Robert,

Thank you very much. But I do not understand what you mean by "Note that you still need to specify some boundary conditions so that the matrix is not singular". Because it seems no boundary condition is constrained in the modified codes, but it still runs, and the result looks correct. Would you please explain more about your saying?
 
Regards
Ronghai


在 2014年11月15日星期六UTC+1下午7时02分31秒,Robert Cimrman写道:
On 11/14/2014 06:02 PM, Ronghai Wu wrote:
> Hi Robert,
>
> I have some explanations of fipy and my idea. Please see attached PDF. If
> you need more information, please just ask.

That's fine, thanks. I have modified your file and it runs. Note that you still
need to specify some boundary conditions so that the matrix is not singular.

r.