Info

Cette question est clôturée. Rouvrir pour modifier ou répondre.

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 vue (au cours des 30 derniers jours)
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

Réponses (1)

madhan ravi
madhan ravi le 12 Déc 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 commentaires

Cette question est clôturée.

Community Treasure Hunt

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

Start Hunting!

Translated by