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 Comments

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

Answers (0)

This question is closed.

Tags

Asked:

Ali
on 7 Jul 2020

Closed:

on 20 Aug 2021

Community Treasure Hunt

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

Start Hunting!