jacobian for Bratu
0001 function [c,fu,flam,b]=bratujac(p,u,lam) 0002 % jacobian for Bratu 0003 u=pdeintrp(p.points,p.tria,u); c=1; 0004 fu=-10*(1-lam*exp(u)); flam=10*exp(u); b=0;