function [c,a,f,b]=bratuf(p,u,lam) 
% pde for Bratu 
u=pdeintrp(p.points,p.tria,u); c=1; a=0;f=-10*(u-lam*exp(u)); b=0;