function [ y ] = my_sqrt3( x0,n,a )
%UNTITLED3 Summary of this function goes here
%   Detailed explanation goes here
xnew=x0;
for k=1:n
    xold=xnew;
    xnew=1/3*(2*xold+a/(xold^2));
end
y=xnew;

end