clc;
clear;

% Given constants
c = 8;
h = 18;
o = 0;
airs = 12.5;
temp = (100:25:3000)';  % Temperature range
press = (20:5:120)';    % Pressure range
feq = (0.8:0.02:1.3)';  % Fuel equivalence ratio

% Initialize results matrix
n5nt = zeros(49, 21, 26);
error_matrix = zeros(49, 21, 26);

% Compute using fzero method
for i = 1:49  % Temperature index
    for j = 1:21  % Pressure index
        for b = 1:26  % Fuel equivalence ratio index
            [n5nt(i, j, b), error_matrix(i, j, b)] = mydissfun_fzero(c, h, o, temp(i), feq(b), airs, press(j)); % Using fzero function
        end
    end
end

% Display the first few values
disp('Computed n5nt values using fzero Method:');
disp(n5nt(:,:,1));  % Display n5nt for first equivalence ratio