outfu for chemotax, put ||u-1||_L^1 first on branch
0001 function out=chembra(p,u,lam) 0002 % outfu for chemotax, put ||u-1||_L^1 first on branch 0003 n=triint(abs(1-u(1:p.np)),p.points,p.tria)/p.vol; 0004 out=[n; max(u(1:p.np))];