function [pascal_square] = pascal_square_function(n)
pascal_square(1,1:n)=1;
pascal_square(1:n,1)=1
for i=2:n,
    for j=2:n;
        pascal_square(i,j)=pascal_square(i,j-1)+pascal_square(i-1,j);
    end
end             
end