[MATLAB基础] 【求助】关于matlab微分方程解的存储和调用问题

[matlabsx]

假设我已经用伦格-库塔法算出一组微分方程组的解（3个微分方程）
tspan=[0,20] %起始和终了的时间
ek4(@dy,[0,0,1],1000) %ek4是伦格-库塔子程序，[0,0,1]是初值
plot(t,y(:,1),'-')
并且能够画出其中任意一个解随时间的变化关系
问题：怎么存储任意一个解随时间的变化关系，存储成ASCII的？
好像，save t, y(:,1) /-ascii 不管用

第二个问题：
如果要算 sin(t)*y(t)，怎么办？也就是说要把时间配套起来(sin(t)和y(t)的时间)
tstep=(tspan(2)-tspan(1))/1000
for n=1,1,1000
z=sin(tstep*n)
z1=z*y(:,1)
end

可是程序报错说维数不匹配