## 355 Reputation

12 years, 29 days

## with(Statistics)...

@Carl Love Is it possible to use maple to solve Eq 1

x=Prob(theta[m]>=1/x)

to get the solution

x = phi*(theta[m]-1/x)/theta[m]+(1-phi)*(1/2); ?

Fig_2.pdf

## More meaningful...

To make the picture more clearer, Here I plotted both the roots, as for as the nature of the problem goes.

## My bad for not clearing myself...

Thanks @Markiyan Hirnyk for your prompt response.

If you see the attached maple sheet, I got exactly the same output as yours.

But my worry is the lower part (y-axis, from 0..0.5) of the plot, which is missing from my

output compare to the orignal plot.

Thanks

## Exactly!...

@Preben Alsholm Thanks dear for the help

## 3D plot for f(x,y)...

@Preben Alsholm How to plot f(x,y) in 3D from eq1 only with bcs using the numerical solution?

## It worked!!...

@Preben Alsholm Dear thanks once again for your kind help. Please ignore my private msg.

## y=0?...

@Preben Alsholm Once we sub the numerical values of f(y) into P2 then why we need using y=0?

Is taking y=O has any effect on the results?

## k>0...

@Preben Alsholm Im sorry to forgot to mention that all these paramters are positive.

k>0, epsilon>0.

## accuracy issue...

@Preben Alsholm Thanks for your interest. For the following parameters, I'm facing accuracy problem.

h:=evalf(1+epsilon*sin(1)):

params:={k=0.5,B=1,epsilon=0.5};

## clue...

@Preben Alsholm Thanks dear as always for your helpful and clear response.

As you mentioned about writing a code, can you proivde any hint or any other example where a code (FDM, etc) has been adopted/used in maple to solve elliptic pdes?

Cheers!

## _C2?...

Thanks to all,

If I understood correctly, then

#_C1=0 considering symmetry and all the parameters are positive

convert(diff(eval(%,_C1=0),r), polynom)=0;

-C1/G3+(1/2)*G5*_C2*r = 0

when we sub r=0 using the bcs, then we left with nothing as far as _C2 is concern.

Am I doing something wrong?

## Thanks!...

@Alejandro Jakubi Thanks for your kind help.

## hmm............

Eq1 := (diff(psi(y), y, y, y, y))/qq-M^2*(diff(psi(y), y, y))-b*y = 0;

Eq1 := diff(psi(y), y, y, y, y)-M*(diff(psi(y), y, y))-Gr*b*y = 0;

except new parameters qq, M^2 nothing is new.

## Error, (in match/InvRed1) numeric except...

Here is my_try.mw

## op([1, 1], p1);...

@sharena2 see the attached sheet.

hydro.mw

 4 5 6 7 8 9 10 Last Page 6 of 21
﻿