---

% --- Задание 1 ---

% 1. Задаем диапазон температур (от 1 до 21 °C, как на графике)
t = 1:0.5:21; 

% 2. Задаем давления в МПа (как на графике)
p_MPa = [5, 10, 15, 20, 25]; 

% 3. Переводим МПа в кгс/см^2, так как этого требует формула
p_kgf = p_MPa * 10.197; 

% 4. Считаем растворимость по эмпирической формуле
% Формула: phi = 0.4342 * p^0.64 / (t + 5.278)^0.3882
phi = zeros(length(p_kgf), length(t)); % Создаем пустую матрицу
for i = 1:length(p_kgf)
    phi(i, :) = 0.4342 * p_kgf(i)^0.64 ./ (t + 5.278).^0.3882;
end

% 5. Строим номограмму
figure;
plot(t, phi, 'LineWidth', 2);
grid on;
xlabel('Температура, °C');
ylabel('Содержание растворенного газа, см^3/г');
title('Растворимость метана в воде');
legend('5 МПа', '10 МПа', '15 МПа', '20 МПа', '25 МПа', 'Location', 'best');

% 6. Выводим значения в командное окно для точек, показанных на графике
t_points = [1, 5, 9, 13, 17, 21];
fprintf('\nЗначения растворимости (см^3/г) в точках:\n');
fprintf('T(C)\t5 МПа\t10 МПа\t15 МПа\t20 МПа\t25 МПа\n');
for T = t_points
    vals = 0.4342 * p_kgf.^0.64 ./ (T + 5.278).^0.3882;
    fprintf('%d\t%.2f\t%.2f\t%.2f\t%.2f\t%.2f\n', T, vals);
end