acinit

 init-routine

0001 function p=acinit(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.jac=@acjac;
0005 pre=sprintf('%s',inputname(1)); p=setfn(p,pre); % set filenames
0006 lx=1; ly=0.9; [p.geo, bc]=recdbc1(lx,ly,1e3); p.bcf=@(p,u,lam) bc;
0007 nx=20;ny=nx;p=stanmesh(p,nx,ny);p=setbmesh(p); % geometry, BC, and mesh
0008 p.u=0*ones(p.np,1); p.lam=1; p.ds=0.1;p.nsteps=30; %
0009 p.tau=zeros(p.neq*p.np+1,1); p.tau(p.neq*p.np+1)=1; p.xi=1/p.np; 
0010 p.amod=20; p.dlammax=0.2; p.lamdtol=0.8; 
0011