t=linspace(0,1,25)
for k=1:length(t)
A=(a1^2*P)/(i*mu*L*alpha^2)
B1=i^1.5*alpha*r/a1
B2=i^1.5*alpha
B11=besselj(0,B1)
B22=besselj(0,B2)
B=1-(B11/B22)
C=exp(i*w*t(k))
u=real(A*B*C)
plot(u,r,u,-r)
drawnow
end
используйте drawnow
после вашей команды заговора. это заставит matlab немедленно нарисовать график, что позволит вам увидеть ход симуляции.
EDIT:
Я не видел последний бит. Чтобы установить ось x & y в положение от 0 до 1, выполните xlim([0 1])
и ylim([0 1])