Recent content by HoosierDaddy

Hello all, I'm trying to write a bit of Matlab to solve the Blasius Equation f*f'' + f''' = 0, where at eta = 0, f' = 0, and at eta = infinity, f' = 1. What I have so far is below...I'm a bit rusty. Two specific questions: 1. I'm trying to drive y2 to the b.c. of 1, however my loop...

Thanks D H, I believe we are required to use Heun's method. Can I just draw a parallel to Heun's method with my values? i.e. xj = tj, xj+1 = tj+1 = tj + \Deltat? If this is the case, what would my yj be? Would I have a separate yj for my p and q equations?

Not sure exactly whick RK approach it is, we will be incrementing time such that ... yj+1 = y*j+1 + \frac{h}{2}(k2 - k1) Where y*j+1 = yj + hk1 and k1 = f(xj, yj) and k2 = f(xj+1, y*j+1) Does this identify which RK2 method I am attempting to use?

Hello All, I need to utilize a Runge Kutta second order approach to solve two coupled first order DE's simultaneously given some initial conditions and a conservation relationship. The DE's are as follows: \frac{dp}{dt} = aq - bp \frac{dq}{dt} = -aq + bp Where a and b are...