Solving nonlinear equation inside for loop

6 次查看（过去 30 天）
Kim Ibrahim
Kim Ibrahim 2021-11-23
评论： Kim Ibrahim 2021-12-2
I want to solve for p in the following equation: a1*P^0.5+a2*p+a3*p^1.5=k
if I need to iterate K from 1 to 10000 this is the code:
a1=1.4*10^25;
a2=2.7*10^24;
a3=1.311*10^23;
k=linspace(1,10000,10000);
syms p
for i=1:10000
p(i)=solve((a1*p(i)^0.5)+(a2*p(i))+(a3*p(i)^1.5)==k(i))
end
however I get the following error: Index exceeds the number of array elements (1).
Thank you.
  3 个评论
显示 1更早的评论
Dyuman Joshi
Dyuman Joshi 2021-11-24
You are expecting 3 solutions for each p and k, right?
Kim Ibrahim
Kim Ibrahim 2021-11-30
yes, this equation gives 3 solutions for p when I give k a value at each itertion. the question is i can`t make the dimentions right beause for the 10000 iteration eah give out 3 solutions. That`s the error I get
"Unable to perform assignment because the left and right sides have a different number of elements."

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Matt J
Matt J 2021-11-30
编辑：Matt J 2021-12-1
Your equations are really polynomials in p^0.5. Therefore, you can do,
a1=1.4*10^25;
a2=2.7*10^24;
a3=1.311*10^23;
k=linspace(1,10000,10000);
p=cell(1,numel(k));
for i=1:10000
coeffs=[a3,a2,a1,-k(i)]; %polynomial coefficients
p{i}=sqrt(roots(coeffs));
end

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Linear Algebra 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by