这个 MATLAB 代码有什么问题？向量化 For 循环

Question

我想计算1/((1^2)(3^2))+ 1/((3^2)(5^2))+1/((5^2)(7^2))+ ...的总和

以下代码是我在不循环的情况下计算总和的尝试：

n = 1000;
t0 = clock;    %vectorizing we get 
x1 = (1./(1:2:n)).^2;
y1 = (1/(3:2:n)).^2y1.*x1;
t = sum(y1.*x1)etime(clock,t0)

但 MATLAB 给出向量大小不一致的错误。

如何解决此问题以获得正确答案？

Answer 1

y1=(1/(3:2:n)).^2y1.*x1;

应该是：

y1=(1/(3:2:n+2)).^2;