aoakindele

still not giving me what i want......

Answered: aoakindele 15

December 17 2020

0 0

@Carl Love ...i have two programs not running well.....1)blesssin.mw 2) work_...ok.mw

Thank you...

Answered: aoakindele 15

December 16 2020

0 0

@Carl Love Am really grateful for your comtributions but am unable to get it... i dont know how to use method= bvp[midrich] or method= bvp[middefer] in the dsolve command....please kindly attach it affer using it on this blesssin.mw......Thank u sir

Thank u for your reply ,,,i have modify ...

Answered: aoakindele 15

November 15 2020

0 0

@tomleslie ...kindly help me outwork_please.mwwork_please.mw

yes please...

Answered: aoakindele 15

October 29 2020

0 0

@tomleslie yes i want it to evaluate to zeroyes

This have been disturbing me...

Answered: aoakindele 15

October 28 2020

0 0

work_...ok.mw...please i need your assistace on this work , dont understand the meaning of the error message neither do i understand how to fix it..........kindly find attachment and assist........THANK YOU SOOOO MUCH @Carl Love @tomleslie

Thank u for this.....helpful...

Answered: aoakindele 15

September 15 2020

0 0

@Carl Love Now i get ......thank u sir

something bothering my mind......

Answered: aoakindele 15

September 14 2020

0 0

@tomleslie Thanks so much for the solution you provided but sir how can can maple show the numerical method which it has used to solve the problem or which numerical method is used in the above problem??(is it runge kutta with shooting, shooting or what??)......Thanks

thank you for this beautiful contributio...

Answered: aoakindele 15

September 03 2020

0 0

@tomleslie it really help

waoooo...thank u sir...

Answered: aoakindele 15

September 03 2020

0 0

@Carl Love pls can u assist me on this tooerror_5.mw

ok ..thanks still having error here...ca...

Answered: aoakindele 15

September 02 2020

0 0

error.mw @Carl Love

(second derivative) F(0)...

Answered: aoakindele 15

September 02 2020

0 0

@acer F!!(0).......F prime prime(0)....not sure if thats how to write d code

it is f(eta).........error that occur wh...

Answered: aoakindele 15

September 02 2020

0 0

@acer error.mw

corrected version...

Answered: aoakindele 15

September 02 2020

0 0

>

restart;

>	with(plots): kernelopts(version);

>

#

>	# Define the ODE system

>

#

>

odeSys:= { (diff(F(eta), eta, eta, eta))*(1+epsilon-alpha*((diff(F(eta), eta, eta))^2))+F(eta)*(diff(F(eta), eta, eta))+S*(diff(F(eta), eta))-(1/2)*S*eta*(diff(F(eta), eta, eta))-(diff(F(eta), eta))^2-M*(diff(F(eta), eta)), (diff(theta(eta), eta, eta))*(1+R)-delta*(diff(F(eta), eta))^2-Pr*((3/2)*S*theta(eta)+(1/2)*S*eta*(diff(theta(eta), eta))-2*(diff(F(eta), eta))*theta(eta)+F*(diff(theta(eta), eta)))};

>

#

>	# Define the first set of boundary conditions

>

#

>	bcs1:= { F(0) = 0, (D(F))(0) = 1+lambda(1+epsilon)(D@@(F))(0) , theta(0) = 1,(1+epsilon)(D@@(F))(inf)=gammatheta(inf),F(inf)=(inf*S)/2, theta(inf) = 0

>

}:

>

>	epsilonVals:=[0.018, 0.5, 3]:

>	for k from 1 by 1 to numelems(epsilonVals) do

>	pList:=[ R = 0.5, M = 0.5,lambda=0.2,gamma=1, S = 1.5, delta = 0.3, Pr = 1.5, alpha = 0.4, epsilon = epsilonVals[k],inf=1]:

>	sol1[k]:= dsolve( eval

>	( `union`( odeSys, bcs1),

>

pList

>

),

>

numeric

>

);

>

    od:
colors:=[red,green, blue]:
  display
  ( [ seq
      ( odeplot
        ( sol1[i],
          [eta, F(eta)],
          eta=0..1,
          color=colors[i]
        ),
        i=1..numelems(epsilonVals)
      )
    ],
    color = [red, green, blue],
    title = typeset( F(eta), " =influence of epsilon on velocity profile "),
    titlefont = [times, bold, 20]
  );

>

${(diff(diff(theta(eta), eta), eta))*(1+R)-delta*(diff(F(eta), eta))^2-Pr*((3/2)*S*theta(eta)+(1/2)*S*eta*(diff(theta(eta), eta))-2*(diff(F(eta), eta))*theta(eta)+F*(diff(theta(eta), eta))), (diff(diff(diff(F(eta), eta), eta), eta))*(1+epsilon-alpha*(diff(diff(F(eta), eta), eta))^2)+F(eta)*(diff(diff(F(eta), eta), eta))+S*(diff(F(eta), eta))-(1/2)*S*eta*(diff(diff(F(eta), eta), eta))-(diff(F(eta), eta))^2-M*(diff(F(eta), eta))}$

Error, (in dsolve/numeric/bvp/convertsys) F(eta) and F cannot both appear in the given ODE

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

Download error.mw

@Carl Love

Appreciate sir .....just need fixing on ...

Answered: aoakindele 15

September 02 2020

0 0

>

restart;

>	with(plots): kernelopts(version);

>

#

>	# Define the ODE system

>

#

>

odeSys:= { (diff(F(eta), eta, eta, eta))*(1+epsilon-alpha((diff(F(eta), eta, eta))^2))+F(eta)*(diff(F(eta), eta, eta))+S*(diff(F(eta), eta))-(1/2)*S*eta*(diff(F(eta), eta, eta))-(diff(F(eta), eta))^2-M*(diff(F(eta), eta)), (diff(theta(eta), eta, eta))*(1+R)-delta*(diff(F(eta), eta))^2-Pr((3/2)*S*theta(eta)+(1/2)*S*eta*(diff(theta(eta), eta))-2*(diff(F(eta), eta))*theta(eta)+F*(diff(theta(eta), eta)))}:

>

#

>	# Define the first set of boundary conditions

>

#

>	bcs1:= { F(0) = 0, (D(F))(0) = 1+lambda(1+epsilon)(D@@(F))(0) , theta(0) = 1,(1+epsilon)(D@@(F))(inf)=gammatheta(inf),F(inf)=(inf*S)/2, theta(inf) = 0

>

}:

>

>	epsilonVals:=[0.018, 0.5, 3]:

>	for k from 1 by 1 to numelems(epsilonVals) do

>	pList:=[ R = 0.5, M = 0.5,lambda=0.2,gamma=1, S = 1.5, delta = 0.3, Pr = 1.5, alpha = 0.4, epsilon = epsilonVals[k],inf=1]:

>	sol1[k]:= dsolve( eval

>	( `union`( odeSys, bcs1),

>

pList

>

),

>

numeric

>

);

>

    od:
colors:=[red,green, blue]:
  display
  ( [ seq
      ( odeplot
        ( sol1[i],
          [eta, F(eta)],
          eta=0..1,
          color=colors[i]
        ),
        i=1..numelems(epsilonVals)
      )
    ],
    color = [red, green, blue],
    title = typeset( F(eta), " =influence of epsilon on velocity profile "),
    titlefont = [times, bold, 20]
  );