MaplePrimes Questions

Hello!

How can I make MAPLE to create the solution of the following system?

Hi, 

I'm currently solving an equations where the boundary conditions is at infinity. I'm trying to solve it by using dsolve but i can't seem to find a solutions. Here is my equations:

ode1 := diff(f(eta), eta$3)+(diff(f(eta), eta$2))*f(eta)-(diff(f(eta), eta))^2-M . (diff(f(eta), eta))-A . (diff(f(eta), eta)+(1/2)*(eta . (diff(f(eta), eta$2)))) = 0;

ode2 := diff(theta(eta), eta$2)+Pr*(f(eta) . (diff(theta(eta), eta))-(diff(f(eta), eta)) . theta(eta)-A . (theta(eta)+1/2 . eta . (diff(theta(eta), eta)))) = 0;

and my boundary conditions are:

bcs := f(0) = 0, (D(f))(0) = 1, ((D@@2)(f))(inf) = 0, theta(0) = 1, theta(inf) = 0;

The value of Pr=7, M=1 and A=[0,1,2,4]

I really need your help, please. 

Thank you :)

The solution to the logistic map .The solution now oscillates but doesn't appear to show any discernible pattern. The value of Xn seems to "jump around". This  called chaotic.

For convenience, I am looking to extract a sequence of numbers that is generated by a simple procedure.

The attached shows such an example.

In this (simple) case, the output I require is [1,4,9,16,25].

Can anyone suggest a way to do this?

Thank you all ...

MaplePrimes_Example.mw

E-_Original_RSA.pdf

 

If I want to calculate and know the exact time duration of running step 8, do the current algorithm, 

st := time();

elapsed_time := (time() - st)*sec

is good enough? precise enough? or accurate enough? 

Any other recommendations? 

https://math.stackexchange.com/questions/301068/how-do-you-find-a-corresponding-recurrence-relation-for-some-random-algorithm/301709

according to this link, how to parse or walk through the algorithm in maple to generate recurrence relation formula?

How to remove all punctuation marks in a string read from a text file 

and then store in a list of words according to the ordering of the original text?

Ahead of the upcoming deadline 14 jan 2020 for switching from Windows 7 to Windows 10,  I am investigating whether Maple 2019 will run on an HP Stream laptop: Celeron CPU N3060 @ 1.6GHz, RAM 2GB, 64 bit.running Windows 10.

My guess is that the above spec is too low to successfully run Maple 2019.  But if anyone has had success running Maple on a similar spec device please let me know.

Melvin

I am trying to do a small procedure which converts an input of degrees into the sin function into its corresponding unit-circle value. I can get this to work. 

Sinus:= proc(g::unit)   
return evalf(sin(g*Pi/180));
end proc; 

 

but is there a more elegant way to do it? Using maybe Units package ? Thanks in advance :)

Hello!

How can I make MAPLE to put out the numerical solution of the following system?
(a =0.12, c = 47.04)  Neither solve nor fsolve does the job!

Thank for your help!

>restart:
>S:=t->c*exp(-a*t)+18;

                      S := t -> c*exp(-a*t)+18

> sys:={S(2)=55,S(8)=36};

        sys := {c*exp(-8*a)+18 = 36, c*exp(-2*a)+18 = 55}

> solve(sys);

         {a = -1/2*ln(18/37*RootOf(18*_Z^3-37)^2), c = 37*RootOf(18*_Z^3-37)}

> fsolve(sys);

         fsolve({c*exp(-8*a)+18 = 36, c*exp(-2*a)+18 = 55},{a, c})

 

I have a differential equation involving several functions of the following form:

diff(h,z) = iAf + iBg,

where h, f and g are functions of the Cartesian coordinates x, y and R and the third coordinate corresponds to z = R for some fixed constant value R.  The derivative is then with respect to the coordinate z and A and B are constants, with i the usual imaginary unit.  Is there some way this equation could be solved explicitly with Maple?

I have a system of nonlinear differential equations with initial conditions:


 

restart

