【求助】 subs 结果出现NAN值

wrp19880903

计算结果中
ss的前470项为正常数值，之后的全部为NAN结果，程序改过别的方法结果还是一样，
求教应对代码做那些改动啊！！！！

程序代码
syms t y;
w=[-(cos((1751642673090177*y)/1125899906842624) + cosh((1751642673090177*y)/1125899906842624) - (7776840453599041*sin((1751642673090177*y)/1125899906842624))/288230376151711744 - (2451593436110691*sinh((1751642673090177*y)/1125899906842624))/2251799813685248)*(0.000000015303395673689924359312878260763*cos((8466132921799293*t)/70368744177664) - 0.0000013908183652339644845780398696329*cos(120.31098495127839953511283965781*t) + 0.0000035702965949742052776345927972124*sin(120.31098495127839953511283965781*t) + (0.00000090530817636565783318086692702176*57.0^(371.19821450474597668289788998663*t))/200.0^(371.19821450474597668289788998663*t) + 0.00000048551018886830665139717294261118*0.285^(12.659496983011933934903936460614*t)), (cos((4311625730013215*y)/1125899906842624) + cosh((4311625730013215*y)/1125899906842624) - (8372372071790843*sin((4311625730013215*y)/1125899906842624))/9007199254740992 - (281675502569541*sinh((4311625730013215*y)/1125899906842624))/281474976710656)*(0.000000041773140401032146699563087539053*cos((400743225491141*t)/549755813888) - 0.00000080132541498326576986609405731528*cos(728.94768071117323415819555521011*t) + 0.0000004895496691306274156693072141733*sin(728.94768071117323415819555521011*t) + (0.00000076461151627112811593430431238547*57.0^(371.19821450474597668289788998663*t))/200.0^(371.19821450474597668289788998663*t) + 0.000000036713898712137653931789744929817*0.285^(12.659496983011933934903936460614*t)), -(cos((1520975453922935*y)/281474976710656) + cosh((1520975453922935*y)/281474976710656) - (8778660859438457*sin((1520975453922935*y)/281474976710656))/9007199254740992 - (9006905136293105*sinh((1520975453922935*y)/281474976710656))/9007199254740992)*(0.0000000048910494591570518238497395192078*cos((1595800028228143*t)/1099511627776) - 0.000000029922672990943306078751688396495*cos(1451.3716707625853814533911645412*t) + 0.0000000092700901151303315125277797514866*sin(1451.3716707625853814533911645412*t) + (0.000000028837934247896044746824631714408*57.0^(371.19821450474597668289788998663*t))/200.0^(371.19821450474597668289788998663*t) + 0.0000000010847387430472613319270566820876*0.285^(12.659496983011933934903936460614*t)), (cos((7834202817773005*y)/1125899906842624) + cosh((7834202817773005*y)/1125899906842624) - (4450158003905025*sin((7834202817773005*y)/1125899906842624))/4503599627370496 - (2251803038868975*sinh((7834202817773005*y)/1125899906842624))/2251799813685248)*(0.0000000034295570018874148957047054310808*cos((5292168850630165*t)/2199023255552) - 0.000000008095990755154777591373999283846*cos(2406.5997652679311613610479980707*t) + 0.0000000015157626918510235815115914525539*sin(2406.5997652679311613610479980707*t) + (0.0000000078193323993627822895847959037324*57.0^(371.19821450474597668289788998663*t))/200.0^(371.19821450474597668289788998663*t) + 0.00000000027665835579199530178920338011357*0.285^(12.659496983011933934903936460614*t)), ((8948734245472809*sin((1197461490573837*y)/140737488355328))/9007199254740992 + (9007198694100851*sinh((1197461490573837*y)/140737488355328))/9007199254740992 - cos(1197461490573837/140737488355328*y) - cosh(1197461490573837/140737488355328*y))*(0.00000000076222374647172140049377353625877*cos((7913109150916167*t)/2199023255552) - 0.00000000082063972017434245036928661876307*cos(3598.465423654564347089035436511*t) + 0.00000000010282198352628437867157106667534*sin(3598.465423654564347089035436511*t) + (0.00000000079313724517636248178832848094935*57.0^(371.19821450474597668289788998663*t))/200.0^(371.19821450474597668289788998663*t) + 0.000000000027502474997979968580958137813721*0.285^(12.659496983011933934903936460614*t))];
S=0;
for i=1:3
    S=S+w(i);
end
S
y=1.015;
S=subs(S);
for i=1:1001
    tt(i)=1*(i-1)/1000;
    t=tt(i);
    SS(i)=subs(S);
end计算结果中
ss的前470项为正常数值，之后的全部为NAN结果，程序改过别的方法结果还是一样，
求教应对代码做那些改动啊！！！！