function p=vteamparam(ver)

if ver==1
  p.Ron=100;     p.Roff=1000;
  p.R=200;
  p.woff=10e-9;  p.won=0;
  p.aoff=3;      p.aon=3;
  p.kon=-8e-6;   p.koff=4e-6;
  p.voff=0.15;   p.von=-0.20;
  p.forceinside=0;
  p.winit=6e-9; %0
  p.Region=3; % if the initial condition wstart is changed then p.Region must be set properly
  p.Umax=0.35;  p.T=0.01;     p.omega=2*pi/p.T;
  p.u=@(t)(p.Umax*sin(p.omega*t));
end

if ver==2
  p.Ron=100;     p.Roff=1000;
  p.R=200;
  p.woff=10e-9;  p.won=0;
  p.aoff=3;      p.aon=3;
  p.kon=-8e-6;   p.koff=4e-6;
  p.voff=0.15;   p.von=-0.20;
  p.forceinside=1;
  p.winit=6e-9; %0
  p.Region=3; % if the initial condition wstart is changed then p.Region must be set properly
  p.Umax=0.35;  p.T=0.01;     p.omega=2*pi/p.T;
  p.u=@(t)(p.Umax*sin(p.omega*t));
end
