for loop for t

t=0.75
a1=2*t*x/(1+t)+t+k;
ad=t/a1;
ap=k/a1;
p(1)=1.01325*10^5;
for i=1:14;
p(i+1)=p(i)+20000;
end
sd=(Ad*p)/(2*rho*g*Veo);
sp=(Ap*p)/(2*rho*g*Vco);
kdd=sd*(ad-q*keb);
kdp=sd*ap;
kpp=sd*o*(q*kcb+ap);
kpd=sd*o*ad;
B=[1];
B(1:15)=B;
C=2*(rw*ep+ed);
C(1:15)=C;
D=(1-kdd)+(rw^2+kpp)+4*rw*ep*ed;
E=2*ed*(rw^2+kpp)+rw*ep*(1-kdd);
F=(rw^2+kpp).*(1-kdd)+kpd.*kdp;
B=B';
C=C';
D=D';
E=E';
F=F';
w=[B C D E F];
[m,n] = size(w) ;
r = zeros(m,n-1) ;
for i = 1:m
r(i,:)= roots(w(i,:))
end
r=real(r)
r=unique(r)

4 commentaires

Rik
Rik le 7 Juil 2020
This time I edited your question for you. Next time, please use the tools explained on this page to make your question more readable.
Have a read here and here. It will greatly improve your chances of getting an answer. And are you also going to click the accept button on the answer you got from KSSV? And delete the 'answer' in which you posted this question previously?
Ali
Ali le 7 Juil 2020
Thanks
Aquatris
Aquatris le 7 Juil 2020
what is the question?
Image Analyst
Image Analyst le 7 Juil 2020
No question, just an announcement.
I tried to run it but after entering random numbers for lots of missing variables, I finally gave up. It's pretty clear that Ali never read the links Rik sent.

Réponses (0)

Cette question est clôturée.

Tags

Question posée :

Ali
le 7 Juil 2020

Clôturé :

le 20 Août 2021

Community Treasure Hunt

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

Start Hunting!

Translated by