% Este es un primer archivo para introducirse en el uso del Octave
% 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 
% >>'C:\Users\....\ejem_1_1.m')

% Es recomendable utilizar los siguientes tres comandos de Octave
clear all; % Borra todas las variables de la memoria ( workspace)
close all; % Cierra todas las ventas de gráfico
clc; % Limpia la consola

% 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("Temperatura (ºC) para 40% de fraccion caliente =  "); disp(T40); % La función disp muestra el resultado en pantalla

pause % puede resultar útil para interrumpir momentaneamente la ejecución, para seguir la ejecución se debe apretar cualquier botón del teclado.

% Las funciones también se pueden definir en un archivo independiente grabado como función.m que puede ser ejecutado 
% mediante la llamada al nombre del archivo donde está guardada la función, por ejemplo 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)')