2 years, 83 days

## can any one solve this?...

Maple 2019

 >
 > restart:
 >
 (1)
 > eq11 := (10-9*lambda1)*n*diff(f(eta),eta,eta,eta)*(-diff(f(eta),eta,eta))^(n-1)-m*diff(f(eta),eta)^2+M^2*diff(f(eta),eta)+(m*(2*n-1)+1)/(n+1)*f(eta)*diff(f(eta),eta,eta)=0;
 (2)
 > eq21 := (10-9*lambda1)*n*diff(theta(eta),eta,eta)*(-diff(theta(eta),eta))^(n-1)+(m*(2*n-1)+1)/(n+1)*f(eta)*diff(theta(eta),eta)-m*diff(f(eta), eta)*theta(eta)=0;
 (3)
 > bc21:= f(0)=0, D(f)(0)=1+alpha*(D@@2)(f)(0), D(f)(N)=0, theta(0)=1+b*D(theta)(0), theta(N)=0;
 (4)
 >
 > n:=1.2: b:=0: M:=0: alpha:=0:
 > N:=6:
 >
 > m:=0.8:
 > A11:=dsolve({eq11,eq21,bc21},numeric,continuation=lambda1,maxmesh=15000, output=array([seq( i, i=0..N,0.01 )])):
 > m:=1.2:
 > A12:=dsolve({eq11,eq21,bc21},numeric,continuation=lambda1,maxmesh=15000, output=array([seq( i, i=0..N,0.01 )])):
 > m:=2:
 > A13:=dsolve({eq11,eq21,bc21},numeric,continuation=lambda1,maxmesh=15000, output=array([seq( i, i=0..N,0.01 )])):
 >
 >
 > with(plots):
 > p1:=odeplot(A11, [[eta,f(eta),color=red,linestyle=2]],0..N): p2:=odeplot(A12, [[eta,f(eta),color=blue,linestyle=2]],0..N): p3:=odeplot(A13, [[eta,f(eta),color=green,linestyle=2]],0..N):
 >
 > display({p1,p2,p3}, axes=boxed,  title=`velerature Profiles:`);
 >
 >
 > for j from 1 to 100*(N-1) do eta[j]:=A11[2,1][j,1] end do:
 > for j from 1 to 100*(N-1) do vel1[j]:=A11[2,1][j,2] end do:
 > for j from 1 to 100*(N-1) do vel2[j]:=A12[2,1][j,2] end do:
 > for j from 1 to 100*(N-1) do vel3[j]:=A13[2,1][j,2] end do:
 >
 > with(LinearAlgebra): E1:= <|||>;
 (5)
 >
 > #ExportMatrix("F:/D/Drive D/Local/Rashad/Paper 9/vel1.txt", E1, format = rectangular):
 >

 > with(plots):
 > p1:=odeplot(A11, [[eta,theta(eta), color=red,linestyle=2]],0..N): p2:=odeplot(A12, [[eta,theta(eta), color=blue,linestyle=2]],0..N): p3:=odeplot(A13, [[eta,theta(eta), color=green,linestyle=2]],0..N):
 >
 > display({p1,p2,p3}, axes=boxed,  title=`temperature Profiles:`);
 >
 >

 > for j from 1 to 100*(N-1) do eta[j]:=A11[2,1][j,1] end do:
 > for j from 1 to 100*(N-1) do temp1[j]:=A11[2,1][j,5] end do:
 > for j from 1 to 100*(N-1) do temp2[j]:=A12[2,1][j,5] end do:
 > for j from 1 to 100*(N-1) do temp3[j]:=A13[2,1][j,5] end do:
 >
 > with(LinearAlgebra): E2:= <|||>;
 (6)
 > #ExportMatrix("F:/D/Drive D/Local/Rashad/Paper 9/temp1.txt", E2, format = rectangular):
 >
 >