de0 := {(1-p)*(diff(f(x), x, x))+p*(beta*[5]*f(x)*(diff(f(x), x))^2+beta[11]*f(x)^7+beta[6]*(diff(f(x), x))^3+beta[7]*f(x)^5+f(x)+diff(f(x), x, x)+beta[3]*f(x)^3+beta[1]*f(x)+beta[2]*(diff(f(x), x))+mu[1]*(diff(f(x), x))+(3/4)*rho*f(x)^3+beta[9]*f(x)^2*(diff(f(x), x))^3+beta[12]*f(x)^4*(diff(f(x), x))^3+beta[13]*f(x)*(diff(f(x), x))^6+beta[10]*f(x)*(diff(f(x), x))^4+beta[8]*f(x)^4*(diff(f(x), x))+beta[4]*f(x)^2*(diff(f(x), x))-(1/4)*rho*(g(x)^3-3*f(x)^2*g(x))*sin(4*Omega*x)-alpha[2]*f(x)^3*cos(2*Omega*x)-(1/4)*rho*(3*f(x)*g(x)^2-f(x)^3)*cos(4*Omega*x)-alpha[1]*g(x)*sin(2*Omega*x)-(1/2)*alpha[2]*(f(x)^3+3*f(x)^2*g(x))*sin(2*Omega*x)+(3/4)*rho*f(x)*g(x)^2-alpha[1]*f(x)*cos(2*Omega*x)-F*Omega^2*cos(Omega*x+theta)), (1-p)*(diff(g(x), x, x))+p*(diff(g(x), x, x)+beta[11]*g(x)^7+beta[7]*g(x)^5+beta[6]*(diff(g(x), x))^3+beta[1]*g(x)+beta[2]*(diff(g(x), x))+beta[3]*g(x)^3+mu[2]*(diff(g(x), x))+(3/4)*rho*g(x)^3+g(x)+beta[9]*g(x)^2*(diff(g(x), x))^3+beta[13]*g(x)*(diff(g(x), x))^6+beta[8]*g(x)^4*(diff(g(x), x))+beta[10]*g(x)*(diff(g(x), x))^4+beta[4]*g(x)^2*(diff(g(x), x))+beta[12]*g(x)^4*(diff(g(x), x))^3+beta[5]*g(x)*(diff(g(x), x))^2-(1/4)*rho*(3*f(x)*g(x)^2-f(x)^3)*sin(4*Omega*x)+(3/4)*rho*g(x)*f(x)^2-F*Omega^2*sin(Omega*x+theta)-(1/2)*alpha[2]*(f(x)^3+3*f(x)*g(x)^2)*sin(2*Omega*x)-alpha[2]*g(x)^3*cos(2*Omega*x)-(1/4)*rho*(3*f(x)^2*g(x)-g(x)^3)*cos(4*Omega*x)-alpha[1]*f(x)*sin(2*Omega*x)+alpha[1]*g(x)*cos(2*Omega*x))}:

ibvc0 := {f(0), g(0), (D(f))(0), (D(g))(0)}:

n := 3:

F := unapply(add(b[k](x)*p^k, k = 0 .. n), x):

de := map(series, eval(de0, {f = F, g = G}), p = 0, n+1):

for k from 0 to n do if k = 0 then ibvc := expand(eval[recurse](ibvc0, {f = F, g = G, p = 0})) else ibvc := {b[k](0)-1, (D(b[k]))(0)-1, (D(c[k]))(0)-1, ((D@@2)(b[k]))(0)-1, ((D@@2)(c[k]))(0)-1, c[k](0) = 1} end if; sys := `union`(simplify(map(coeff, de, p, k)), ibvc); soln := dsolve(sys); b[k] := unapply(eval(b[k](x), soln), x); c[k] := unapply(eval(c[k](x), soln), x) end do; 'F(x)' = F(x); 'G(x)' = G(x)

Error, (in dsolve) invalid input: `PDEtools/NumerDenom` expects its 1st argument, ee, to be of type algebraic, but received [0]

 

F(x) = b[1](x)*p+b[2](x)*p^2+b[3](x)*p^3

 

G(x) = c[1](x)*p+c[2](x)*p^2+c[3](x)*p^3

(1)

``


How would I deal with this error message

Download 6_11_2019_nonlinear_d_e.mw
 

>  

 

 

 

 

 

Basically it spits out the subset of values for which a division by zero error will occur for the function you specify on  range you specify for each of it's arguments, but I get an ambigous error when ever exponentiation features in the function I specify, which of course dramatically reduces the application of the calculator. Division,addition,substraction and multiplication are currently the only available arithmetic operators availble for the function window that I know the error will not occur.

If some one can help it is much appriciated

 

DIVISION_BY_ZERO_CALCULATOR.mw

Dear all

I have an first-order ode

y'(x) =0, for x in [0,1] but x different to 1/2

and we define a function y(x)= 2   if 0<=x<0.5  and y(x)=0 for 0.5<=x <=1

The attached code define these function ( maybe piecewise function is not welle defined) and how can we show that y is solution or not of the differential equation

ode.mw

many thanks

 

I was playing around with some of the maple examples(specifically the knots) and all the plots seem difficult to visualize in 3D... well, more difficult than they should.

In some cases the plots look more like projections and 2D because the graphical representation does seem to project in to the depth of the screen in a natural way.

 

Mainly it happens when rotating when some of the curve's appear to stay in front due to the highlighting and size issues.

 

I'm not sure exactly what the problem is but I imagine it definitely can be improved. It may be simply that there are no other depth cues in the plot to make things more real.

 

1. I'd like anti-aliased edges. The edges of the plot are hard and not correct coloring wise(do not fade like in real life).

2. Some type of depth perspective where the size changes with depth in to the screen in a meaningful way.

3. Proper coloring of this. Maple seems to handle the lighting correctly but it is far too weak to give the impresson of depth

 

I think maple essentially does all this(except probably AA, although there is probably a setting somewhere) but it seems they are not properly tweaked for maximum realism. I played with the lighting and projection and it seems to help a little but not much.

 

Any ideas on how to improve this? Ideally something that doesn't have to be set every plot.

First 622 623 624 625 626 627 628 Last Page 624 of 2427