Tekanan Uap Gas Propana dan Butana

Baru-baru ini, gue dapet tugas Fisika Material Komputasi dari Prof. Kamsul Abraha sebagai pengganti Ujian Mid Semester Ganjil 2011/2012. Tugasnya ntu disuruh menentukan tekanan uap gas propana dan gas butana, yang sudah diketahui tekanan kritis dan temperatur kritisnya. Cara menentukannya adalah dengan metode komputasi. Akhirnya gue pake program Matlab 2009 buat menentukan tekanan uapnya.

Matlab, kita harus membuat listing program dan me-run-kan program tersebut. Ini dia listing programnya:

clear; clc;

disp('*>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<*');
disp('| FISIKA MATERIAL KOMPUTASI |');
disp('| MENENTUKAN TEKANAN UAP GAS PROPANA DAN BUTANA |');
disp('| ARIEF GUNAWAN |');
disp('| 07/256754/PA/11785 |');
disp('*>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<*');
disp(' ');
%Variabel Masukan
Gas=input('Tuliskan nama gasnya : ','s');
Tc=input('Tuliskan suhu kritis gas dalam derajat celcius : ');
Pc=input('Tuliskan harga tekanan kritis gas dalam atm : ');
%T=input('Tuliskan suhu gas dalam Celcius : ');
n=input('Tuliskan jumlah mol gas dalam mol : ');
T=[0,20,50,70,100];
Tc=Tc+273;
T=T+273;
R=0.082056;
%mencari nilai Vc
a=27*(R^2)*(Tc^2)/(64*Pc);
b=R*Tc/(8*Pc);
Vc=3*b;
fPV=@(V) ((n*R .*T)./(V-n*b))-((a*n^2)./V^2);
Pm=fPV(Vc);
%output
disp(' ');
disp(['Untuk gas ',Gas,' diperoleh :']);
disp(['a : ',num2str(a),' atm.liter^2/mol^2']);
disp(['b : ',num2str(b),' liter/mol']);
disp(['Vc : ',num2str(Vc),' liter']);
%disp(['Pm : ',num2str(Pm),' atm']);
Hasil=[T-273;Pm];
Hasil=Hasil';
disp(' ');
disp(['Tabel Variasi Gas ',Gas]);
disp('-----------------------');
disp('Temperatur Tekanan Uap ');
disp('-----------------------');
disp(Hasil);
disp('>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<');


Setelah di-run, ini dia hasilnya.

Propana

Butana

Semoga bermanfaat. :)