方程为
1./2.*(820000.*z.*9./0.05415).*[y.^2.*(x-y./3)-4./3.*((2000000-1000000)./(1820000-1000000)).*x.^3-(((2.*x-y).*y).^0.5).^2.*((-2.*((((2.*x-y).*y).^0.5)./x).*((((2.*x-y).*y).^0.5)./x.*cos(5.*pi./6)+((x.^2-(((2.*x-y).*y).^0.5).^2).^0.5./x).*sin(5.*pi./6)).*0.05415./(1820000.*x.^2.*z.*9)-(4./3).*2000000./1820000+(2./3-((x.^2-(((2.*x-y).*y).^0.5).^2).^0.5./x)+(1./3).*((x.^2-(((2.*x-y).*y).^0.5).^2).^0.5./x).^3)+1000000./1820000.*(2./3+((x.^2-(((2.*x-y).*y).^0.5).^2).^0.5./x)-(1./3).*(((x.^2-(((2.*x-y).*y).^0.5).^2).^0.5./x)).^3)).*x./[(1000000./1820000-1).*((((2.*x-y).*y).^0.5)./x).^2])]-((((2.*x-y).*y).^0.5).*sin(acos(y./x-1)-(5.*pi./6)))=0
方程里面含x,y,z三个变量,我想算出y=f(x,z) 这种形式
但是用solve说方程复杂了,算不出结果,方程如下
syms x y z
solve(1./2.*(820000.*z.*9./0.05415).*[y.^2.*(x-y./3)-4./3.*((2000000-1000000)./(1820000-1000000)).*x.^3-(((2.*x-y).*y).^0.5).^2.*((-2.*((((2.*x-y).*y).^0.5)./x).*((((2.*x-y).*y).^0.5)./x.*cos(5.*pi./6)+((x.^2-(((2.*x-y).*y).^0.5).^2).^0.5./x).*sin(5.*pi./6)).*0.05415./(1820000.*x.^2.*z.*9)-(4./3).*2000000./1820000+(2./3-((x.^2-(((2.*x-y).*y).^0.5).^2).^0.5./x)+(1./3).*((x.^2-(((2.*x-y).*y).^0.5).^2).^0.5./x).^3)+1000000./1820000.*(2./3+((x.^2-(((2.*x-y).*y).^0.5).^2).^0.5./x)-(1./3).*(((x.^2-(((2.*x-y).*y).^0.5).^2).^0.5./x)).^3)).*x./[(1000000./1820000-1).*((((2.*x-y).*y).^0.5)./x).^2])]-((((2.*x-y).*y).^0.5).*sin(acos(y./x-1)-(5.*pi./6)))==0,y,x>0,z>0)
好像可以算出z=f(x,y),但我想算出y=f(x,z) ,求高手指点怎么列程序算出y=f(x,z)
x y z
可以把限制在一定区间内,看能得到结果不
x区间 0,0.000001
y区间 0,0.000001
z区间 0,50000,