@tomleslie when i change the another value for another parameter. It still have error, i try my method, and your method. 

@Preben Alsholm thanks for that.. I just change position B=0 and M=L[k] and L :=[.1, .3, .6, .9, 1.2] and its have same error again..

@Preben Alsholm new_probs.mw here

@Preben Alsholm 

@tomleslie I just change the value for S=-.5, but i have this error 

"Error, (in dsolve/numeric/bvp) initial Newton iteration is not converging"

@Carl Love when I put B*.. this error happen 

Error, (in dsolve/numeric/bvp) Newton iteration is not converging

@Carl Love but its different with value that my lecturer get.. he get [-1.0001,-1.0519,-1.1019] for 1st 3 value

@Preben Alsholm I did like what you said, remove the line N:=1, but after i assign q=....... I get graph like y=3 for example..

@Preben Alsholm what means do not assign to N and Pr? i not understand..

@Preben Alsholm I already get the graph, thanks for helping.. Do you have idea how to plot graph theta (0) vs N, when Pr is fix to 1 and 7..

@Preben Alsholm I just try command that you give..but its error..can you give me the worksheet?

@tomleslie  value theta (0) and I want smooth graph

@tomleslie I want graph Y1 vs L

@Preben Alsholm thanks for your help.. I got the graph, but its not smooths.. like graph theta against eta

@Preben Alsholm Im sorry, I want to plot graph theta against Pr.. Pr will be x-axis and theta(0) will be y-axis..