array = ConstantArray[0, {10, 1000}]

array[[i, 1000]] = some_expression

Do[
  h = j*k;
  density = expression1;
  phi(r) = expression2;
  Delta = expression3;
  bin = Table[(i - 0.5)*Delta, {i, 1000}];
  generalP = ConstantArray[0, {10, 1000}];
  generalP[[j, 1000]] = some_expression_with_bin[[j]]_&_h,
  {j, 10, 1, -1}]

generalP[[1,1000]]=some_expression_with_bin_&_h(1)
generalP[[2,1000]]=0
generalP[[3,1000]]=0
...
generalP[[10,1000]]=0

Do[generalP[[j, 1000]] = j - 1, {j, l, 1, -1}]