function testsqroot()
    a = [1,1,-1;1,2,0;-1,0,5];
    e = sqroot(a);
    disp(e)
endfunction