SciLab to Matlab transfer impedance code issue

I was provided he attached Scilab code for a continuing ed lesson and have gone back and forth with the instructor as to why my Matlab ode wont work to find the transfer impedance like his provided scilab code does:

k_1=2000;

k_2=5000;

k_3=1000;

m_1=2.5;

m_2=5.0;

f=linspace(1,25,100);

omega= 2*pi*f;

while i<20

R_m= input('R_m for both elastic elements =');

Z_1=R_m+i*(omega*m_1-(k_1+k_3)./omega);

Z_2=R_m+i*(omega*m_2-(k_2+k_3)./omega);

M=i*k_3./omega;

Z_12=(Z_1.*Z_2-M.^2)./M;

Level=20*log10(abs(Z_12));

plot(f,Level)

i=i+1;

An=input('Another?, Yes=1, No=0');

if An==0

abort

end

end

Answer 1

[Sean de Wolski]

abort

should be

break

[Phillip Rauscher]

Thank you for the quick reply- but it continues with the same error. If it is selected to try a new R_m value the chart continues to distort without resetting like the scilab code allows for.

[Sean de Wolski]

What's the error?