Write a computer code in Matlab to generate a phase diagram for $n-$dodecane. Your code should do the following:

Use the Peng$-$Robinson $($PR$)$ equation of state $($EOS$)$ to compute the PVT properties of the component assigned.

Implement a cubic root$-$solving algorithm to return the liquid and vapor roots $($if applicable$).$

Generate a phase diagram for $n-$dodecane at the following set of temperatures: $,0\mathrm{.}8,0\mathrm{.}9,1\mathrm{.}0,1\mathrm{.}1,1\mathrm{.}25,2\mathrm{.}0.$