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;