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