%?sub?module?for?diffusion?step?in?simpltrans
for?i?=?2:N-1
????c2（i）?=?c1（i）?+?Neumann*（c1（i-1）-2*c1（i）+c1（