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 comentarios

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

Respuestas (0)

La pregunta está cerrada.

Etiquetas

Preguntada:

Ali
el 7 de Jul. de 2020

Cerrada:

el 20 de Ag. de 2021

Community Treasure Hunt

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

Start Hunting!

Translated by