x=[-1:.001:1];
t= x.*(1./(1+x.^2))+(.5)*(x.^2).*(-2./((1+x.^2).^2).*x)+(1/6)*(x.^3).*(8./(1+x.^2).^3.*x.^2-2./(1+x.^2).^2);
m=.97238588.*x-.19193797.*x.^3;
MaxErrorm = atan(x)-m;
MaxErrort = atan(x)-t;
plot(x,MaxErrorm,x,MaxErrort);