function [TM]=EnergyPairParamArbitrary(pair)

% Retrieve parameters of Entalpy and Entropy for any given nucleotide pair
% Currently only works for DNA, need verification to ensure RNA pairs hold
% the same parameters.


% NOTE: Will crash if any other letter is presented. Needs error coding.


% Constants Original
%               Pairs    Entalpy ΔH Entropy ΔS
%                         kcal/mol   cal/molK
% EntalpyPairs={'AA',	'TT',	 -7.9,	 -22.2;
%               'AT',	'AT',	 -7.2,	 -20.4;
%               'TA',	'TA',	 -7.2,	 -21.3;
%               'CA',	'TG',	 -8.5,	 -22.7;
%               'GT',	'AC',	 -8.4,	 -22.4;
%               'CT',	'AG',	 -7.8,	 -21;
%               'GA',	'TC',	 -8.2,	 -22.2;
%               'CG',	'CG',	-10.6,	 -27.2;
%               'GC',	'GC',	 -9.8,	 -24.4;
%               'GG',	'CC',	 -8,     -19.9};
          
% Experimental Constants        
%               Pairs        New MT (arbitrary)       Original MT 
%                                Aprox diff 2 
EntalpyPairs={  'AA',           14;               % 61.9183 
                'TT',           6;                    % 61.9183 
                'AT',		 	8;               % 57.4697
                'TA',           7;               % 44.3483
                'CA',           13;               % 79.8799 
                'TG',           5;                    % 79.8799  
                'GT',           12;               % 80.1284 
                'AC',           4;                    % 80.1284 
                'CT',           11;               % 75.4178 
                'AG',           3;                    % 75.4178 
                'GA',           10;               % 74.6424
                'TC',           2;                    % 74.6424
                'CG',           15;               % 97.7740 
                'GC',           0;              % 107.0297 
                'GG',           9;              % 102.8378 
                'CC',           1};                  % 102.8378

          
% Compare and determine parameter row
PairFound=strcmp(EntalpyPairs,pair);
[Index,C]=find(PairFound);


% Retrieve data
% Entalpy=EntalpyPairs{Index,3};
% Entropy=EntalpyPairs{Index,4};
TM=EntalpyPairs{Index,2};