Info

此问题已关闭。 请重新打开它进行编辑或回答。

why my first for loop in the following code is not working ?i'm able to run the program and i got the results for the last iteration i.e '650' .how can i get the result for all the iterations at once.

1 次查看(过去 30 天)
for T=50:30:650
S=0.005;
r=90;
Q=1;
t=6.94*10^-5:0.0299:0.75;
u=(r^2)*S./(4.*T.*t);
if u<1
w=-(0.5772)-(log(u))+u-(u.^2/(2*factorial(2)))+(u.^3/(3*factorial(3)))-(u.^4/(4*factorial(4)))+(u.^5/(5*factorial(5)));
else
w=(2.718.^-u./(u)).*[(u.^2+2.334733.*u+0.250621)./(u.^2+ 3.330657.*u+1.681534)];
end
s=(Q*w./(4*3.14.*T));
d=s'
u1=u';
end

回答(1 个)

madhan ravi
madhan ravi 2018-12-12
Your code works fine: (loop is superfluos though)
T=50:30:650;
S=0.005;
r=90;
Q=1;
t=6.94*10^-5:0.0299:0.75;
u=(r^2)*S./(4.*T.*t);
if u<1
w=-(0.5772)-(log(u))+u-(u.^2/(2*factorial(2)))+(u.^3/(3*factorial(3)))-(u.^4/(4*factorial(4)))+(u.^5/(5*factorial(5)));
else
w=(2.718.^-u./(u)).*[(u.^2+2.334733.*u+0.250621)./(u.^2+ 3.330657.*u+1.681534)];
end
s=(Q*w./(4*3.14.*T));
d=s.
u1=u.';
  6 个评论

此问题已关闭。

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by