function y = vol_to_lev(x)
V = 1000*x; % volume yang masuk ke silo [cm3]
hk = 12.5; % 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 <= Vk
   y = ((3*(ho^2)*(V + ((1/3)*ho*pi*r^2))/(pi*r^2))^(1/3)) - ho; % konversi V ke y;
else
   
         y = ((v1-Vk)/(pi*R^2))) + hk; % konversi v1 ke y, y1>hk
end