function y = vol_to_lev(x)
V = 1000*x; % volume yang masuk ke silo [cm3]
hk = 13; % tinggi kerucut terpotong [cm]
r = 2.5;  % jari-jari kecil [cm]
R = 13; % jari-jari besar [cm]
ho = (r/(R-r))*hk; % tinggi kerucut kecil [cm]
Vk = ((1/3)*(hk+ho)*pi*R^2) -((1/3)*ho*pi*r^2); % volume potongan kerucut [cm]
if V <= 0
    y = 0;
elseif V <= Vk
   y = ((3*(ho^2)*(V + ((1/3)*ho*pi*r^2))/(pi*r^2))^(1/3)) - ho; % konversi V ke y untuk V <= Vk;
else
   y = ((V-Vk)/(pi*R^2)) + hk; % konversi V ke y, untuk V> Vk;        
end