actually my code is taking too long to run.i am numerically
integrating equation using quad.i think reason behind code taking
so much time is that there are lot of singularities coming while
integrating the expression.So while going thro' these
discontinuities ODE45(i also used) is taking very small steps and
hence it is taking hours to run the code.
possible solution to thses kind of stiff equations is to use ODE15s.I
tried that too but its accuracy is far too low and hence it is not
showing appropriate results.
If u have any other solution to this problem ..plz suggest me.