function [ma,mb,mc,Drab,Drbc,Drac,lrab,lrbc,lrac,Brab,Brbc,Brac,mina,maxa,minb,maxb,cont,atma,atmb,atmc] = params(a,b,c)

% Gets the parameters for any atom set or returns error message if no
% parameters exist.

% labels:
% H F Cl D I O
% 1 2 3  4 5 6

% one body terms

if (a==1)
    ma = 1;
    atma = imread('hydrogen.png');
elseif (a==2)
    ma = 19;
    atma = imread('fluorine.png');
elseif (a==3)
    ma = 35.457;
    atma = imread('chlorine.png');
elseif (a==4)
    ma = 2;
    atma = imread('deuterium.png');
elseif (a==5)
    ma = 126.91;
    atma = imread('iodine.png');
elseif (a==6)
    ma = 16;
    atma = imread('oxygen.png');
end

if (b==1)
    mb = 1;
    atmb = imread('hydrogen.png');
elseif (b==2)
    mb = 19;
    atmb = imread('fluorine.png');
elseif (b==3)
    mb = 35.457;
    atmb = imread('chlorine.png');
elseif (b==4)
    mb = 2;
    atmb = imread('deuterium.png');
elseif (b==5)
    mb = 126.91;
    atmb = imread('iodine.png');
elseif (b==6)
    mb = 16;
    atmb = imread('oxygen.png');
end

if (c==1)
    mc = 1;
    atmc = imread('hydrogen.png');
elseif (c==2)
    mc = 19;
    atmc = imread('fluorine.png');
elseif (c==3)
    mc = 35.457;
    atmc = imread('chlorine.png');
elseif (c==4)
    mc = 2;
    atmc = imread('deuterium.png');
elseif (c==5)
    mc = 126.91;
    atmc = imread('iodine.png');
elseif (c==6)
    mc = 16;
    atmc = imread('oxygen.png');
end

% two body terms

if (a == 1) && (b == 1) || (a == 1) && (b == 4) || (a == 4) && (b ==1) || (a == 4) && (b == 4)
    Drab = 104.0;
    lrab = 0.74;
    Brab = 1.99;
elseif (a == 1) && (b == 2) || (a == 2) && (b == 1) || (a == 4) && (b == 2) || (a == 2) && (b == 4)
    Drab = 134.0;
    lrab = 0.92;
    Brab = 2.27;
elseif (a == 1) && (b == 3) || (a == 3) && (b == 1) || (a == 4) && (b == 3) || (a == 3) && (b == 4)
    Drab = 106.5;
    lrab = 1.28;
    Brab = 1.85;
elseif (a == 2) && (b == 2)
    Drab = 36.0;
    lrab = 1.42;
    Brab = 1.60;
elseif (a == 2) && (b == 3) || (a == 3) && (b == 2)
    Drab = 498.0;
    lrab = 1.63;
    Brab = 2.12;
elseif (a == 3) && (b == 3)
    Drab = 58.0;
    lrab = 1.99;
    Brab = 2.03;
elseif (a == 1) && (b == 5) || (a == 5) && (b == 1) || (a == 4) && (b == 5) || (a == 5) && (b == 4)
    Drab = 74.0;
    lrab = 1.60;
    Brab = 1.75;
elseif (a == 2) && (b == 5) || (a == 5) && (b == 2)
    errordlg('Parameters for this atom combination do not exist','Error','modal')
elseif (a == 3) && (b == 5) || (a == 5) && (b == 3)
    errordlg('Parameters for this atom combination do not exist','Error','modal')
elseif (a == 5) && (b == 5)
    Drab = 36.0;
    lrab = 2.67;
    Brab = 1.85;
elseif (a == 1) && (b == 6) || (a == 6) && (b == 1) || (a == 4) && (b == 6) || (a == 6) && (b == 4)
    Drab = 102.0;
    lrab = 0.96;
    Brab = 2.26;
elseif (a == 2) && (b == 6) || (a == 6) && (b == 2)
    errordlg('Parameters for this atom combination do not exist','Error','modal')
elseif (a == 3) && (b == 6) || (a == 6) && (b == 3)
    errordlg('Parameters for this atom combination do not exist','Error','modal')
elseif (a == 5) && (b == 6) || (a == 6) && (b == 5)
    errordlg('Parameters for this atom combination do not exist','Error','modal')
elseif (a == 6) && (b == 6)
    Drab = 54.0;
    lrab = 1.32;
    Brab = 2.32;
end    
    
if (b == 1) && (c == 1) || (b == 1) && (c == 4) || (b == 4) && (c ==1) || (b == 4) && (c == 4)
    Drbc = 104.0;
    lrbc = 0.74;
    Brbc = 1.99;
