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

n = [1:1:100]';

Error_1=.99639478.^(n-1);
Error_2=.95.^(n-1);
Transpath=0.2*exp(Error_1-Error_2);


g = vertcat(0.2 , Transpath, 0.4 );
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.2-Transpath, 0.4);