Problem 44360. Pentagonal Numbers

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]