elseif (b == 1) && (c == 2) || (b == 2) && (c == 1) || (b == 4) && (c == 2) || (b == 2) && (c == 4)
    Drbc = 134.0;
    lrbc = 0.92;
    Brbc = 2.27;
elseif (b == 1) && (c == 3) || (b == 3) && (c == 1) || (b == 4) && (c == 3) || (b == 3) && (c == 4)
    Drbc = 106.5;
    lrbc = 1.28;
    Brbc = 1.85;
elseif (b == 2) && (c == 2)
    Drbc = 36.0;
    lrbc = 1.42;
    Brbc = 1.60;
elseif (b == 2) && (c == 3) || (b == 3) && (c == 2)
    Drbc = 498.0;
    lrbc = 1.63;
    Brbc = 2.12;
elseif (b == 3) && (c == 3)
    Drbc = 58.0;
    lrbc = 1.99;
    Brbc = 2.03;
elseif (b == 1) && (c == 5) || (b == 5) && (c == 1) || (b == 4) && (c == 5) || (b == 5) && (c == 4)
    Drbc = 74.0;
    lrbc = 1.60;
    Brbc = 1.75;
elseif (b == 2) && (c == 5) || (b == 5) && (c == 2)
    errordlg('Parameters for this atom combination do not exist','Error','modal')
elseif (b == 3) && (c == 5) || (b == 5) && (c == 3)
    errordlg('Parameters for this atom combination do not exist','Error','modal')
elseif (b == 5) && (c == 5)
    Drbc = 36.0;
    lrbc = 2.67;
    Brbc = 1.85;
elseif (b == 1) && (c == 6) || (b == 6) && (c == 1) || (b == 4) && (c == 6) || (b == 6) && (c == 4)
    Drbc = 102.0;
    lrbc = 0.96;
    Brbc = 2.26;
elseif (b == 2) && (c == 6) || (b == 6) && (c == 2)
    errordlg('Parameters for this atom combination do not exist','Error','modal')
elseif (b == 3) && (c == 6) || (b == 6) && (c == 3)
    errordlg('Parameters for this atom combination do not exist','Error','modal')
elseif (b == 5) && (c == 6) || (b == 6) && (c == 5)
    errordlg('Parameters for this atom combination do not exist','Error','modal')
elseif (b == 6) && (c == 6)
    Drbc = 54.0;
    lrbc = 1.32;
    Brbc = 2.32;
end  

if (a == 1) && (c == 1) || (a == 1) && (c == 4) || (a == 4) && (c ==1) || (a == 4) && (c == 4)
    Drac = 104.0;
    lrac = 0.74;
    Brac = 1.99;
elseif (a == 1) && (c == 2) || (a == 2) && (c == 1) || (a == 4) && (c == 2) || (a == 2) && (c == 4)
    Drac = 134.0;
    lrac = 0.92;
    Brac = 2.27;
elseif (a == 1) && (c == 3) || (a == 3) && (c == 1) || (a == 4) && (c == 3) || (a == 3) && (c == 4)
    Drac = 106.5;
    lrac = 1.28;
    Brac = 1.85;
elseif (a == 2) && (c == 2)
    Drac = 36.0;
    lrac = 1.42;
    Brac = 1.6;
elseif (a == 2) && (c == 3) || (a == 3) && (c == 2)
    Drac = 498.0;
    lrac = 1.63;
    Brac = 2.12;
elseif (a == 3) && (c == 3)
    Drac = 58.0;
    lrac = 1.99;
    Brac = 2.03;
elseif (a == 1) && (c == 5) || (a == 5) && (c == 1) || (a == 4) && (c == 5) || (a == 5) && (c == 4)
    Drac = 74.0;
    lrac = 1.60;
    Brac = 1.75;
elseif (a == 2) && (c == 5) || (a == 5) && (c == 2)
    errordlg('Parameters for this atom combination do not exist','Error','modal')
elseif (a == 3) && (c == 5) || (a == 5) && (c == 3)
    errordlg('Parameters for this atom combination do not exist','Error','modal')
elseif (a == 5) && (c == 5)
    Drac = 36.0;
    lrac = 2.67;
    Brac = 1.85;
elseif (a == 1) && (c == 6) || (a == 6) && (c == 1) || (a == 4) && (c == 6) || (a == 6) && (c == 4)
    Drac = 102.0;
    lrac = 0.96;
    Brac = 2.26;
elseif (a == 2) && (c == 6) || (a == 6) && (c == 2)
    errordlg('Parameters for this atom combination do not exist','Error','modal')
elseif (a == 3) && (c == 6) || (a == 6) && (c == 3)
    errordlg('Parameters for this atom combination do not exist','Error','modal')
