Tayler

n=input('n=');

x=input('x=');

a=input('a=');

sum=sqrt(a)+(x-a)*1/(2*sqrt(a));

mul=1;

number=1;

for i=2:n

mul=mul*number;

sum=sum+(x-a)^i/factorial(i)*(-1)^(i+1)*(1/2)^i*a^(1/2-i)*mul;

number=number+2;

end

sum

error=abs(sum-sqrt(x))