function [y] = myexpinf(x,n)
%y=1+x+x^2/2 +......x^n/factorial(n)
y=1;z=1;
for k=1:n
    z=z*x/k;
    y=y+z;
end
end