package body destilacion_auxiliar is
   -- Ecuacion (2)
   function temperatura(st1, st2 : Float) return Float is
   begin
      return (st1 + st2)/2.00;
   end temperatura;

   -- Ecuación (1)
   function temperatura_campo(st1, sr1, tem, st4 : Float) return Float is
       beta:Float:=0.13;
       leq:Float:=15.00;
       H:Float:=4.00;
       c:Float:=1080000.00;
       Cp:Float:= 4190.00;
       p:Float:=975.00;
       sc1:Float:=20.00;
   begin
      return st1 + ((beta*leq*sr1) - (H*(tem-st4))) * (c/(sc1*Cp*p));
   end temperatura_campo;

   -- Ecuación (3)
   function sensor_entrada(st2 : Float) return Float is
   begin
       return st2 - 11.00;
   end sensor_entrada;

   -- Ecuación (4)
   function destilado(st2, st3 : Float) return Float is
       sc2:Float:=400.00;
   begin
       return 24.00*(0.135 + 0.003 * st2 - 0.0204 * st3 - 0.001 * sc2
              + 0.00004 * st2 * sc2);
   end destilado;

   -- Parte 3

   -- Despejar SC1 de la ecuación (1)
   function sensor_caudal_uno(sr1, st4, tem, st2, st1_anterior : Float) return Float is
   begin
      return (72.00*(39.00*sr1+80.00*(st4-tem)))/(5447.00*(st2-st1_anterior));
   end sensor_caudal_uno;

   -- Despejar SC2 de la ecuación (4)
   function sensor_caudal_dos(st2_anterior, st3, sd1 : Float) return Float is
   begin
      return (5.00*(306.00*st3+625.00*sd1-45.00*(st2_anterior+45.00)))/(3.00*(st2_anterior-25.00));
   end sensor_caudal_dos;
end destilacion_auxiliar;