clear?all;
clc;

pfa1=1e-3;
pfa2=1e-6;
pfa3=1e-9;
pfa4=1e-12;
pfa5=1e-14;

pd1=0.5;
pd2=0.6;
pd3=0.7;
pd4=0.8;
pd5=0.9;

index=0;

for?np=1:1:10e2
????index=index+1;
????I1（index）=improv_fac（nppfa1pd1）;
????I2（index）=improv_fac（nppfa2pd2）;
????I3（index）=improv_fac（nppfa3pd3）;
????I4（index）=improv_fac（n