编写的plot语句报错了。

sharkN

  g' [6 z, Q1 l1 {& y6 G+ Z
我用一个调用函数来计算一个区间上的函数值并画图。此调用函数的输出是x = [x1;  x2;  x3];我想要把三个值的图像都画出来，编写的plot语句却都错了。错误信息：
错误使用 plot。
数据必须为可转换为双精度值的数值、日期时间、持续时间或数组。
" _. }# B- _8 u8 _5 f8 ?* E
出错 Untitled3 (line 3)
plot(y,c(:,1));
) S* M' \; ^6 S& ?4 f0 R* a' a
( h0 g/ a9 V: g! }! r. H( M以下是程序，使用的64位R2018。
- s6 h! E* i: W- B! {4 xy=9.9:0.01:10.1;
c=Solve3Polynomial(1,2.*k-20,power(k-10,2)+0.00005,-0.000001);
plot(y,c(:,1));
5 |' J' u# ?1 C" U; E" v. iplot(y,c(:,2));
plot(y,c(:,3));
function x = Solve3Polynomial(a,b,c,d)
" E: t. d/ y0 \: Z% 范盛金. 一元三次方程的新求根公式与新判别法[J]. 海南师范学院学报,1989,2(2):91-98.
. x8 w1 H1 t* ssyms x1 x2 x3;
3 {  J& V  h# s# V- l6 h) v* oA = b.*b - 3.*a.*c;   if abs(A) < 1e-14;    A = 0;  end
/ M+ k' e( }3 \1 W. ^7 w+ CB = b.*c - 9.*a.*d;   if abs(B) < 1e-14;    B = 0;  end
C = c.*c - 3.*b.*d;   if abs(C) < 1e-14;    C = 0;  end  
7 D" ?' j" g4 o4 `( W, g* BDET = B.*B - 4.*A.*C; if abs(DET) < 1e-14;  DET = 0;  end  
8 E' L6 f) }7 C3 Wif (A==0)&(B==0)
! y1 V& D5 @7 U& z+ [    x1 = -c./b;      x2 = x1 ;    x3 = x1;
end
' X( f7 ^; Z1 [: \1 \( }if DET > 0
* s, h+ N/ u  G- D% b4 w; {    Y1 = A.*b + 1.5.*a.*(-B + sqrt(DET));
0 @8 c; v( |) u+ d3 S  N    Y2 = A.*b + 1.5.*a.*(-B - sqrt(DET));
    y1 = nthroot(Y1,3);  y2 = nthroot(Y2,3);
    x1 = (-b-y1-y2)./(3.*a);
    vec1 = (-b + 0.5.*(y1 + y2))./(3.*a);  
4 n/ A" ?" O  w! e4 b    vec2 = 0.5.*sqrt(3).*(y1 - y2)./(3.*a);
% s. z; `; e# Q1 S    x2 = 0;
3 V+ A( a! w% F* z    x3 = 0;
    clear Y1 Y2 y1 y2 vec1 vec2;
! I0 A. {6 C; z" D% P8 Aend
if DET == 0 & (A ~= 0) & (B ~= 0)
$ y; R( A7 e3 d5 \    K = (b.*c-9.*a.*d)./(b.*b - 3.*a.*c); K = round(K,14);
    x1 = -b./a + K;   x2 = -0.5.*K;   x3 = x2;
9 e, \; V0 \/ Dend
if DET < 0
    sqA = sqrt(A);
( `7 J6 s- v& r& a    T = (A.*b - 1.5.*a.*B)./(A.*sqA);
, m' C9 q3 T& i    theta = acos(T);
. Z/ R! ?, p2 ?% d; o2 a    csth  = cos(theta./3);
    sn3th = sqrt(3).*sin(theta./3);
    x1 = (-b - 2.*sqA.*csth)./(3.*a);
$ ?9 u/ B" k/ d8 Q) z    x2 = (-b + sqA.*(csth + sn3th))./(3.*a);
    x3 = (-b + sqA.*(csth - sn3th))./(3.*a);
# v& h9 G1 Y) U4 ~9 U: Q    clear sqA T theta csth sn3th;
0 k6 d) H; w5 P$ a! D: @( n% Gend
, Q5 K* Z6 S9 X& v, d% D, m6 Vx = [x1;  x2;  x3];
end

麻烦大神帮忙指点一下，非常感谢！

yxlk

提示说，数据必须为可转换为双精度值的数值、日期时间、持续时间或数组。 所以你的数据是什么类型？