Clear Filters
Clear Filters

integration with known boundary an 2 variables

1 view (last 30 days)
Hi guys I have this equations that I have been trying to solve for different values of "n"
dz= - dh / [ (1 + v/k ) + (a * v/k * abs ( h ) ^ n ) ] where k and a are constants z range from 0 to a known value let's say Z and h range from a value h1 to h2 that are also known values.
Is there a way to write a matlab code that integrate this equations in term of h and then solve for the value of v
Thanks for the help in advance.
  2 Comments
Walter Roberson
Walter Roberson on 12 Feb 2020
Are dz and dh representing derivatives? If so with what independent variable? It is especially important if dz is z(h) derivative with respect to h
What is known about n? If it is symbolic you are not going to be able to solve this in the general case. You might be able to solve for n being positive integer.
Abdullah Azzam
Abdullah Azzam on 12 Feb 2020
ya dz and dh represent derivatives. Where Z represent depth and since it is goint to be integrated from 0 to desired depth I don't think it will be a problem you can consider dz as z(h). n is a constant 1,2,3,4...etc as users input.

Sign in to comment.

Accepted Answer

Walter Roberson
Walter Roberson on 12 Feb 2020
"Is there a way to write a matlab code that integrate this equations in term of h and then solve for the value of v"
No. Even for n=1 there is no closed form for the integral that can be found at all easily, so you cannot solve for v except numerically given all of the constants
There is no closed form integral for this applicable to all positive integers n. Not even if you restrict h so that it does not cross 0 to avoid the question what the integral of abs() is.
  2 Comments
Abdullah Azzam
Abdullah Azzam on 12 Feb 2020
Thanks. If possible what would be a good Numrical solution for that problem?
Walter Roberson
Walter Roberson on 12 Feb 2020
balance = @(z) integral(@(h) - 1 ./ () (1 + v./k ) + (a .* v./k * abs(h).^n), h1, h2) - Z;
v0 = rand; %or other good starting value
best_v = fsolve(balance, v0);

Sign in to comment.

More Answers (0)

Products


Release

R2017a

Community Treasure Hunt

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

Start Hunting!