@Rouben Rostamian  

Dear sir the parameter values are given in the image above the graph. And the related pdf is given below once check it sir. What are the points you are asking may included in pdf I think. pde_paper_.pdf

Please give the maple code how to implement it .

It is mentioned as homework.

 I am ready to type the code . How to implement for pde I m not getting. Due to that I have posted image file. 

Please update the maple code to solve .

 It will be helpful for me to generate the solution. 

After long back I have posted the question. I have familiar with ode equation. But pde little bit different. 

Please help me to solve these Questions. 

Thank you 

@acer all cf columns together and nu columns together sir

[|Gr=0.1, eta= 0.001,0.002..220, cf, nu |,|Gr=0.5, eta= 0.001,0.002..220, cf, nu |,|Gr=1, eta= 0.001,0.002..220, cf, nu |] in this manner.

Fistly for one value of Gr like Gr=0.1 , all corresponding values of cf and nu for eta= (0..20,0.01)

In next columns for the value of Gr=0.5 all corresponding values. In the similar way remaing values.

Download Demo_paper_work.mw

Dear sir,

 Here i m attaching a similar post  Here i want similar type of answers 

in the given post, the system of Odes, but in my case, i have Pdes. How do I extract similar data for my problem? 

@SHIVAS 

Try to produce one - one values of Kr by using the above answer and combine indivisual plots. I have also did like that only for this type of plots. 

May be for animation plots looping will be work or not i guess.

If it would have worked the experts would have replied for your post.

@mmcdara 

Thank you so much it is working. 

But for all values of M

how to find numerical values of diff(h(Theta,Theta) and diff(f(Theta,Theta) at Theta= Pi/4

@mmcdara 

Thank you for your answer

Actually, my problem is three-point boundary

f(-(1/4)*Pi) = 0, g(-(1/4)*Pi) = 0, h(-(1/4)*Pi) = 1, H(-(1/4)*Pi) = 1, f(0) = 1, h((1/4)*Pi) = 1, f((1/4)*Pi) = 0, g((1/4)*Pi) = 0, H((1/4)*Pi) = 1

please look at the pdf  which type of plots i need. new_paper_two_phase.pdf

i want to draw similar plots Ode_New_TWO_phase.mw

starting question i have seperated the boundary i have plotted the graphs and merged both but i m not getting similar graphs there is variations in the answers and plots.

i m giving the reference solved in matlab the boundary conditions are specified in PDF

f(±1)=0,g(±1)=0,f(0)=1,h(±1)=1,H(±1)=1

new_paper_two_phase.pdf

I have changed the boundary then graphs are comming but in base paper graphs having two phases -1 to 1 range

Ode_New_TWO_phase.mw

@C_R

This is the example i want to solve

@acer

How to draw Qz1a[j] versus gamma. At Y=0.5

I can able to draw 

Qz1a[j] versus Y for any value of gamma but reverse i am not able to draw.

Please hlep me to draw with other parameters also i will check once.

Thank you.

with(plots):
cols := [red, blue, black,green]:
 plotAA:= display
  ( [ seq
      ( odeplot
        ( Ans1[k],[gamma,Qz1a[k]],
          gamma=0..1,Y=0.5,
          color=cols[k]
        ),
        k=1..numelems(Q1Vals)
      )
    ],
    'axes'= 'boxed','linestyle' = 'solid',labels=[gamma,'Q'],size=[900,700]
  ):
if i kept like this getting error not able to evaluate

@acer

Thank you for your responce sir,

I have another  question

for any value of Y i.e, Y=0.5,

How to write plot commond for Qz1a[j] versus gamma( range of gama from 0 to 1) for different values of Q1

And Thetabz1a[j] versus gamma.

ode_Plots_1_ac.mw

 @SHIVAS 

substituted all the equations, then also no change in the final expressions. 

If any related questions please tag i will check how to impliment they implimented

I have some more equation in the worksheet how to evaluate Q and Thetab terms.

can we substitute P1=(D1*ga+C1)/(1+ga) and P2= (D1-C1)/(1+ga), Z1= -(N*(C1-D1)^2 )/(1+ga)^2,

because the solution is getting long. to reduce the length of equation by substitution with appendixed term

MaplePrimes_ODE_problem.mw

@C_R

I have checked with those posts sir.

There is no answers for this type of questions.

Thank you