clc;
clear;

% Input parameters
c = 8;
h = 18;
o = 0;
airs = 12.5;
temp = (100:25:3000)';
press = (20:5:120)';
feq = (0.8:0.02:1.3)';

% Initialize result storage
n5nt = zeros(49, 4, 26);
iterations = zeros(49, 4, 26);
computationTimes = zeros(49, 4, 26);

% Loop through temperature, pressure, and fuel equivalence ratio
for i = 1:length(temp)
    for j = 1:length(press)
        for b = 1:length(feq)
            [n5nt(i, j, b), iterations(i, j, b), computationTimes(i, j, b)] = ...
                Improved_mydissfuniii(c, h, o, temp(i), feq(b), airs, press(j));
        end
    end
end

% Display average performance metrics
avgIterations = mean(iterations(:));
avgComputationTime = mean(computationTimes(:));

fprintf('Average Iterations: %.2f\n', avgIterations);
fprintf('Average Computation Time: %.4f seconds\n', avgComputationTime);

% Save results to a .mat file for later analysis
save('results.mat', 'n5nt', 'iterations', 'computationTimes');