Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

I am having touble plotting the equation, diff(y(t), t$3)+3*(diff(y(t), t$2))+4*(diff(y(t), t))+12*y(t) = 0. With the intional conditons y'(0)=y''(0)=0 and y(0)=3.

Currently I am getting this error, "Error, (in dsolve/numeric/type_check) insufficient initial/boundary value information for procedure defined problem."

My current eqution is DEplot(ode, y(t), t = -1 .. 5, y(t) = -5 .. 5, [ivp]), with ivp := [((D@D)(y(t)))(0) = 0, (D(y(t)))(0) = 0, (y(t))(0) = 3].

HI

please help me for dsolve this nonlinear differential equations

thanls...

HAB.mw
 

restart; Digite := 100; Phi0 := 5; A := b*h; g13 := 31/250000000; g1 := 113/500000; f13 := 1/1000000000; c1 := 226000000000000; b := 10*10^(-9); J := (1/12)*b*h^3; h := 15*10^(-9); L := 100*10^(-9); E1 := (339/10000000000000000000000)*(diff(u(x), x, x, x, x))+(1017/10000000000000000000000)*(diff(w(x), x, x))*(diff(w(x), x, x, x))+(339/10000000000000000000000)*(diff(w(x), x, x, x, x))*(diff(w(x), x))-(339/10000)*(diff(u(x), x, x))-(339/10000)*(diff(w(x), x))*(diff(w(x), x, x)) = 0

E2 := -(1017/1600000000000000000000000000000000000000)*(diff(w(x), x, x, x, x, x, x))+(1589109/2500000000000000000000000)*(diff(w(x), x, x, x, x))-(339/10000*(diff(u(x), x, x)+(diff(w(x), x))*(diff(w(x), x, x))))*(diff(w(x), x))-(diff(w(x), x, x))*((339/10000)*(diff(u(x), x))+(339/20000)*(diff(w(x), x))^2+0.5824000000e-4)+(339/10000000000000000000000)*(diff(w(x), x, x))*(diff(u(x), x, x, x)+(diff(w(x), x, x))^2+(diff(w(x), x, x, x))*(diff(w(x), x)))+(339/10000000000000000000000)*(diff(w(x), x))*(diff(u(x), x, x, x, x)+3*(diff(w(x), x, x))*(diff(w(x), x, x, x))+(diff(w(x), x, x, x, x))*(diff(w(x), x)))-2 = 0:
 

E3 := -5.385803274*10^(-17)*(diff(Phi(x), x, x))+2.659881780*Phi(x)-5.125107476*10^(-20)*(diff(psi(x), x, x))+1.146681319*psi(x)+3.300000000*10^(-8)*(diff(w(x), x, x)) = 0:

E4 := -5.125107476*10^(-20)*(diff(Phi(x), x, x))+1.146681319*Phi(x)+(891/100000000000000)*(diff(psi(x), x, x))/Pi+34976.39822*psi(x)+0.4351500000e-5*(diff(w(x), x, x)) = 0:

dsys3 := {EQ1, EQ2, EQ3, EQ4, c1*J*((D@@2)(w))(0)+A*g13*((D@@2)(w))(0)-2*b*f13*Phi0-g1*J*((D@@4)(w))(0)+g1*A*(((D@@2)(u))(0)+((D@@1)(w))(0)*((D@@2)(w))(0))*((D@@1)(w))(0) = 0, c1*J*((D@@2)(w))(L)+A*g13*((D@@2)(w))(L)-2*b*f13*Phi0-g1*J*((D@@4)(w))(L)+g1*A*(((D@@2)(u))(L)+((D@@1)(w))(L)*((D@@2)(w))(L))*((D@@1)(w))(L) = 0, Phi(0) = 0, Phi(L) = 0, psi(0) = 0, psi(L) = 0, u(0) = 0, u(L) = 0, w(0) = 0, w(L) = 0, (D(u))(0) = 0, (D(u))(L) = 0, ((D@@2)(w))(0) = 0, ((D@@2)(w))(L) = 0}; dsolve(dsys3, numeric, initmesh = 3024, abserr = 0.1e-4)

Error, (in dsolve/numeric/process_input) system must be entered as a set/list of expressions/equations

 

``


 

Download HAB.mw

 

