jacobian for AC, 'text' syntax, only here as a template
0001 function [c, fu, flam,b]=acjact(p,u,lam) 0002 % jacobian for AC, 'text' syntax, only here as a template 0003 slam=mat2str(lam); c=0.2; fu=char([slam '+3*u.^2-5*u.^4']); flam=char('u'); 0004 [c,fu,flam]=pdetxpd(p.points,p.tria,u,c,fu,flam);b=0;