%parameter N 10 100;


%matlab
disp(u);
%end

for j = 1:1:2*N+2,
  for i = 1:1:2*N+2,
    [A(j,i)] = _ReadMatrix_Zeros_64x64();
  end
end

for l = 1:1:1,
    [x(1)] = _ReadMatrix_X();
end

for k= 1 : 1 : N,
    for j=k+1: 1 :N+1,
        for i=k: 1 :2*N+1,
            if i<=k, 
	        if i>=k, 
		    [phi]=atan(A(j,i),A(i,i));
                    [A(i,i)]=sqrt(A(i,i),A(j,i));
                    [A(j,i)]=init(x(1));
%%		else
%%                    [temp_A(1)]=funcA(A(k,i), phi, A(j,i));
%%                    [A(j,i)]=funcB(A(k,i), phi, A(j,i));
%%                    [A(k,i)]=init(temp_A(1));
                end
            else
                [temp_A(1)]=funcA(A(k,i), phi, A(j,i));
                [A(j,i)]=funcB(A(k,i), phi, A(j,i));
                [A(k,i)]=init(temp_A(1));
            end
        end
    end
end

for j = 1:1:N,
  for i = 1:1:N,
    [ ] = _WriteMatrix_Rout( A(i,j) );
  end
end