Please proofread for a newbie :)
显示 更早的评论
Hi, I am having so much trouble with matlab, this is what I have been trying:
function P=PR(V);
%Peng-Robinson EOS
%used to calculate vg and vl
R=8.314;
T=313.15;
b=5.63e-5;
a=1.1188;
P=1.3e6;
P=((R*T)./(V-b))-(a./(V*(V+b)+(b*(V-b))));
end
and then i do:
fsolve(@PR,1)
as I am trying to solve for V
it gives an answer that I know is not right, as I had the m-file working at uni today and now at home it just does not work. Any help or advice would be great, thanks
1 个评论
Jan
2011-4-1
Some advice: 1. Use code formatting in your question to improve readability (see: http://www.mathworks.com/matlabcentral/answers/help/markup ). Readers are more likely to answer, if they see, that you spend time to formulate the question as clear as possible. 2. If you know, that result is not correct, please post the correct value in addition! How should we find the error without knowing the correct result?! 3. "It does not work" is a magic keyword: Whenever this appears in a question, it is most likely impossible to create a meaningful answer. We cannot guess, what the problem is, if you do not describe it with all necessary details -- except for walter, who obviously owns a magic crystal ball :-)
采纳的回答
Jan
2011-4-1
1 个投票
"fsolve(@PR, 1)" solves the equation "P=((R*T)./(V-b))-(a./(V*(V+b)+(b*(V-b))))" such that P==0 using the initial value V==1. Do you want to do this?
It worked in the University and now it does not calculate the correct answer? This would be very helpful: Then the error is inside the modifications you have inserted in between! Or the command FSOLVE of one of the computers is broken - but this is less likely.
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 MATLAB 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
