numbers = linspace(1, 100, 100);
primearray = [];

for i = 1:length(numbers)
    if isprime(i) == 1
        primearray = [primearray i];
    end
end

for j = 1:1:10000
    msize = numel(numbers);
    start = numbers(randperm(msize, 1));
    chain = []
    if start > 50 && isprime(start) == 1
        %get a new starting value
    else
        %start building chain
        chain = [chain start];
        next = numbers(randperm(msize, 1));
        if ismember(next, chain) ~= 0
            
            
        
    end
end