Question: not executing for three values

http://www.sciencedirect.com/science/article/pii/S100757041300508X> restart;
> l := 1; p := 1; A := .5; B := .5; pr := 1; n := [.5, 1, 1.5]; M := 0; b := .5; L := 0; s := .5; K := [blue, green];
                                      1
                                [0.5, 1, 1.5]
                                [blue, green]
> 
> for j to nops(n) do R1 := 2*n[j]/(n[j]+1); R2 := 2*p/(n[j]+1); R3 := 2/(n[j]+1);

sol1 := dsolve([diff(diff(diff(f(eta), eta), eta), eta)+f(eta)*(diff(diff(f(eta), eta), eta))+R1*(1-(diff(f(eta), eta))^2)-M*(diff(f(eta), eta)) = 0, diff(diff(theta(eta), eta), eta)+pr*f(eta)*(diff(theta(eta), eta))-R2*pr*(diff(f(eta), eta))*theta(eta)+R3*(A*(diff(f(eta), eta))+B*theta(eta)) = 0, f(0) = 1, (D(f))(0) = L+b*((D@@2)(f))(0), (D(f))(7) = 1, theta(0) = 1+s*(D(theta))(0), theta(7) = 0], numeric, method = bvp);

plots[odeplot](sol1, [eta, ((D@@2)(f))(eta)], color = red);

fplt[j] := plots[odeplot](sol1, [eta, diff(diff(f(eta), eta), eta)], color = K[j], axes = boxed); fplt[j] := plots[odeplot](sol1, [eta, f(eta)], color = K[j], axes = boxed);

tplt[j] := plots[odeplot](sol1, [[eta, theta(eta)]], color = K[j], axes = normal) end do;

plots:-display([seq(fplt[j], j = 1 .. nops(n))]);

plots:-display([seq(fplt[j], j = 1 .. nops(n))]);

plots:-display([seq(tplt[j], j = 1 .. nops(n))]);
http://www.sciencedirect.com/science/article/pii/S100757041300508X

Dear sir 

I am trying to plot the following link paper graphs for practice but I getting the plots for only one set of values here in this paper they plotted many so if you dont muned can help in this case. For example in this first graph named as Fig.1. please can you do this favour... and the paper link is  http://www.sciencedirect.com/science/article/pii/S100757041300508X