请教:SAS矩阵非正定-经管之家官网！

在用SAS做矩阵运算，程序如下，但是提示有错误，麻烦各位帮忙指正。
赋值
%let c11=15440898; %let c12=1214053;%let c13=28629633; %let c14=10549932; %let c15=2781941;%let c16=59966;
%let c21=1214053;%let c22=7074056;%let c23=-1.456e8; %let c24=-5720831; %let c25=-208604;%let c26=-181805;
%let c31=28629633; %let c32=-1.456e8; %let c33=3.4837e8; %let c34=2467840;%let c35=-1624589; %let c36=6741274;
%let c41=10549932; %let c42=-5720831; %let c43=2467840;%let c44=34762183; %let c45=-1816935; %let c46=2096463;
%let c51=2781941;%let c52=-208604;%let c53=-1624589; %let c54=-1816935; %let c55=3112999;%let c56=-756615;
%let c61=59966; %let c62=-181805;%let c63=6741274;%let c64=2096463;%let c65=-756615;%let c66=1684279;
程序
proc iml;
reset print;
use rescaled_u;
read all var{c_u0 c_u1 c_u2 c_u3 c_u4 c_u5} into c_u;
s=((c_u`) * c_u)/298;
ls=(root(s))`;
sigma={&c11 &c12 &c13 &c14 &c15 &c16,
&c21 &c22 &c23 &c24 &c25 &c26,
&c31 &c32 &c33 &c34 &c35 &c36,
&c41 &c42 &c43 &c44 &c45 &c46,
&c51 &c52 &c53 &c54 &c55 &c56,
&c61 &c62 &c63 &c64 &c65 &c66};
lsigma=(root(sigma))`;
sigma2=lsigma*lsigma`;
a=(lsigma*inv(ls))`;
newu=c_u*a;
s1=((newu`)*newu)/298;
col={"new_u0" "new_u1" "new_u2" "new_u3" "new_u4" "new_u5"};
create sasuser.new_levlel_2_resd from newu [colname=col];
append from newu;
quit;
proc means data=sasuser.new_level_2_resd;
run;
错误提示
lsigma=(root(sigma))`;
ERROR: (execution) Matrix should be positive definite.