% INPUT:
% ------
% t    = current time = scalar;
% Y    = is the  temperature in Celsius deg. 
%
% OUTPUT:
% -------
% DY   = is the initial value problem : right side f(Y(t);t)
%
% This function is called by:   Main_SLS_distribute.m
%---------------------------------------------------------------------------------
function [DY] = Standard_Linear_Model_SLS(t, Y)
global E_0 E_infty tau
global t_TAB EPSILON_TAB dot_EPSILON_TAB

%-----------------------------------
% the right side f(Y(t);t):   Begin
%--------------------------
Sigma    = Y(1);
tim      = Y(2); % t := Y(2)
 
epsilon     = interp1(t_TAB, EPSILON_TAB, tim); 
dot_epsilon = interp1(t_TAB, dot_EPSILON_TAB, tim); 

f_right_side  = + E_0 * dot_epsilon + E_infty/tau * epsilon + ...
                - Sigma / tau;

DY(1)  =  f_right_side ;
DY(2)  =  1.; % 

DY = DY';
%-----------------------------------
% the right side f(Y(t);t):   End
%--------------------------