function y = armstrong(a,x)
% Armstrong and Frederick model for the center of the yield surface
% Used for Example 9.7
% a(1) = C
% a(2) = gamma
y = a(1).*(1-exp(-a(2).*x));