bd output for GP
0001 function out=gpbra(p,u,lam) 0002 % bd output for GP 0003 imax=max(u(p.np+1:2*p.np)); rmax=max(u(1:p.np)); 0004 out=[rmax; imax; imax/rmax; triint(u(1:p.np).^2,p.points,p.tria)]; 0005 end