for i=1:3
    for j=1:3
        A(i,j)=1/(i+j-1);
    end
end
D=det(A)
B=inv(A)

function y=essai10(x)
    y=x+1;
endfunction