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