Empty Sym when using vpasolve-operator.

e = [0.05:0.05:1];
syms tp
eqn = 900-e.*5.67*10^(-8)*tp^4-20*(tp-293)==0;
S = vpasolve(eqn,tp)
I get the following message: Empty Sym: 0-by-1 when running this. Any changes I can make?

回答(1 个)

e = [0.05:0.05:1];
syms tp
for i=1:numel(e)
eqn = 900-e(i)*5.67*10^(-8)*tp^4-20*(tp-293)==0;
S{i} = vpasolve(eqn,tp);
end

2 个评论

Thank you! I want to plot S vs e (e at x-axis), and use the second element of every cell in S to do so. Any help?
e = [0.05:0.05:1];
syms tp
for i=1:numel(e)
eqn = 900-e(i)*5.67*10^(-8)*tp^4-20*(tp-293)==0;
S(i) = double(vpasolve(eqn,tp));
end
plot(e,S)
end
Maybe you should use "roots" instead of "vpasolve" because your polynomial equation has 4 zeros.

请先登录,再进行评论。

类别

帮助中心File Exchange 中查找有关 Symbolic Math Toolbox 的更多信息

标签

编辑:

2022-1-20

Community Treasure Hunt

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

Start Hunting!

Translated by