Untitled
%Plottare la funzione f(V) R=0.08206; A0=16.6037; B0=0.2354; a=0.11171; b=0.07697; c=3*1e6; T=408; P=36; beta=R*T*B0-A0-(R*c)/T^2; delta=(R*B0*b*c)/T^2; gamma=-R*T*B0*b+a*A0-(R*c*B0)/T^2; fun=@(V) ((R*T + beta/V + gamma/V^2 + delta/V^3)*1/P -V); V0=2; sol = fsolve(fun, V0); fprintf('V= %.4f L\n', sol); V_vals = linspace(1, 5, 100); % Genera 100 valori di V tra 1 e 5 F_vals = arrayfun(fun, V_vals); % Calcola i valori della funzione per ogni V figure; plot(V_vals, F_vals, 'b-', 'LineWidth', 1.5); hold on; xlabel('Volume (L)'); ylabel('F(V) '); title('Tank Temperature'); legend('Variazione volume'); grid on;
Leave a Comment