%% Entrada datos-Resolucion
clear all;
close all;
% Parametros
a=0;
b=1.22;
L=b-a;
h=1e-2;
N=L/h-1;
% Condiciones de contorno
uaT=0.382;
ubT=0;
uaD=0;
ubD=0.179;
uaM=0;
ubM=0.161;
uaF=0;
ubF=0.761;
% datos del flujo
delta=0.0013556; % eta velocidad conveccion
muT=2.66e-5; % viscosidad trigliceridos
gamma1=0.00168; % sigma absorbcion
muD=3.44e-5; % viscosidad digliceridos
gamma2=0.0422; % sigma absorbcion
muM=4.79e-5; % viscosidad monogliceridos
gamma3=0.0516; % sigma absorbcion
muF=4.94e-5; % viscosidad FAME
gamma4=0; % sigma absorbcion
% forzamiento trigliceridos
bvpfunT=inline('0.*x','x');
% resolucion trigliceridos
[xhT,uhT]= bvp(a, b, N, muT, delta, gamma1, bvpfunT, uaT, ubT)
% forzamiento digliceridos
bvpfunD=gamma1.*uhT(2:end-1);
% resolucion digliceridos
[xhD,uhD]= bvpmod(a, b, N, muD, delta, gamma2, bvpfunD, uaD, ubD)
% forzamiento monogliceridosgliceridos
bvpfunM=gamma2.*uhD(2:end-1);
% resolucion monogliceridos
[xhM,uhM]= bvpmod(a, b, N, muM, delta, gamma3, bvpfunM, uaM, ubM)
% forzamiento FAME
bvpfunF=gamma1.*uhT(2:end-1)+gamma2.*uhD(2:end-1)+gamma3.*uhM(2:end-1);
% resolucion FAME
[xhF,uhF]= bvpmod(a, b, N, muF, delta, gamma4, bvpfunF, uaF, ubF)
% Representacion grafica de los resultados
figure
subplot(2,2,1)
plot(xhT,uhT)
title('Concentracion Trigliceridos')
legend('Trigliceridos')
subplot(2,2,2)
plot(xhD,uhD)
title('Concentracion Digliiceridos')
legend('Digliceridos')
subplot(2,2,3)
plot(xhM,uhM)
title('Concentracion Monogliiceridos')
legend('Monogliceridos')
subplot(2,2,4)
plot(xhF,uhF)
title('Concentracion FAME')
legend('FAME')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Guardar resultados
print -dpng grafico.png            
print -dpdf grafico.pdf            
salida = [xhT uhT xhD uhD xhM uhM xhF uhF];
save -ascii resultados.txt salida