I'm working on a time-dependent equation from fluid mechanics for the velocity (u) in a loop. The equation I'm trying to solve is:
Where c1, c2, c3, and c4 are all positive, known constants. I've tried pretty much all of the ODE solvers in MATLAB and they all blow up to infinity and return the following error: "Warning: Failure at t=3.237651e+00. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (1.150245e-14) at time t."
Can anyone give me some insights into other solvers that might work for this type of problem, or modifications to the solvers I've tried that can resolve the integration tolerance issue?