@Carl Love I have no idea where I can change it. 

Incrs:= [.01, .1, 0.1, 0.1]: #Increments in a, b, c and d, respectively

As you mentioned the above is increments in a, b, c, d

nIncrs:= [2,2,2,2]:  #Numbers of increments

This is the upper bound.

Now where I can change the lower bound?

@Carl Love You mean that we cannot choose the lower bound freely?

@Preben Alsholm, It seems that when you increase no of iterations (instead of 8, if we use 12) then there is no two solutions.

@Carl Love I want the parameters to vary in the interval [0, 2], where the lower bound is 0 and the upper bound is 2 with increment 0.01

One thing more, for a huge data, the loop take long time to execute, anyway to make it fast.  

@Carl Love Thx 

Its working fine in M17.

I have another question, say if we want to vary a,b,c,d from 0 to 2 then?

@Carl Love 

If there is no way around then I have to install M17, right?

@Preben Alsholm Thx

The next step will be to loop sol for a, b, c, d?

@Carl Love 

Thx for your response but I am getting this Error for the loop or I missunderstood something?

 for row while not It[finished] do 

 V:= It[nextvalue]();
 A[row,..]:= < V[], eval(diff(f(eta),eta$2), Sol(V[])(0)) >
 end do:

Error, (in dsolve/numeric/bvp) invalid orientation specification

I suspect that there is (no) exact solution possible for your ode?

Thx @Preben Alsholm for the correction.

Why don't you try the numeric?

restart:with(DEtools):

ODE := diff(y(t),t) = t*(1-0.3*t)-t*y(t)/(1+0.6*t);

ics := [y(0) = 1];

DEplot(ODE,y(t),t=0..5,ics,arrows=medium,linecolor=black);

@Carl Love Whats the possible interpretation of such results?

Thanks @rlopez. Have you had the chance to look to my "possible explanation"?

@rlopez How to extract the intersecting point upto 5 decimal places from

plot([R,C],.1..2);?

@taro yamada 

v := int(q(i)^2,i = 0 .. M);

    =  Q(M)-Q(0) which is a constant

How you can diff v w.r.t q(i)?

@taro yamada