Advanced Search

Arif Ullah khan

25 Reputation

4 Badges

8 years, 208 days
Contact Arif Ullah khan

MaplePrimes Activity


These are questions asked by Arif Ullah khan

How to get convergent solution of the sh...

Arif Ullah khan 25 Maple 2020
ode dsolve numeric
August 19 2020
2 20

Hi,

Need your help to find the convergent solution of the attached problem.

Here is some information about parameters

(0 <= gamma <= 10,      0 <= rho & nu <= 200)

 

Arif_Ullah.mw
 

 

> 

restart:
> 

with(plots):
> 

N:=6:
> 

 

> 

eq1 :=3*diff(f(eta),eta,eta,eta)+2*f(eta)*diff(f(eta),eta,eta)-(diff(f(eta),eta))^2;

3*(diff(diff(diff(f(eta), eta), eta), eta))+2*f(eta)*(diff(diff(f(eta), eta), eta))-(diff(f(eta), eta))^2

 (1)
> 

NULL
> 

eq2 := 3*nu*diff(g(eta),eta,eta,eta)+2*g(eta)*diff(g(eta),eta,eta)-(diff(g(eta),eta))^2;

3*nu*(diff(diff(diff(g(eta), eta), eta), eta))+2*g(eta)*(diff(diff(g(eta), eta), eta))-(diff(g(eta), eta))^2

 (2)
> 

 

> 

bc:=f(0) = 0, g(0) = 0, D(f)(0) = D(g)(0), D(D(f))(0) = -rho*nu*D(D(g))(0), D(D(f))(6) = 1, D(D(g))(N) = gamma;

f(0) = 0, g(0) = 0, (D(f))(0) = (D(g))(0), ((D@@2)(f))(0) = -rho*nu*((D@@2)(g))(0), ((D@@2)(f))(6) = 1, ((D@@2)(g))(6) = gamma

 (3)
> 

 

> 

para:={nu=3, gamma=10};

{gamma = 10, nu = 3}

 (4)
> 

 

> 

A1 := dsolve(subs(para,rho=2, {bc, eq1, eq2}), numeric,method = bvp[midrich],maxmesh=12500, output=array([seq( 0.01*i, i=0..100*N)])):
> 

 

Error, (in dsolve/numeric/bvp) Newton iteration is not converging

  
> 

pf1 := odeplot(A1, [[eta,diff(f(eta),eta), linestyle = 1, color = blue]], 0 .. N):

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

  
> 

 

> 

pf2 := odeplot(A1, [[eta,diff(g(eta),eta), linestyle = 3, color = red]], 0 .. N):

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

  
> 

display({pf1,pf2}, axes = boxed,thickness=3,labels = [eta, "f' '"],labelfont = ["ROMAN", 22,Bold,Italic],axesfont = ["ROMAN", "ROMAN", 14,Bold],axis=[thickness=3]);

Error, (in plots:-display) expecting plot structures but received: {pf1, pf2}

  

NULL


 

Download Arif_Ullah.mw

 

Find the integration of nonlinear integr...

Arif Ullah khan 25 Maple 2017
integration
January 08 2018
1 0

Hi 

I need yours help to solve the given integral 
 

> 

restart
> 

A = 3*(q-1)^(-1/(1+r))*GAMMA(q-3/(2+2*r))*GAMMA(3/(2+2*r))/(GAMMA(5/(2+2*r))*GAMMA(q-5/(2+2*r)));
> 

g := int(`&epsilon;`*(1+(`&epsilon;`/(A*t))^(1+r)/(q-1))^(-q)*c*exp(-`&epsilon;`/a), `&epsilon;` = 0 .. infinity);
> 

``
> 

NULL
> 

NULL

q,c,a and r are parameters.
Γ is a gamma function.

Download Karim.mw
 

 

 

7-th order polynomial...

Arif Ullah khan 25 Maple
polynomial solve roots
December 10 2017
4 3

Hello Maple users,
Help me to find all positive real solutions of the polynomial containing a small parameter "c" as:

(solution must be in term of c).


 

How to add random values in Maple?...

Arif Ullah khan 25
statistics sum random
November 03 2017
1 1

I have been using random numbers in Maple and want to generate a sum of that random number. I have consider the following

with(stats):

x = random[uniform](10);

and now want 

sum(x_i)^a

where (x_i) is the ith values in x (random array) and "a" is any real number.

how to handle the error "initial Newton ...

Arif Ullah khan 25
differential-equation dsolve numeric
October 18 2017
2 5

Dear Maple users,

need your help to resolve the error in attach file.
thanks

code_for_eigen_value.mw

1 2 Page 1 of 2