% task_3.m
clear; close all; clc;

% Define the constants and parameters
c = 8; 
h = 18; 
o = 0; 
airs = 12.5; 
press = 105; 

% Temperature and feq parameters
temp = 400;
feq = 0.98; 
% Call the function for solving
try
    [n5nt, e, iterations] = mydissfuniii_fzero(c, h, o, temp, feq, airs, press);
    disp(['Solution n5nt: ', num2str(n5nt)]);
    disp(['Energy e: ', num2str(e)]);
    disp(['Iterations: ', num2str(iterations)]);
catch ME
    disp('Error occurred:');
    disp(ME.message);
end