a = [2;3;4;5;2;2;2;2;2;2;2;2;2];
g0 = [2;2;2;2;3;4;5;2;2;2;2;2;2];
L = 400*ones(size(a));
L(8) = 300;
L(9) = 200;
L(10) = 100;

W= 400*ones(size(a));
W(11) = 300;
W(12) = 200;
W(13) = 100;

a = 1e-6*a;
g0 = 1e-6*g0;
L = 1e-6*L;
W= 1e-6*W;

da = -.28*1e-6;
dg = 0.49*1e-6;
sigmav = 0.5;
dv = sigmav*rand(size(a));
v = zeros(size(a));
for i=1:13
   v(i) = pullin(a(i)+da, g0(i)+dg, L(i), W(i), 1.68e11);
end

[a,g0,L,W,v+dv]
