pde for AC
0001 function [c,a,f,b]=acf(p,u,lam) 0002 % pde for AC 0003 u=pdeintrp(p.points,p.tria,u); c=0.25; a=0; b=0; f=lam*u+u.^3-u.^5;