function op=stress(ii,q,T);

// T is an array 1 X n
// there are n calls to simulate 2017, each with T(i) days
//
// ii are the input heads in a 1 X m array. 
// recall that in ii, the first and last entries are not
// changed, i.e., they are constant head.
// q is an array the same size as iii and q(j) gives the 
// infiltration in meters per day

[ig, n]=size(T);

op=ii;

for i=1:n
 deltah=simulate2017(ii,q,T(i));
 ii=ii+deltah; // update of heads after time T(i)
 op=[op;ii]; // storing it as rows of matrix op
end;

return;

