[hwasin]

m=1:50;
y=[y,fun1(m)];
这样写会报错的吧
m的长度是50，y的长度不是，并且y没有初值。