please do it on matlab

Exercise 4: Using the function you create for the Peng Robinson equation of state, what is the pressure (in Pa) of refrigerant R-134a at T = 380 K and densities of 3 mol/ mand 7,000 mol/m?? Create a function to represent the Peng Robinson equation below where the inputs into the function are the temperature (K), molar volume (m /mol), and the output is pressure (Pa) RT a(T) P= V-b V (V +b)+(V-b) where R?T? a(T)=0.45724 RT a(T), b=0.07780 . P K=0.37464 +1.54220-0.2699202 o is the Accentric Factor a(n)=[^{ - eq 5 You function should start something like: function f = peng_rob (T,V); Hints The units of R, P, and V must be consistent with each other. Stick with SI units: mol/m, K, and Pa. The critical constants and the accentric factor for R-134a are: Te = 374.21 K, Pc = 4.06x106 Pa, @= 0.327. R = 8.3145 J/mol.K Check your equation at ideal gas conditions and against data from NIST. You have been given densities in the question but specific volumes are the input! If you get an error like "Subscript indices must either be real positive integers or logicals." then you have forgotten to put a multiplication symbol in somewhere. For example: b (v-2) + b* (v-2). Exercise 4: Using the function you create for the Peng Robinson equation of state, what is the pressure (in Pa) of refrigerant R-134a at T = 380 K and densities of 3 mol/ mand 7,000 mol/m?? Create a function to represent the Peng Robinson equation below where the inputs into the function are the temperature (K), molar volume (m /mol), and the output is pressure (Pa) RT a(T) P= V-b V (V +b)+(V-b) where R?T? a(T)=0.45724 RT a(T), b=0.07780 . P K=0.37464 +1.54220-0.2699202 o is the Accentric Factor a(n)=[^{ - eq 5 You function should start something like: function f = peng_rob (T,V); Hints The units of R, P, and V must be consistent with each other. Stick with SI units: mol/m, K, and Pa. The critical constants and the accentric factor for R-134a are: Te = 374.21 K, Pc = 4.06x106 Pa, @= 0.327. R = 8.3145 J/mol.K Check your equation at ideal gas conditions and against data from NIST. You have been given densities in the question but specific volumes are the input! If you get an error like "Subscript indices must either be real positive integers or logicals." then you have forgotten to put a multiplication symbol in somewhere. For example: b (v-2) + b* (v-2)