function ss = BCSN_steady(num_x,num_y,N)
%

ss = zeros(num_x,num_y);
del_x = 1/(num_x-1);
del_y = 1/(num_y-1);
x = [0:del_x:1];
y = [0:del_y:1];
for ind_x = 1:num_x
    x_i = x(ind_x);
    for ind_y = 1:num_y
        y_i = y(ind_y);
        ss(ind_x,ind_y) = BCSN_temp_analytic(x_i,y_i,N);
    end
end
