chemft

 pde for chemotaxis model, text syntax, another template

0001 function [c,a,f,b]=chemft(p,u,lam)
0002 % pde for chemotaxis model, text syntax, another template
0003 lams=mat2str(lam);  a=0; b=0; D='0.25'; r='1.52';
0004 c1111=D; c1112='0'; c1121='0'; c1122=D; 
0005 c1211=['-' lams '*u(1,:)']; c1212='0'; c1221='0';c1222=['-' lams '*u(1,:)'];
0006 c2111='0'; c2112='0'; c2121='0'; c2122='0'; 
0007 c2211='1'; c2212='0'; c2221='0'; c2222='1'; 
0008 c=char({c1111 c1121 c1112 c1122 c2111 c2121 c2112 c2122 ... 
0009     c1211 c1221 c1212 c1222 c2211 c2221 c2212 c2222});  
0010 f1=[r '*u(1,:).*(1-u(1,:))']; f2='u(1,:)./(1+u(1,:))-u(2,:)'; f=char({f1 f2});
0011 [c,a,f]=pdetxpd(p.points,p.tria,u,c,a,f); % evaluate on triangles