Untitled

% Parametri
k = 0.005; %[W/m*K]
d = 0.3; %[m]
T0 = 253; %[K]
Qs_giorno = 1200; %[W/m^2]
Qs_notte = 0; %[W/m^2]
function y = temperatura(T, Qs, T0, k, d)
    % Funzione per calcolare l'equazione energetica
    Qr = 0.89 * 5.67e-8 * T^4; %[W/m^2]
    Qc = (T - T0) * k / d; %[W/m^2]
    y = Qs - Qr - Qc; % Equazione da risolvere
end
% Risoluzione per il giorno
T0_giorno = 300; % Stima iniziale
T_giorno = fsolve(@(T) temperatura(T, Qs_giorno, T0, k, d), T0_giorno);
% Risoluzione per la notte
T0_notte = 100; % Stima iniziale
T_notte = fsolve(@(T) temperatura(T, Qs_notte, T0, k, d), T0_notte);
% Risultati
fprintf('T_giorno= %.4f K\n', T_giorno);
fprintf('T_notte= %.4f K\n', T_notte);