Problems: Mtalab the overlapping part a function of what function may be combined image (FIG. 1), The accompanying renderings (FIG. 2) and program code desirable. Great God please help.
Matlab limited by the skill, quite cumbersome procedures, can run on R2017b version, but also look great God is not like do not spray.
Figure 1 below 
the following FIG. 2 
of Matlab program code:
%以下 x;[%Al],y:[%Mg]，MgO-二阶+[%O]=[0.0003 0.0005 0.001 0.0015]
z=0.001
f = @(x,y) -280.3.*y-1.29.*x-426.05.*z+log10(y)+log10(z)+350137.*z.^2-20050.*y.^2-0.01.*x.^2-150.*x.*y+401462.*y.*z+72.05.*x.*z++6.79
fp=fimplicit(f,[0 2 0 0.02])
fp.Color = 'r';
fp.LineStyle = '-';
fp.LineWidth = 2;
hold on
z=0.0005
f = @(x,y) -280.3.*y-1.29.*x-426.05.*z+log10(y)+log10(z)+350137.*z.^2-20050.*y.^2-0.01.*x.^2-150.*x.*y+401462.*y.*z+72.05.*x.*z++6.79
fp=fimplicit(f,[0 2 0 0.02])
fp.Color = 'g';
fp.LineStyle = '-';
fp.LineWidth = 2;
hold on
z=0.0015
f = @(x,y) -280.3.*y-1.29.*x-426.05.*z+log10(y)+log10(z)+350137.*z.^2-20050.*y.^2-0.01.*x.^2-150.*x.*y+401462.*y.*z+72.05.*x.*z++6.79
fp=fimplicit(f,[0 2 0 0.02])
fp.Color = 'black';
fp.LineStyle = '-';
fp.LineWidth = 2;
hold on
z=0.0003
f = @(x,y) -280.3.*y-1.29.*x-426.05.*z+log10(y)+log10(z)+350137.*z.^2-20050.*y.^2-0.01.*x.^2-150.*x.*y+401462.*y.*z+72.05.*x.*z++6.79
fp=fimplicit(f,[0 2 0 0.02])
fp.Color = 'blue';
fp.LineStyle = '-';
fp.LineWidth = 2;
hold on
%以下 x;[%Al],y:[%Mg]，MgAl2O4-二阶+[%O]=[0.0003 0.0005 0.001 0.0015]
z=0.001
f = @(x,y) -1124.46.*y-4.72.*x-430.52.*z-80200.*y.^2-0.04.*x.^2+350216.64.*z.^2+178826.*y.*z+978.14.*x.*z-600.*x.*y+log10(y)+2.*log10(x)+4.*log10(z)+20.54
fp=fimplicit(f,[0 2 0 0.02])
fp.Color = 'r';
fp.LineStyle = '-';
fp.LineWidth = 2;
hold on
z=0.0005
f = @(x,y) -1124.46.*y-4.72.*x-430.52.*z-80200.*y.^2-0.04.*x.^2+350216.64.*z.^2+178826.*y.*z+978.14.*x.*z-600.*x.*y+log10(y)+2.*log10(x)+4.*log10(z)+20.54
fp=fimplicit(f,[0 2 0 0.02])
fp.Color = 'g';
fp.LineStyle = '-';
fp.LineWidth = 2;
hold on
z=0.0015
f = @(x,y) -1124.46.*y-4.72.*x-430.52.*z-80200.*y.^2-0.04.*x.^2+350216.64.*z.^2+178826.*y.*z+978.14.*x.*z-600.*x.*y+log10(y)+2.*log10(x)+4.*log10(z)+20.54
fp=fimplicit(f,[0 2 0 0.02])
fp.Color = 'black';
fp.LineStyle = '-';
fp.LineWidth = 2;
hold on
z=0.0003
f = @(x,y) -1124.46.*y-4.72.*x-430.52.*z-80200.*y.^2-0.04.*x.^2+350216.64.*z.^2+178826.*y.*z+978.14.*x.*z-600.*x.*y+log10(y)+2.*log10(x)+4.*log10(z)+20.54
fp=fimplicit(f,[0 2 0 0.02])
fp.Color = 'blue';
fp.LineStyle = '-';
fp.LineWidth = 2;
hold on