elseif (a == 5) && (c == 6) || (a == 6) && (c == 5)
    errordlg('Parameters for this atom combination do not exist','Error','modal')
elseif (a == 6) && (c == 6)
    Drac = 54.0;
    lrac = 1.32;
    Brac = 2.32;
end 

% three body terms 

if (a==1) && (b==1) && (c==1) || (a==4) && (b==4) && (c==4) || (a==4) && (b==1) && (c==1) || (a==1) && (b==4) && (c==1) || (a==1) && (b==1) && (c==4) || (a==1) && (b==4) && (c==4) || (a==4) && (b==1) && (c==4) || (a==4) && (b==4) && (c==1)
    mina = 0.4;
    maxa = 2.5;
    minb = 0.4;
    maxb = 2.5;
    cont = 18;
elseif (a==1) && (b==1) && (c==2) || (a==4) && (b==4) && (c==2)
    mina = 0.5;
    maxa = 2.5;
    minb = 0.5;
    maxb = 2.5;
    cont = 18;
elseif (a==2) && (b==1) && (c==1) || (a==2) && (b==4) && (c==4)
    minb = 0.5;
    maxb = 2.5;
    mina = 0.5;
    maxa = 2.5;
    cont = 18;
elseif (a==1) && (b==1) && (c==3) || (a==4) && (b==4) && (c==3) 
    mina = 0.35;
    maxa = 3.5;
    minb = 0.9;
    maxb = 4.0;
    cont = 15;
elseif (a==3) && (b==1) && (c==1) || (a==3) && (b==4) && (c==4)
    minb = 0.35;
    maxb = 3.5;
    mina = 0.9;
    maxa = 4.0;
    cont = 15;
elseif (a==1) && (b==2) && (c==1) || (a==4) && (b==2) && (c==4)
    mina = 0.55;
    maxa = 3.5;
    minb = 0.55;
    maxb = 3.5;
    cont = 20;
elseif (a==1) && (b==3) && (c==1) || (a==4) && (b==3) && (c==4)
    mina = 0.85;
    maxa = 3.5;
    minb = 0.85;
    maxb = 3.5;
    cont = 23;
elseif (a==2) && (b==2) && (c==2)
    mina = 0.95;
    maxa = 4.0;
    minb = 0.95;
    maxb = 4.0;
    cont = 20;
elseif (a==2) && (b==2) && (c==1) || (a==2) && (b==2) && (c==4)
    mina = 0.8;
    maxa = 3.5;
    minb = 0.5;
    maxb = 4.0;
    cont = 28;
elseif (a==1) && (b==2) && (c==2) || (a==4) && (b==2) && (c==2)
    minb = 0.8;
    maxb = 3.5;
    mina = 0.5;
    maxa = 4.0;
    cont = 28;
elseif (a==2) && (b==2) && (c==3)
    mina = 0.6;
    maxa = 3.5;
    minb = 1.2;
    maxb = 5.5;
    cont = 30;
elseif (a==3) && (b==2) && (c==2)
    minb = 0.6;
    maxb = 3.5;
    mina = 1.2;
    maxa = 5.5;
    cont = 30;
elseif (a==2) && (b==1) && (c==2) || (a==2) && (b==4) && (c==2)
    mina = 0.55;
    maxa = 3.5;
    minb = 0.55;
    maxb = 3.5;
    cont = 17;
elseif (a==2) && (b==3) && (c==2)
    mina = 1.3;
    maxa = 4.5;
    minb = 1.3;
    maxb = 4.5;
    cont = 13;
elseif (a==3) && (b==3) && (c==3)
    mina = 1.6;
    maxa = 4.6;
    minb = 1.6;
    maxb = 4.6;
    cont = 17;
elseif (a==3) && (b==3) && (c==1) || (a==3) && (b==3) && (c==4)
    mina = 1.6;
    maxa = 4.5;
    minb = 0.8;
    maxb = 3.5;
    cont = 23;
elseif (a==1) && (b==3) && (c==3) || (a==4) && (b==3) && (c==3)
    minb = 1.6;
    maxb = 4.5;
    mina = 0.8;
    maxa = 3.5;
    cont = 23;
elseif (a==3) && (b==3) && (c==2)  
    mina = 1.45;
    maxa = 3.5;
    minb = 1.2;
    maxb = 4.5;
    cont = 23;
elseif (a==2) && (b==3) && (c==3)
    minb = 1.45;
    maxb = 3.5;
    mina = 1.2;
    maxa = 4.5;
    cont = 23;
