Question: Numerical Stability

Hi all,

 I am having a very hard time with numerical stability. I am solving system of ode's (7-coupled ode's) using dsolve(stiff) and then using spline function for interpolation and finally solving system of pde's (2-coupled pde's) using pdsolve for one time step and solving all again for the next step. the solution is not stable and it requires very fine/small time step. Is there any procedure/method to improve the stability? 

 What is the stability criteria of dsolve(method=stiff(numeric)) and pdsolve(numeric)?

 

Thanks,

 Manohar 

Please Wait...