Clear Filters
Clear Filters

Info

This question is closed. Reopen it to edit or answer.

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 view (last 30 days)
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

Answers (1)

madhan ravi
madhan ravi on 12 Dec 2018
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 Comments

This question is closed.

Community Treasure Hunt

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

Start Hunting!