Hi. I want to solve this system of equations by varying the value of n. I managed to aolve and plot for n=1, 1.1 and 1.2 but it happens to be a problem when I let n=1.3.

>restart;

>Digits := 15;

>with(plots):n:=1.3: mu(eta):=(diff(U(eta),eta)^(2)+diff(V(eta),eta)^(2))^((n-1)/(2)):

>Eqn1 := 2*U(eta)+(1-n)*eta*(diff(U(eta), eta))/(n+1)+diff(W(eta), eta) = 0;

>Eqn2 := U(eta)^2-(V(eta)+1)^2+(W(eta)+(1-n)*eta*U(eta)/(n+1))*(diff(U(eta), eta))-mu(eta)*(diff(U(eta), eta, eta))-(diff(U(eta), eta))*(diff(mu(eta), eta)) = 0;

>Eqn3 := 2*U(eta)*(V(eta)+1)+(W(eta)+(1-n)*eta*U(eta)/(n+1))*(diff(V(eta), eta))-mu(eta)*(diff(V(eta), eta, eta))-(diff(V(eta), eta))*(diff(mu(eta), eta)) = 0;

>bcs1 := U(0) = 0, V(0) = 0, W(0) = 0;

>bcs2 := U(20) = 0, V(20) = -1;

>R1 := dsolve({Eqn1, Eqn2, Eqn3, bcs1, bcs2}, {U(eta), V(eta), W(eta)}, initmesh = 30000, output = listprocedure, numeric);

**Error, (in dsolve/numeric/bvp) precision is insufficient for required absolute error, suggest increasing Digits to approximately 23 for this problem**

>for l from 0 by 2 to 20 do R1(l) end do;

>plot1 := odeplot(R1, [eta, U(eta)], 0 .. 20, numpoints = 2000, color = red);

**Error, (in plots/odeplot) input is not a valid dsolve/numeric solution**

I have tried increasing the Digits as suggested to 23, 25, 30, 31 up untill 500 yet still same error occur suggesting to increase the Digits. Is there any other way to solve this kind of error? Can someone help me? Thank you in advance.