What is wrong with this line of script?

1 次查看（过去 30 天）
Alice
Alice 2014-4-5
评论： Alice 2014-4-5
I've written: [t,h] = ode45('cylinder_tank_drain', [0, 200], h0)
it keeps coming up with an error message what just says there is an error in this line. what is wrong? :(

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

采纳的回答

Youssef Khmou
Youssef Khmou 2014-4-5
Change r0 to :
r0=linspace(0.01,0.1,length(r));

更多回答（3 个）

Walter Roberson
Walter Roberson 2014-4-5
There is no obvious error in that line. Perhaps there is an error in the previous line that is affecting it. Please post more of your script.
Also, I suggest converting to function handle form:
[t, h] = ode45(@cylinder_tank_drain, [0, 200], h0);
Youssef Khmou
Youssef Khmou 2014-4-5
The problem resides in the function "cylinder_tank_drain", if it is possible you post the function to study its input and outputs.
Alice
Alice 2014-4-5
编辑：Walter Roberson 2014-4-5
function [ dhdt ] = cylinder_tank_drain( t, h)
%draining of a uniform cross sectional area tank (cylinder)
r = [1:0.5:6];
A = pi.*r.^2;
r0 = [0.01:0.01:0.1];
a0 = pi.*r0.^2;
g = 9.8;
dhdt = (-a0.*sqrt(2*g.*h))./A;
  2 个评论
Walter Roberson
Walter Roberson 2014-4-5
Your "r" is length 11. Your "r0" is length 10. You are trying to combine two vectors of different lengths.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by