jacobian for AC
0001 function [c,fu,flam,b]=acjac(p,u,lam) 0002 % jacobian for AC 0003 u=pdeintrp(p.points,p.tria,u); c=0.25; fu=lam+3*u.^2-5*u.^4; flam=u; b=0;