Effacer les filtres
Effacer les filtres

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