这个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;