function lineareqoiler
    f=dsolve('(x-2)^2*D2y-3*(x-2)*Dy+4*y=x','y(3)=0','Dy(3)=2','x');
    s=-5:0.01:5;
    g=subs(f,'x',s);
    plot(s,g)
    M=max(g)
    m=min(g)

end

