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 Kommentare

Rik
Rik am 7 Jul. 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 am 7 Jul. 2020
Thanks
Aquatris
Aquatris am 7 Jul. 2020
what is the question?
Image Analyst
Image Analyst am 7 Jul. 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.

Antworten (0)

Diese Frage ist geschlossen.

Tags

Gefragt:

Ali
am 7 Jul. 2020

Geschlossen:

am 20 Aug. 2021

Community Treasure Hunt

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

Start Hunting!

Translated by