restart; with(Statistics):
X := RandomVariable(Normal(0, 1)): Y := RandomVariable(Uniform(-2, 2)):
Probability(X*Y < 0);

crashes my comp in approximately 600 s. Mma produces 1/2 on my comp in 0.078125 s.

Why Deplot3d does not work in Maple 2016?

For example, 

with(DEtools);

DEplot3d({diff(y(x), x) = y(x)-z(x), diff(z(x), x) = z(x)-2*y(x)}, {y(x), z(x)}, x = 0 .. 3, [[y(0) = 1.638, z(0) = 2.31]], y = 0 .. 2, z = -4 .. 4, scene = [x, z(x), y(x)], linecolor = COLOR(HUE, .5));

There is nothing shown in Maple 2016.

Hi all,

 how to calculate this integral:

int(sqrt(a^2+cos(x)),x=0..Pi);

 

I uploaded this file AnimationTest.mw to the Maple Cloud. One a 3d animation and one just a 3d plot. When I open the file in the Maple Cloud from my browser,  I cannot rotate either plot with the mouse. Is there a way to change the worksheets so rotation from the browser will be possible?

I Could Not Write An If Then Or Ifelse Statement. Please Help Me.

f := unapply(x^2-2, x); a := 1; b := 2; n := 10; Digits := 10;
      2    
x -> x  - 2
                               1
                               2
                               10
                               10
c := evalf(eval((a*f(b)-b*f(a))/(f(b)-f(a))));
                          1.333333333
if  f(c)*f(a)<0 then ;
          "      k:=evalf(eval(|(f(c))/(b-c)|)) and "

                          /(1 + k) a f(b) - b f(a)\
             x[i] := evalf|-----------------------|
                          \  (1 + k) f(b) - f(a)  /
            "     elif f(x[i])*f(a)<0 then b:=x[i]"
                 "     else b:=c and a:=x[i] "
                  "     if f(c)*f(a)>0 then "
                 "      k:=|(f(c))/(b-c)|and "

                          /a f(b) - b f(a) (1 + k)\
             x[i] := evalf|-----------------------|
                          \  f(b) - f(a) (1 + k)  /
            "     elif f(x[i])*f(a)>0 then a:=x[i]"
              "     else a:=c and b:=x[i] end if"

Error, unterminated 'if' statement
     Typesetting:-mambiguous(Typesetting:-mambiguous(

       if fApplyFunction(c)sdotfApplyFunction(a)lt0 then , 

       Typesetting:-merror("unterminated 'if' statement")))

how to solve  Fourier Series on maple ??

Could you tell about manual (book) at maple which tells how to make calculations in quantum field theory (Grassmann algebra, a Lie algebra, producing functional, fermionic determinants) and high-temperature quantum field theory (partition functions, thermodynamic potentials)?

Let us consider

with(Statistics):
X1 := RandomVariable(Normal(0, 1)):
X2 := RandomVariable(Normal(0, 1)):
X3 := RandomVariable(Uniform(0, 1)): 
X4 := RandomVariable(Uniform(0, 1)):
Z := max(X1, X2, X3, X4); CDF(Z, t);

int((1/2)*(_t0*Heaviside(_t0-1)-_t0*Heaviside(_t0)-Heaviside(1-_t0)*Heaviside(-_t0)+Heaviside(-_t0)+Heaviside(1-_t0)-1)*(1+erf((1/2)*_t0*2^(1/2)))*(2^(1/2)*Heaviside(_t0-1)*exp(-(1/2)*_t0^2)*_t0-2^(1/2)*Heaviside(_t0)*exp(-(1/2)*_t0^2)*_t0-2^(1/2)*Heaviside(-_t0)*Heaviside(1-_t0)*exp(-(1/2)*_t0^2)-Pi^(1/2)*undefined*erf((1/2)*_t0*2^(1/2))*Dirac(_t0)-Pi^(1/2)*undefined*erf((1/2)*_t0*2^(1/2))*Dirac(_t0-1)+2^(1/2)*Heaviside(-_t0)*exp(-(1/2)*_t0^2)+2^(1/2)*Heaviside(1-_t0)*exp(-(1/2)*_t0^2)-Pi^(1/2)*undefined*Dirac(_t0)-Pi^(1/2)*undefined*Dirac(_t0-1)+Pi^(1/2)*Heaviside(_t0-1)*erf((1/2)*_t0*2^(1/2))-Pi^(1/2)*Heaviside(_t0)*erf((1/2)*_t0*2^(1/2))-exp(-(1/2)*_t0^2)*2^(1/2)+Pi^(1/2)*Heaviside(_t0-1)-Pi^(1/2)*Heaviside(_t0))/Pi^(1/2), _t0 = -infinity .. t)

