%function lanza(FIG,I)
%
%
% Script de lanzamiento de temperatura
% Descomentar la primera linea para convertirlo en función
% 
% 
%
%

% ---------------------------------------------------------------------%
% Datos experimentales enviados por Ratis
%
% Parametros F'\isicos Dependen de FIG n
%
if FIG==0
    Q       = [5.89e7, 5.89e7    ];
    LAMBDA  = [13      ,13       ];
    A       = [4.23e-6,4.23e-6   ];
    
    V_DELTA = [0.53 0.53         ];
    DELTA   = [2.663e-3,2.663e-3 ];
    
    % Tiempos
    TP      = [1.272e-3,1.272e-3 ];
    TC      = [1.522e-3,1.522e-3 ]; 
    TS      = [0      , 1.79e-4  ];
    ALPHA_S = [0      , 27.29e4  ];
    ALPHA   = [5.207e4, 5.207e4  ];  

elseif FIG == 1
    % -----------   FIG 1 ------------------------------- 
    Q       = [28.52e7  25.92e7     21.94e7 ];
    
    LAMBDA  = [28.5     28.5        28.5 ];
    A       = [7.29e-6  7.29e-6     7.29e-6];
    
    V_DELTA = [0.53 0.53 0.53];
    DELTA   = [0.001204 0.001204 0.001204];
    
    % Tiempos
    TP      = [3.10e-3, 1.104e-3, 1.104e-3 ];
    TC      = [3.14e-3, 1.57e-3 , 1.57e-3]; % cambiado
    TS      = [0      , 3.79e-4 , 5.3e-4  ];
    ALPHA_S = [0      , 27.29e4 , 27.29e4];
    ALPHA   = [5.207e4, 5.207e4 , 5.207e4];
    
elseif FIG == 2    
    % -----------   FIG 2 ------------------------------- 
    Q       = [13.79e7  11.05e7 ];
    LAMBDA  = [29.9     29.9    ];
    A       = [6.39e-6  6.39e-6];
    
    V_DELTA = [0.45     0.45];
    DELTA   = [2.916e-3 2.916e-3];
    
    % Tiempos
    TP      = [1.046e-2 1.426e-3];
    TC      = [1.046e-2 1.745e-3];
    TS      = [0        2.83e-4];
    ALPHA_S = [0        27.29e4];
    ALPHA   = [5.207e4  5.207e4];
elseif FIG == 3
     % -----------   FIG 3 ------------------------------- 
    Q       = [8.62e7   8.48e7];
    LAMBDA  = [8.8      8.8   ];
    A       = [3.63e-6  3.63e-6];
    
    V_DELTA = [0.53     0.53   ];
    DELTA   = [2.66e-3  2.66e-3];
    
    % Tiempos
    TP      = [9.132e-3 1.272e-3];
    TC      = [9.132e-3 1.522e-3];
    TS      = [0        2.5e-4];
    ALPHA_S = [0        27.29e4];
    ALPHA   = [5.207e4  5.207e4];   
else
    error('check FIG')
end
% nt = 150;
% m  = 50;
% M  = 70;
global q lambda a v_delta delta factor tc tp ts alpha alpha_s

flag = 0;
nexp = 50;
i=2; 
q = Q(i);
lambda = LAMBDA(i);
a = A(i);
v_delta = V_DELTA(i);
delta =  DELTA(i);

KK=q/v_delta;

tc = TC(i);
tp = TP(i);
ts = TS(i);

%NT = 400;
%MM = [30 80 100 500];% 1000:1000:10000 100000];
%MMM = 50:10:200;
x =linspace(-delta,1*delta,50);
%x=0;
y=0;
m = 50;
ind = 0;
V = [0.1 0.3 0.53 0.7];
salida = cell(3,2);
nt =100
%load Barrido_completo_velocidades_x_0.mat -mat
%V = [(v_delta+0.01):.01:2];
for v_delta = V
	AAAA = clock;
   q=KK*v_delta;
   factor = q/2/lambda*a*sqrt(1/pi);
   [sn,T0] = tempy0(.035,x,nt,100);
   ind = ind +1;
   salida{1,ind} = sn;
   salida{2,ind} = T0;
   archivo = ['figura8.cuaprote.mat'];
   save (archivo,'-mat');      
   B = etime(clock,AAAA);
	fprintf('Tiempo total transcurrido: %8.6f seg \n',B);
end