elseif (a==3) && (b==1) && (c==3) || (a==3) && (b==4) && (c==3)
    maxa = 4.5;
    minb = 0.9;
    maxb = 4.5;
    cont = 18;
elseif (a==3) && (b==2) && (c==3)
    mina = 1.3;
    maxa = 3.5;
    minb = 1.3;
    maxb = 3.5;
    cont = 20;
elseif (a==1) && (b==2) && (c==3) || (a==4) && (b==2) && (c==3)
    mina = 0.6;
    maxa = 3.0;
    minb = 1.3;
    maxb = 4.5;
    cont = 23;
elseif (a==3) && (b==2) && (c==1) || (a==3) && (b==2) && (c==4)
    minb = 0.6;
    maxb = 3.0;
    mina = 1.3;
    maxa = 4.5;
    cont = 23;
elseif (a==2) && (b==1) && (c==3) || (a==2) && (b==4) && (c==3)
    mina = 0.5;
    maxa = 3.5;
    minb = 0.75;
    maxb = 3.5;
    cont = 23;
elseif (a==3) && (b==1) && (c==2) || (a==3) && (b==4) && (c==2)
    minb = 0.5;
    maxb = 3.5;
    mina = 0.75;
    maxa = 3.5;
    cont = 23;
elseif (a==1) && (b==3) && (c==2) || (a==4) && (b==3) && (c==2)
    mina = 0.8;
    maxa = 3.5;
    minb = 1.3;
    maxb = 4.5;
    cont = 23;
elseif (a==2) && (b==3) && (c==1) || (a==2) && (b==3) && (c==4)
    minb = 0.8;
    maxb = 3.5;
    mina = 1.3;
    maxa = 4.5;
    cont = 23;
elseif (a==1) && (b==1) && (c==5) || (a==4) && (b==4) && (c==5)
    minb = 1.2;
    maxb = 4.5;
    mina = 0.4;
    maxa = 4.5;
    cont = 18;
elseif (a==1) && (b==5) && (c==1) || (a==4) && (b==5) && (c==4)
    minb = 1.2;
    maxb = 4.5;
    mina = 1.2;
    maxa = 4.5;
    cont = 18;
elseif (a==5) && (b==1) && (c==1) || (a==5) && (b==4) && (c==4)
    mina = 1.2;
    maxa = 4.5;
    minb = 0.4;
    maxb = 4.5;
    cont = 18;
elseif (a==1) && (b==5) && (c==5) || (a==4) && (b==5) && (c==5)
    minb = 2.3;
    maxb = 5.0;
    mina = 1.2;
    maxa = 4.0;
    cont = 23;
elseif (a==5) && (b==1) && (c==5) || (a==5) && (b==4) && (c==5)
    minb = 1.2;
    maxb = 4.0;
    mina = 1.2;
    maxa = 4.0;
    cont = 21;
elseif (a==5) && (b==5) && (c==1) || (a==5) && (b==5) && (c==4)
    mina = 2.3;
    maxa = 5.0;
    minb = 1.2;
    maxb = 4.0;
    cont = 23;
elseif (a==5) && (b==5) && (c==5)
    minb = 2.3;
    maxb = 5.5;
    mina = 2.3;
    maxa = 5.5;
    cont = 23;
elseif (a==1) && (b==1) && (c==6) || (a==4) && (b==4) && (c==6)
    minb = 0.65;
    maxb = 3.5;
    mina = 0.4;
    maxa = 3.5;
    cont = 18;
elseif (a==1) && (b==6) && (c==1) || (a==4) && (b==6) && (c==4)
    minb = 0.65;
    maxb = 3.5;
    mina = 0.65;
    maxa = 3.5;
    cont = 18;
elseif (a==6) && (b==1) && (c==1) || (a==6) && (b==4) && (c==4)
    mina = 0.65;
    maxa = 3.5;
    minb = 0.4;
    maxb = 3.5;
    cont = 18;
elseif (a==1) && (b==6) && (c==6) || (a==4) && (b==6) && (c==6)
    minb = 1.0;
    maxb = 3.5;
    mina = 0.65;
    maxa = 3.5;
    cont = 23;
elseif (a==6) && (b==1) && (c==6) || (a==6) && (b==4) && (c==6)
    minb = 0.65;
    maxb = 3.5;
    mina = 0.65;
    maxa = 3.5;
    cont = 23;
elseif (a==6) && (b==6) && (c==1) || (a==6) && (b==6) && (c==4)
    mina = 1.0;
    maxa = 3.5;
    minb = 0.65;
    maxb = 3.5;
    cont = 23;
elseif (a==6) && (b==6) && (c==6)
    minb = 1.0;
    maxb = 3.5;
    mina = 1.0;
    maxa = 3.5;
    cont = 23;
end

