Finding Corresponding values for graph/matrix.

I want to find the value of zxt when x=3 for different speeds(v), I've using x=3 and using a second loop but the numbers don't match up!

if true % code

ib=ones(101,1);

ia=[1:1:101];

i=ib*ia;

l=6;

xa=[0:1:100];

x2=reshape(xa,101,1);

x1=ones(1,101)*l/100;

x=x2*x1;

r=98000;

u=13875;

ei=1676500000;

e=0.0;

b=e/sqrt(1-e^2);

f1=15.16;

w1=2*pi*f1;

wb=e*w1/sqrt(1-e^2);

wi=2*pi^2*i.*i*sqrt(ei/u)/(2*l^2);

wi1=sqrt(wi.^2-wb^2);

for v=1:160;

t=x/(v*0.25);

w=pi*v*0.25/l;

a=w/w1;

zxta=2*r*l^3/(ei*pi^4);

zxti=sin(i.*x.*pi./l);

zxtii=i.^2.*(i.^2.*(i.^2-a^2)+4*a^2*b^2);

zxtb=zxti./zxtii;

part1=i.^2.*(i.^2-a^2).*sin(t.*i.*w);

part2=(-a*i.*(i.^2.*(i.^2-a^2)-2*b^2)./(i.^4-b^2).^0.5).*exp(-t.*wb).*sin(t.*wi1);

part3=-2*a*b*i.*(cos(t.*i.*w)-exp(-t.*wb).*cos(t.*wi1));

zxtc=part1+part2+part3;

zxtd=zxtb.*zxtc;

zxte=sum(zxtd,2);

zxt=zxta*zxte;

temp_fdaf=max(zxt);

fdaf(v)=temp_fdaf;

end

V=[0.25:0.25:40];

plot(V,fdaf)

if true
  % code
end

Thanks!

[Image Analyst]

Don't double space. Just highlight the code and click the {}Code button. Don't click the code button first, or if code is not highlighted. Please read this: http://www.mathworks.com/matlabcentral/answers/13205-tutorial-how-to-format-your-question-with-markup