[yisanming]

x=linspace(eps,pi,100);
y=acot(x);
plot(x,y);