%Simulasi Nilai dari P Load
clc
clear

ro=1.2;                         %kg/m^3
cd=0.64;                        %drag coefficient
A=2.99;                         %frontal area
m=2160:-40:1840;                %massa kendaraan 1
v=((70)*1000)/3600;             %input kecepatan awal kendaraan km/jam
vt=((120)*1000)/3600;           %input kecepatan akhir kendaraan km/jam
cr=0.015;                       %koefisien rolling
g=9.81;                         %percepatan gravitasi
t=3600;                         %waktu perubahan kecepatan
a=(vt-v)/t;                     %percepatan dalam m/s^2
teta=(30*pi)./180;              %sudut tanjakan
H=9.9;                          %diesel kWH/l
ef=95/100;                      %efficiency
s=70*(t/3600)+0.5*(50)*(1.^2); %jarak dalam km
Pload=(0.5*ro*cd*A*(v.^3)+m*(cr*g*v+0.5*a*v+g*v*(sin(teta))))./(3600);
Q=Pload/(H*ef);
y=2160-m;


subplot(2,1,2);
plot (y,Q,'r-o');
xlabel ('Vehicle Mass Reduction (Kg)');
ylabel ('Fuel Consumption (liter/h)');
title  ('Reduce Mass vs Fuel Consumption');
grid on

x=2160:-40:1840;
subplot(2,1,1);
plot (x,Q,'ko');
xlabel ('Mass Vehicle(Kg)');
ylabel ('Fuel Consumption (liter/h)');
title  ('Mass Vehicle vs Fuel Consumption');