// Este es un primer archivo para introducirse en el uso del Scilab
// Lo que aparece a la derecha de la doble barra inclinada // y que se muestra en color verde es un comentario y no se toma en cuenta para hacer los cálculos.
// Si se quiere que al ejecutar no aparezcan en la pantalla las líneas completas debe ponerse al final punto y coma ;
// No obstante en la ejecución de un script mediante  el botón de play el comando por defecto es 
// -->exec('C:\Users\....\ejem 1.1.sce', -1)
// el -1 que aparece al final indica que no se va a mostrar en la consola ninguna línea del programa al ejecutarlo.

// Supongamos un tanque aislado donde se mezcla en forma perfecta una corriente caliente a 60ºC con otra fría a 10ºC.
// Predecir la temperatura de la corriente de salida en función de la fracción de la corriente caliente (asumiendo estado estacionario).

// Si llamamos  x  a dicha fracción   Ts = x*Tcal + (1-x)*Tfr

// Podemos definir una función  Ts = f(x) de la siguiente manera

function y = Ts(x)
Tcal = 60 //ºC , corriente caliente
Tfr  = 10 //ºC , corriente fría
y = x*Tcal + (1-x)*Tfr //ºC , corriente de salida
endfunction

// por ejemplo si deseamos calcular Ts para x = 40%
T40 = Ts(0.4)
disp(T40,"Temperatura (ºC) para 40% de fracción caliente =  ") // La función disp muestra el resultado en pantalla

pause // puede resultar útil para interrumpir momentaneamente la ejecución, se sale escribiendo "resume" en la consola.

// Las funciones también se pueden definir en un archivo independiente grabado como función.sci que puede ser cargado mediante el comando
//  getf('función.sci')  y ejecutado en el ambiente en el que se está trabajando mediante el comando
//  exec('f',x)   donde f es el nombre de la función y x es un valor de la variable al que se le quiere aplicar la función.


// Con las flechas para arriba y para abajo se pueden repetir las últimas líneas sin necesidad de reescribirlas.

// para graficar la función debemos construir un vector de valores de  x  y el correspondiente valor de la función  f(x)

x = 0:.05:1

// esto quiere decir que el vector  x  va de 0 a 1 y sus elementos están espaciados 0.05.
// es lo mismo que escribir   x = [0 0.05 0.10 0.15 ...... 0.95 1]

f = Ts(x) //calcula el vector de Ts

// ahora vamos a graficar usando la función  plot

plot(x,f)
title('temperatura de salida')
xlabel('fracción de corriente caliente')
ylabel('T salida (ºC)')
