taui= repmat(0, 102, 1);
tauc= repmat(0, 102, 1);
tauk= repmat(0, 102, 1);

g = vertcat([ 0.2;0.2;0.2;0.2;0.2;0.2;0.2;0.2;0.2; 0.2 ], repmat(0.4, 92, 1));
k = vertcat([ 1.5;1.5;1.5;1.5;1.5;1.5;1.5;1.5;1.5; 1.5 ], repmat(1.5, 92, 1));
c = vertcat([ 0.6;0.6;0.6;0.6;0.6;0.6;0.6;0.6;0.6; 0.6 ], repmat(0.4, 92, 1));