新手求助：关于微分方程的问题

[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

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

[yisanming]

问题二
求sin(t)*y(t),首先表达出t，在M文件或命令窗口输入sin(t).*y(t)，得到的是t值对应的sin值与y(t)乘机，如果想要得到矩阵，输入sin(t)*y(t)'就行了。