666 jvbasha

javid basha jv

65 Reputation

8 Badges

3 years, 252 days

MaplePrimes Activity


These are questions asked by 666 jvbasha

Dear maple users,

Greetings.


JVB.mw
 

restart; with(plots)

fcns := {f(x), g(x), t(x)}

m := .2; pa := 3.14*(1/3); aa := .1; bb := .3; ta := .2; kt := .4; h2 := 1+m*ax+bb*sin((2*3.14)*(ax-ta)); h1 := -1-m*ax-aa*sin((2*3.14)*(ax-ta)+pa); a2 := kt+aa*sin((2*3.14)*(ax-ta))+bb*sin((2*3.14)*(ax-ta)+pa)

1+.2*ax+.3*sin(6.28*ax-1.256)

 

-1-.2*ax-.1*sin(6.28*ax-.209333333)

 

.4+.1*sin(6.28*ax-1.256)+.3*sin(6.28*ax-.209333333)

(1)

eq1 := 1.5*(diff(f(x), `$`(x, 4)))-.40*(diff((diff(f(x), `$`(x, 2)))^3, `$`(x, 2)))-.20*(diff(f(x), `$`(x, 2)))+.20*(diff(t(x), `$`(x, 1)))+.20*(diff(g(x), `$`(x, 1))) = 0

1.5*(diff(diff(diff(diff(f(x), x), x), x), x))-2.40*(diff(diff(f(x), x), x))*(diff(diff(diff(f(x), x), x), x))^2-1.20*(diff(diff(f(x), x), x))^2*(diff(diff(diff(diff(f(x), x), x), x), x))-.20*(diff(diff(f(x), x), x))+.20*(diff(t(x), x))+.20*(diff(g(x), x)) = 0

(2)

eq2 := 1.75*(diff(t(x), `$`(x, 2)))+(diff(t(x), `$`(x, 1)))*(diff(g(x), `$`(x, 1)))+2*(diff(t(x), `$`(x, 1)))^2+.7*((diff(f(x), `$`(x, 2)))^2-(diff(f(x), `$`(x, 2)))^4)+.1*(diff(f(x), `$`(x, 1)))^2+.1*t(x) = 0

1.75*(diff(diff(t(x), x), x))+(diff(t(x), x))*(diff(g(x), x))+2*(diff(t(x), x))^2+.7*(diff(diff(f(x), x), x))^2-.7*(diff(diff(f(x), x), x))^4+.1*(diff(f(x), x))^2+.1*t(x) = 0

(3)

eq3 := diff(g(x), `$`(x, 2))+2*(diff(t(x), `$`(x, 2)))-.5*g(x) = 0

diff(diff(g(x), x), x)+2*(diff(diff(t(x), x), x))-.5*g(x) = 0

(4)

bc := f(h2) = (1/2)*a2, (D(f))(h2) = 0, f(h1) = -(1/2)*a2, (D(f))(h1) = 0, t(h2) = 1, t(h1) = 0, g(h2) = 1, g(h1) = 0

f(1+.2*ax+.3*sin(6.28*ax-1.256)) = .2000000000+0.5000000000e-1*sin(6.28*ax-1.256)+.1500000000*sin(6.28*ax-.209333333), (D(f))(1+.2*ax+.3*sin(6.28*ax-1.256)) = 0, f(-1-.2*ax-.1*sin(6.28*ax-.209333333)) = -.2000000000-0.5000000000e-1*sin(6.28*ax-1.256)-.1500000000*sin(6.28*ax-.209333333), (D(f))(-1-.2*ax-.1*sin(6.28*ax-.209333333)) = 0, t(1+.2*ax+.3*sin(6.28*ax-1.256)) = 1, t(-1-.2*ax-.1*sin(6.28*ax-.209333333)) = 0, g(1+.2*ax+.3*sin(6.28*ax-1.256)) = 1, g(-1-.2*ax-.1*sin(6.28*ax-.209333333)) = 0

(5)

L := [.1]; AP := NULL; NN := nops(L); for k to NN do R := dsolve(eval({bc, eq1, eq2, eq3}, ax = L[k]), fcns, type = numeric, method = bvp[midrich], maxmesh = 1000, abserr = 10^(-4), AP); AP := approxsoln = R; X1 || k := rhs(R(0)[4]) end do; ba3 := [(X1 || (1 .. NN))(0)]

[.1]

 

1

 

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

 

R(0)

(6)

``


 

Download JVB.mw

 

How to obtain a solution for various values of "ax"

waiting for your reply.

 

 

Dear maple users @acer @Carl Love @Kitonum @Preben Alsholm @dharr @tomleslie

Greeting.

I have solved some PDEs in analytically and numerically.

But the numerical and analytical results are not matching.

I hope there is some problem with an analytical solution, especially in the first order and second-order boundary conditions.

Both the codes have enclosed here, waiting for a reply.

AN.mw

NUM.mw

Note: The PDEs are performed in the Maple 18 version. 

When compiling an analytical solution, in clarify expression click the remember table assignment.

 

Dear maple users,

Greetings.

Q1.How to plot a figure for different values of M?

      like M=1,2,3,4:

Code:Mplot.mw

Q2. How to assign all the expressions are "remember table assignment"

Dear Maple Users @acer @Carl Love @Kitonum @Preben Alsholm ,

Greetings.

How to plot a function "Am" for various values of "kt"(eg: kt=-1..1) at a point x=0.

Am vs kt.

The code has provided below.

Waiting for your response.

ktplot.mw

 

Dear maple users,

Greetings.

How to plot a contour for the below-mentioned function.

f(x):=-0.09465519086 x^3+0.02711194463 x^2+0.3862193003 x-0.00030060626-0.0003613678673 x^6-0.001538973646 x^5-0.01937304057 x^4-3.822344860 10^(-8) x^8-0.000007297718101 x^7

 

1 2 3 4 5 6 7 Page 3 of 9