两方程A(x,y,z)和B(x,y)近似相等,能用matlab求出y=f(x,z)或绘出关于xyz的三维图吗有两个方程描述同一个模型第一个方程为 A(x,y,z)的形式第二个方程为 B(x,y)的形式两方程x,y的取值区间为0-0.000001,z为1

两方程A(x,y,z)和B(x,y)近似相等,能用matlab求出y=f(x,z)或绘出关于xyz的三维图吗有两个方程描述同一个模型第一个方程为 A(x,y,z)的形式第二个方程为 B(x,y)的形式两方程x,y的取值区间为0-0.000001,z为1
两方程A(x,y,z)和B(x,y)近似相等,能用matlab求出y=f(x,z)或绘出关于xyz的三维图吗
有两个方程描述同一个模型
第一个方程为 A(x,y,z)的形式
第二个方程为 B(x,y)的形式
两方程x,y的取值区间为0-0.000001,z为1-50000
A,B两个方程在这区间内是近似相等的,即 A(x,y,z)≈B(x,y)
这种情况下能求出y=f(x,z)或绘出关于xyz的三维图吗
具体方程我贴到百度空间上了,高手们帮我看看,

两方程A(x,y,z)和B(x,y)近似相等,能用matlab求出y=f(x,z)或绘出关于xyz的三维图吗有两个方程描述同一个模型第一个方程为 A(x,y,z)的形式第二个方程为 B(x,y)的形式两方程x,y的取值区间为0-0.000001,z为1
F=@(x,y,z)1./2.*((1820-1000).*z.*9.8./0.05415).*((x-(x.^2- y.^2).^0.5).^2.*...
    (x-(x-(x.^2- y.^2).^0.5)./3)-4./3.*((2000-1000)./(1820-1000)).*x.^3-y.^2.*...
    ((-2.*(y./x).*(y./x.*cos(5.*pi./6)+((x.^2-y.^2).^0.5./x).*sin(5.*pi./6)).*...
    0.05415./(1820.*x.^2.*z.*9.8)-(4./3).*2000./1820+(2./3-((x.^2-y.^2).^0.5./x)+...
    (1./3).*((x.^2-y.^2).^0.5./x).^3)+1000./1820.*(2./3+((x.^2-y.^2).^0.5./x)-...
    (1./3).*(((x.^2-y.^2).^0.5./x)).^3)).*x./((1-1000/1820).*(y./x).^2)))-...
    y.*sin(acos((x-(x.^2- y.^2).^0.5)./x-1)-(5.*pi./6));
implicitmesh(F,[0 0.000001],[0 0.000001],[1 50000],50);
-----------------------------------------------------------
function h=implicitmesh(f,xlimit,ylimit,zlimit,gd)
%implicitmesh(f,span,gd):画隐函数曲面f(x,y,z)=0的网格图,
%                          各坐标范围均限定在span=[lb,ub],
%                          网格数为gd,默认为25
%implicitmesh(f,xspan,yspan,zspan,gd):画隐函数曲面f(x,y,z)=0,
%                          各坐标范围分别限定在xspan,yspan,zspan
%h=implicitmesh(...)：画隐函数曲面并输出句柄
%例一：
%implicitmesh(inline('x.*y+z.^2'),[-5 5])%注意*\^一定要设成点运算
%例二：
%f=@(x,y,z)x.^2+y.^2+0*z-1;%注意如果f中不含某个变量一定要加上诸如0*y的项.
%implicitmesh(f,[-1 1],10)
%例三：
%f=@(x,y,z)(x.^2 + (9/4)*y.^2 + z.^2 - 1).^3 - x.^2.*z.^3 - (9/80)*y.^2.*z.^3;
%g=@(x,y,z)(sqrt(x.^2+y.^2)-2).^2+z.^2-.09;
%implicitmesh(f,[-1.5 1.5],[-.8 .8],[-1.5 1.5],50);
%hold on%可以添加图形
%h=implicitmesh(g,[-2.3,2.3]);
%colormap hsv;set(h,'facecolor','none');%可以设置各种效果
%axis off;axis equal;
if nargin==2
    ylimit=xlimit;zlimit=xlimit;gd=25;
elseif nargin==3
    gd=ylimit;ylimit=xlimit;zlimit=xlimit;
elseif nargin==4
    gd=25;
elseif nargin==5
else
    error('Error in input arguments')
end
x=linspace(xlimit(1),xlimit(2),gd);
y=linspace(ylimit(1),ylimit(2),gd);
z=linspace(zlimit(1),zlimit(2),gd);
[x,y,z]=meshgrid(x,y,z);val=f(x,y,z);
[f,v]=isosurface(x,y,z,val,0);
if isempty(f)
    warning('There is no graph in the range.');
    p=[];
else
    newplot;
    p=patch('Faces',f,'Vertices',v,'CData',v(:,3),'facecolor','w','EdgeColor','flat');
    isonormals(x,y,z,val,p);view(3);grid on
end
if nargout==0
else
    h=p;
end