求助，曲线画不出来。

代码：

clc;

close all;

syms x y;

y=221.9*(2*coth(2*a)-coth(a));

a=(9.27*10^-21*2.22*(10+31519.1*y))/(1.38*10^-16*x);

ezplot(y,[0,3000,0,100000000]);

title 'a';

报错：

。。。。。。。（省略很多这种数字）34193393270570125./166153499473114484112975882535043072))./(1399486262651965.*x)))./91343852333181432387730302044767688728495783936 +34193393270570125./166153499473114484112975882535043072))./(1399486262651965.*x)))./10

错误: {、[ 和 ( 的嵌套深度不能超过 32。

出错 inline/feval (line 33)

INLINE_OUT_ = inlineeval(INLINE_INPUTS_, INLINE_OBJ_.inputExpr, INLINE_OBJ_.expr);

出错 ezplotfeval (line 53)

z = feval(f,x(1),y(1));

出错 ezplot>ezimplicit (line 271)

u = ezplotfeval(f, X, Y);

出错 ezplot (line 167)

hp = ezimplicit(cax, f{1}, vars, labels, args{:});

出错 sym/ezplot (line 78)

h = ezplot(fhandle(f),varargin{:});%#ok<EZPLT>

出错 Untitled3 (line 6)

ezplot(y,[0,3000,0,100000000]);

Answer 1

[百家乐网站|网址【微8785092】]

不知道对不对：

clc;
close all;
f=@(x,y)221.9.*(2.*coth(2.*(9.27.*10.^-21.*2.22.*(10+31519.1.*y))./(1.38.*10.^-16.*x))...
    -coth((9.27.*10.^-21.*2.22.*(10+31519.1.*y))./(1.38.*10.^-16.*x)))-y;
fimplicit(f)
title 'a';