在matlab环境下的向前和向后euler算法价偏微分方程问题

				function ft(N)
				%x(n)为输入信号;N为所取点的个数
				y=zeros(1,N);
				for n=1:N
				    x(n)=cos(0.48*pi*n)+cos(0.52*pi*n);
				end
				for k=1:N
				    for m=1:N
				        y(k)=y(k)+x(m)*exp(-j*2*pi*(k-1)*(m-1)/N);
				        z(k)=norm(y(k));
				    end
				end
				plot(z)