if
for i=1:length(indices)
[Max_mdot_corr(i),k(i)]=max(mdot_corr(indices{1,i}));
[Min_mdot_corr(i),k1(i)]=min(mdot_corr(indices{1,i}));
end
ind=length(k1)
for g=1:ind
indices_speed(g)=indices{1,g}(k1(g))
end
for q= 1:ind
PI_C_tt_speed{q}=(PI_C_tt(indices_speed(q)))
end
for s=1:length(k)
zeilenum_max(s)=indices{1,s}(k(1,s));
PI_C_tt_comp_max(s)=PI_C_tt(zeilenum_max(1,s));
zeilenum_min(s)=indices{1,s}(k1(1,s));
PI_C_tt_comp_min(s)=PI_C_tt(zeilenum_min(1,s));
end
mdot_corr_vec=mdot_corr(:);
PI_C_tt_vec=PI_C_tt(:)
eta_c_tt_is_vec=eta_c_tt_is(:);
hold on
for j=1:length(indices)
plot3(mdot_corr(indices{1,j}(1:end)),PI_C_tt(indices{1,j}(1:end)),eta_c_tt_is_hun(indices{1,j}(1:end)),'-bs','Linewidth',1.5,...
'MarkerEdgeColor','b');
plot(Max_mdot_corr(1,:),PI_C_tt_comp_max(1,:),'linewidth',2);
plot(Min_mdot_corr(1,:),PI_C_tt_comp_min(1,:),'linewidth',2);
Wnodes=40;
PRnodes=40;
smooth=1.5;
[eff_grid,M_grid,P_grid]=gridfit(mdot_corr_vec,PI_C_tt_vec,eta_c_tt_is_vec,Wnodes,PRnodes,'smoothness',smooth);
tx=linspace(min(mdot_corr),max(mdot_corr),40);
ty=linspace(min(PI_C_tt),max(PI_C_tt),40);
eff_grid_text=ceil(eff_grid*100);