Proleblem with empty plot

Hello enginners. I have problem with my mfile script. I want create a plot current against frequency, but when I create it plot is empty, Thank you for answers. Here is script:

h=[1:1:40]
%cinne zataze
R1=2;
R2=1.23;

Zns=((1.1*22e3)/(sqrt(3)*13.1e3))*(0.114+j*h*0.994);

Znsprevod=Zns*((0.4*0.4)/(22*22));

%trafo 22/0.4;

Ztr=j*h*(((400*400)/(1000e3))*(6/100));
%ASM 45kw;
Zam1=(400*400)/(45e3);
Xam1=Zam1/(sqrt(1.01));
Ram1=0.1*Xam1;
Zam1k=Ram1+j*h*Xam1;

Zam2=(400*400)/(45e3);
Xam2=Zam2/(sqrt(1.01));
Ram2=0.1*Xam2;
Zam2k=Ram2+j*h*Xam2;

Zam3=(400*400)/(80e3);
Xam3=Zam3/(sqrt(1.01));
Ram3=0.1*Xam3;
Zam3k=Ram3+j*h*Xam3;

Zam4=(400*400)/(35e3);
Xam4=Zam4/(sqrt(1.01));
Ram4=0.1*Xam4;
Zam4k=Ram4+j*h*Xam4;

Zam5=(400*400)/(50e3);
Xam5=Zam5/(sqrt(1.01));
Ram5=0.1*Xam5;
Zam5k=Ram5+j*h*Xam5;

Zam6=(400*400)/(35e3);
Xam6=Zam6/(sqrt(1.01));
Ram6=0.1*Xam6;
Zam6k=Ram6+j*h*Xam6;

%kompenzacne zariadenie

Xc=9.642./h;

Zc=-j*Xc;

Zpar=(1./Zam1k)+(1./Zam2k)+(1./Zam3k)+(1./Zam4k)+(1./Zam5k)+(1./Zam6k)+(1./Zc)+(1./R1)+(1./R2);

Zpar1=1./Zpar;

IPCC=(72./h).*(Zpar1./(Zpar1+Ztr+Znsprevod));

IC=(IPCC.*(Znsprevod+Ztr))/Zc

ans=abs(IC);

plot(h*50,ans)

Answer 1

[Star Strider]

You forgot to vectorise the division in the ‘IC’ calculation:

IC=(IPCC.*(Znsprevod+Ztr))./Zc
                          ↑ ← VECTORISE HERE

Then the plot works.

[Tomas Schnierer]

Thank you budyy

[Star Strider]

As always, my pleasure.