syms x;
g = @(x) (3 + x - 2*x.^2).^(1/4);

A = -1;
B = 3/2;
figure(1);
X = linspace(A,B,100);
% x = -1 : 0.1 : 3/2;
y = g(X);
plot(X,y,'--*');
hold on
xlim([-1.5 2]);
ylim([-1.5 2]);
line([-1,-1,3/2,3/2,-1], [-1,3/2,3/2,-1,-1]);

hold off

gd = matlabFunction(diff(g(X),X));

figure(2);
line ([-5,5],[1,1]);