Numerically evaluating an equation

[Mad_MechE]

i have the following equation:

e^-p (dp/dx) = v*(h/(h-h0)^3)

i have h0 and values for h for x = -1 to -2 i thought i could use trapz in matlab... but then i have the e^-p... how do i deal with that?

This is the data i am currently working with:
The explicit expression for h = h(x) is extremely hard to integrate and i don't think would be worth the effort. Can i do this with trapz?

x = -2:0.001:-1
h is calculated using an equation that is extremely hard to integrate.
if i run the x values through i get data like:
h = [2e-4, 3e-4, 3.2e-4...]
i have a value of h0
and i know that p(-inf) = 0
I also know the value of v


thanks

 

[Aaron Barnard]

Yes, you can use trapz in MATLAB to solve your equation. You will need to first calculate the values of e^-p at each x value using the known p(-inf) = 0 and the given h values. Then, you can use the trapz function to integrate the e^-p (dp/dx) term with respect to x. Finally, you can solve the equation for p by dividing both sides by the v*(h/(h-h0)^3) term.