MATLAB: How to update matrix values using algorithm based on position

Question

Suppose I've already create a square matrix, A=zeros(n).

Now let's say, for every position, A(i,j), I want to update the value to 1 / (i + j^2), can I do this easily?

Thanks!

Answer 1

The straightforward way is to just use a loop:

A = zeros(n);
for ii = 1:n
   for jj = 1:n
      A(ii,jj) = 1/(ii+jj^2);
   end
end

Here is another way to do it:

B = bsxfun(@(x,y) 1./(x+y.^2),(1:n).',1:n)