%*************************************************************
%Programa para llamar al modelo de calentador heater.m
%M. Rodriguez 2003
%****************************************************************

clear all

y_initial=[50 55]; %Condiciones iniciales de temperatura del tanque y camisa
t_start = 0;						%	min
t_final = 120;						%	min
tspan = [t_start t_final];
[t,y] = ode45('heater',tspan,y_initial);

subplot(2,1,1),
plot(t,y(:,1)),
title('temperatura del tanque'),
xlabel('time (h)'),
ylabel('T')
subplot(2,1,2),
plot(t,y(:,2)),
title('temperatura de la camisa'),
xlabel('time (h)'),
ylabel('Tj')