clc
clear
c=8;
h=18;
o=0;
airs=12.5;
temp=(100:25:3000)';
press=(20:5:120)';
feq=(0.8:0.02:1.3)';
n5nt=zeros(49,4,26);
for i=1:49
    for j=1:21
        for b=1:26
            [n5nt(i,j,b),e(i,j,b)]=mydissfuniii(c,h,o,temp(i,1),feq(b,1),airs,press(j,1));
        end
    end
end
%%