whereas Mma 11 produces the correct piecewise expression (see that here screen15.11.16.docx).

Edit. Mma output.

I'm trying to write some simple code to help plot an approximation of a solution to the heat equation. On one set of axes I want to plot the temperature characteristic of a 1-D bar at different temperatures. I am able to do this easily enough, but the code uses a for loop and my core i7 is not being utilized properly. To rectify this I have tried to use threads, but after many hours of playing around I cannot get it to work fully. The first time I run my code I am faced with the error "Error, (in plot) cannot determine if this expression is true or false: not fHwLibraryInitialized".

Then after this I can run my code a few more times and it will produce a graph, but each time the graph will be different, meaning most of them are incorrect. After running my display command several times I am then faced with the error "kernel connection has been lost", and I need to restary. I have never played with threads before and have very little understanding of how it works.

with(Threads);
with(plots);

para := proc (A, B)
return [op(A), op(B)]; 
 end proc;

HeatP := proc (f, g, n, lambda, X, N, M, T)
 local k, u, i, L, a; 
u := proc (x, t) options operator, arrow; sum((int(sin(k*lambda*x)*g, x = 0 .. 1))*exp(-k*lambda*t)*sin(k*lambda*x), k = 1 .. n) end proc;
 L := [];
 if N-X < 3 then for i from X to N do
 L := [op(L), plot(u(x, (i*T-T)/M), x = 0 .. 1)] 
end do; 
return L;
else a := floor((1/2)*(N-X))+X; 
Threads:-Task:-Continue(para, Task = [HeatP, f, g, n, lambda, X, a, M, T], Task = [HeatP, f, g, n, lambda, a+1, N, M, T])
 end if
 end proc;

g := piecewise(x <= 0, 0, 0 < x and x < 1, exp(-2/(1-(x-1/2)^2)), x >= 0, 0);

display(Threads:-Task:-Start(HeatP, 0, g, 5, Pi, 1, 8, 8, 2));

I cannot see where the problem is, but there is obviously at least one. Any help would be appreciated.

What method of worksheet dissemination maximizes the number of recipients and methods of access?

Our students do have access to Maple on campus. However some students beyond our campus do not. Hence I'm trying maximize the number of recipients. This means that the recipient can accesses my worksheet either via tradition computer OSs such as Mac, Linux, Windows, or via tablet / personal computer OSs: iOS or Android. 

For a static worksheet, I've found PDF is the obvious choice. However, I'm writing more modifiable worksheets. The most obvious options I see are:

* MaplePlayer
* MapleCloud
* Maple WWW, by DigiArea Team

In investigating the options, I have found that:
* Maple Player is useful to those who have no access to Maple, but there appears to be no Android/iOS version of it.
* MapleCloud looks very promising since all OSs read the web, but there appears to be no way a person can join a private group via the MapleCloud website. Am I missing or misreading the help?
* Maple WWW, while pretty, appears to be painfully slow on a tablet (at least the examples I looked at).

Thus the initial purpose of this post is to answer the questions I raise, which I view as far more strategic, with more nuances than most questions that appear on this forum. And the second purpose is to initiate more discussion among the users year, particularly from the educational community, about guiding the direction of Maple.

I look forward to your comments and suggestions.

Let us consider 

J := int(x^n/sqrt(1+x^n), x = 0 .. 1) assuming n > 0;

2*(2^(1/2)-hypergeom([1/2, 1/n], [(n+1)/n], -1))/(2+n)

limit(J,n=infinity);
FAIL
MultiSeries:-limit(J,n=infinity);
FAIL

Mma 11 finds the limit is zero. Hope one feels the difference.

how i

Find the length of the curve in interval of x

Im trying to solve 12 equations with 12 variables but I can't solve. Please help and advise me to solve this problem. Iproject3.mw
project3.mw

 

 

 

First 1031 1032 1033 1034 1035 1036 1037 Last Page 1033 of 2224