% timestep
time = zeros([],1) ; % variable to save time
t=ti;%intitial time step
count=1;
time(count) = t ;
while t < tf
for i=2:(nx-1)%should be 2D
for j=2:(ny-1)
Grid_order =(j-1)*nx +i;
k=Grid_order;
B(k,1)= alpha*Pcurrent(k,1);
B(Prod_Grid_order,1)=Pprod;
end
end
% Solution for Future time pressure
Pnplus1=A\B;
Prod_Grid_order=(6-1)*nx + 6;
Pnplus1(Prod_Grid_order,1)=Pprod;
Pcurrent=Pnplus1;
Psave(:,j)=Pcurrent; % store matrix
t=t+dt;
count = count+1 ;
time(count) = t ;
end
Or else you can create your own time step instead of saving, as the time step is dt always.
time = ti+(1:count)*dt ;
