Unfortunately, a "for" loop is not ideal because it has a fixed number of iterations. If you want to change "h" dynamically, it's better to use a "while" loop, which lets you control how "h" changes and when the loop ends.
Refer to the following code implementation using the "while" loop:
% Parameters
om = 0.4;
f = 0.92e6;
del = 1600;
Ik = 0.952e6;
alfa = 0.00092;
c0 = 2e-3;
kap = 2.175;
kaplam = -1.35;
kapmy = 0.175;
h = 0.0001;
z0 = 0.0;
p0 = 0.1e6;
t0 = 0.0;
%initial conditions
z(1) = z0;
p(1) = p0;
t(1) = t0;
i = 1;
while t(i) < 0.026 % equivalent to 260 steps of h=0.0001
z(i+1) = z(i) + h*(p(i)/Ik);
psi(i+1) = kap*z(i+1) + (kaplam*z(i+1)^2) + (kapmy*z(i+1)^3);
p(i+1) = (f*om*psi(i+1)) / (c0 - om/del*(1 - psi(i+1)) - alfa*om*psi(i+1));
t(i+1) = t(i) + h;
% Check if z(i+1) or psi(i+1) > 1, reduce step
if z(i+1) > 1 || psi(i+1) > 1
h = h / 2; % reduce time step
end
i = i + 1;
end
disp('Max pressure is')
disp(max(p))
plot(t, p, 'r', 'LineWidth', 3);
xlabel('Time [s]', 'FontSize', 15);
ylabel('Pressure [Pa]', 'FontSize', 15);
grid on;
I hope this resolves your query!
Thanks
