function dupr1903
y=dsolve('D3y=2*(D2y-1)*cotg(x)','y(pi/2)=1','Dy(pi/2)=0','D2y(pi/2)=1','x')
x=-5:0.01:5
y =eval(y);
plot(x,y,'r')
end