function[f] = fun_ex22(V)

a = 0.401; %Pa m^6 (coefficient)
b = 42.7*10^(-6); % m^3 (coefficient)

T = 300; %K (temperature)
p = 3.5*10^(7); % Pa (pressure)

N = 1000; % (molecule de CO2)
k = 1.3806503 * 10^(-23); %Joule K^(-1) (Boltzmann constant)

f = [p+a*(N/V)^2]*(V-N*b) - k*N*T;

end