% energyofneighbours.m
% computes the energy of a vector and each of its immediate neighbours
%
% omega = inhibition constant
% inVector = current state

function u = energyofneighbours (omega, inVector)
u = inVector;
n = 9;
i = 0;
while i <= n
    if i == 0
        u
        energy(omega,u);
        i = i+1;
    else
        uNeighbour = u;
        uNeighbour = flipbit(uNeighbour, i)
        energy(omega, uNeighbour);
        i = i+1;
    end
end