x=linspace(-1,3/2,100);
g=(3+x-2*x.^2).^(1/4);
figure(1);
plot(x,g);
hold on
x1=-1;
x2=3/2;
y1=-1;
y2=3/2;
xx=[x1,x2,x2,x1,x1];
yy=[y1,y1,y2,y2,y1];
plot(xx,yy,'r');
hold off;
figure(2)
gderiv=-(4*x-1)./(4*(-2*x.^2+x+3).^(3/4));
plot(x,gderiv);
hold on;
plot([-1,1.5],[-1,-1],'g');
plot([-1,1.5],[1,1],'g');
hold off;