init-routine
0001 function p=ac6init(p) 0002 % init-routine 0003 p=stanparam(p); % set generic parameters to standard, if needed reset below.. 0004 p.neq=1;p.f=@acf;p.bcf=@ac6bcfx; p.jsw=3; 0005 pre=sprintf('%s',inputname(1));p=setfn(p,pre); 0006 lx=0.5;ly=0.5; p.geo=hexgeo(lx,ly); p=stanmesh(p,0.1); p=setbmesh(p); 0007 p.nsteps=50;p.ngen=10; p.lam=0; p.u=0*ones(p.np,1); p.ds=0.05; 0008 p.tau=zeros(p.neq*p.np+1,1); p.tau(p.neq*p.np+1)=1; p.xi=1/p.np; 0009 p.dlammax=0.2; p.lamdtol=0.8; p.lammin=-1; p.lammax=2; 0010