Matlab - Ayuda Urgente

clc; clear all; close all;

disp('Ingrese el valor minimo de la razon entre frecuencias:');

a=input('');

disp('Ingrese el intervalo en que cambia:');

b=input('');

disp('Ingrese el valor maximo de la razon entre frecuencias:');

c=input('');

disp('Ingrese el valor minimo de la razon de amortiguamiento:')

g=input('');

disp('Ingrese el intervalo en que cambia:');

h=input('');

disp('Ingrese el valor maximo de la razon de amortiguamiento:');

i=input('');

for d=(a:b:c); e=(g:h:i);

    j=1./((((1-(d.^2)).^2)+(2*e*d.^2)).^(1/2));

end

f=(a:b:c);

xlim([a c]);

ylim([g i]);

plot(j);

hold on;

plot(f);

title('D vs B');

grid on;

xlabel('D');

ylabel('B');

clear all

beta=0.01:0.01:3;

s=0.01:0.01:1;

hold on

for i=1:length(s)

D=1./sqrt((1-beta.^2).^2+(2*s(i)*beta).^2);

plot(D,beta)

end

hold off