What is wrong with this line of script?

1 view (last 30 days)
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? :(

Accepted Answer

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

More Answers (3)

Walter Roberson
Walter Roberson on 5 Apr 2014
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 on 5 Apr 2014
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 on 5 Apr 2014
Edited: Walter Roberson on 5 Apr 2014
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 Comments
Walter Roberson
Walter Roberson on 5 Apr 2014
Your "r" is length 11. Your "r0" is length 10. You are trying to combine two vectors of different lengths.

Sign in to comment.

Community Treasure Hunt

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

Start Hunting!