


pde for AC with lam=mu=Diff coeff.,
set p.up1=old_lambda ("User Parameter 1") somewhere

0001 function [c,a,f,b]=acfmu(p,u,lam) 0002 % pde for AC with lam=mu=Diff coeff., 0003 % set p.up1=old_lambda ("User Parameter 1") somewhere 0004 u=pdeintrp(p.points,p.tria,u); c=lam; a=0; f=p.up1*u+u.^3-u.^5; b=0;