jacobian for Bratu, text expression syntax, just a template
0001 function [c,fu,flam,b]=bratujact(p,u,lam) 0002 % jacobian for Bratu, text expression syntax, just a template 0003 lam=mat2str(lam); c=1; fu=char(['-10*(1-' lam '.*exp(u))']); b=0; 0004 flam=char('10*exp(u)'); [c,fu,flam]=pdetxpd(p.points,p.tria,u,c,fu,flam);