function testsqroot()
    a = [1,2,3;4,5,6;7,8,9];
    e = sqroot(a);
    disp(e)
endfunction