Hi Irfan,
I understand that you want to iterate over columns of a matrix, checking each column's values to see any value of x is less than 1. And that once a value in the current column reaches 1, the iteration should move to the next column and repeat the process until all columns have been processed.
The below code outlines how you can structure the required logic:
% Assuming n and other utilized functions are pre-defined
itn = 1;
for i = 1:n-1
x = 0; % Initialize x for each column
while any(x < 1) % Check if any value of x is less than 1
itn = itn + 1;
% Update Rhn, qhn, and Vh based on provided code
Rhn(1:n, i) = ((pi * rh(1:n, i).^4) / (8 * nuD * L)) .* (1 ./ (x + (1 - x) * M));
qhn(1:n, i) = qh(1:n, i) .* Rh(1:n, i) ./ Rhn(1:n, i);
Vh(1:n, i) = qhn(1:n, i) ./ Ah(1:n, i);
% Determine delXmax based on the value of x
if any(x <= 0.5)
delXmax = 0.1;
else
delXmax = 0.07;
end
% Calculate deltat and update x
deltat = min(min(delXmax ./ Vh(1:n, i)));
x = x + Vh(1:n, i) * deltat;
if any(x >= 1)
break; %move to the next column
end
end
end
Hope this helps!
