function E = energy(FED)
switch FED
    case 'AA',
        E = -0.9;
    case 'AU',
        E = -0.9;
    case 'UA',
        E = -1.1;
    case 'UU',
        E = -0.9;
    case 'AG',
        E = -2.3;
    case 'AC',
        E = -1.8;
    case 'UG',
        E = -2.1;
    case 'UC',
        E = -1.7;
    case 'GA',
        E = -2.3;
    case 'GU',
        E = -2.1;
    case 'CA',
        E = -1.8;
    case 'CU',
        E = -1.7;
    case 'GG',
        E = -2.9;
    case 'GC',
        E = -3.4;
    case 'CG',
        E = -3.4;
    otherwise
        E = -2.9;
end

    