求解隐式方程

guohongwei2002

a=5;
b=10;
c=0.2;
x=0:0.1:5;
y^2*((1-x^2+b*(1-2/pi*asin(1/y)+1/pi*sin(2*asin(1/y))))^2+(c*x)^2)=a^2
请问如何求解对应的y值，多谢。

liuyalong008

1. a=5;
   
2. b=10;
   
3. c=0.2;
   
4. x=0:0.1:5;
   
5. s=arrayfun(@(x)fsolve( @(y)y.^2.*((1-x.^2+b.*(1-2./pi.*asin(1./y)+...
   
6.     1./pi.*sin(2.*asin(1./y)))).^2+(c.*x).^2)-a.^2,2),x)

复制代码

布衣

solve或者fsolve函数；或者编个程序，做一维搜索