Your function will receive a lower and upper bound. It should return all pentagonal numbers within that inclusive range in ascending order. Additionally, it should return an array that indicates those numbers that are divisible by 5. For example,
[p,d] = pentagonal_numbers(10,40)
should return
p = [12,22,35] d = [ 0, 0, 1]
good
cool
so good solutions