MaplePrimes Questions

Hello.

Can I ask maple to give me a vector parametrization of any arbituary surface? Lets say I am given some surface in the form z=(g(x,y) or f(x,y,z).

In some cases finding a vector parametrization can be tidous and hard. I have looked at the algcurves help page, it is helpful in som cases but not that much for surfaces.

Lets for the purpose of this exampl say I wanted to obtain a parametrization for the torus (C-sqrt(x^2+y^2))^2+z^2=a^2

Is maple able to help me with this process beside making the plot and start the manual visualization and work?

Hello everyone.

I was wondering if anyone knew of a way to have a supscript+superscript layout on the left hand side of the equation in 1D math?

When I try to use (A)[b]^(c) from the layout palette it does not allow me to use it in 1D Math, and when using it in any of the math types (1D and both of the 2D options) it is not accepted on the left hand side of the equation.

I am writing a project where we use transformation matrices, and I seem to be unable to define them as they are supposed to be. I am left with T__40 when I want to write (T)[4)^(0).

Help would be very much appreciated.

I  encountered a non-integrable integral in the process of solving the following process, . How to achieve its numerical solution? Such as in a looping   code:

#######
pa[i] := pa[i-1]-(Int(subs(t = tau, Lpa[i-1]+Na1[i-1]-Na2[i-1]), tau = 0 .. t)); 

pw[i] := pw[i-1]-(Int(subs(t = tau, Lpw[i-1]+Nw1[i-1]-Nw2[i-1]), tau = 0 .. t)); u[i] := u[i-1]-(Int(subs(t = tau, Lu[i-1]+Nu1[i-1]+Nu2[i-1]), tau = 0 .. t));

######
Detailed code see annexBC2.mw

what does it means and what it will do. Can some one help me for solving this

Shootlib := "C:/Shoot9"; libname := Shootlib, libname; with(Shoot);

while i m receiving the following message:

"Error, invalid input: with expects its 1st argument, pname, to be of type {`module`, package}, but received Shoot "

Full program is :

restart; Shootlib := "C:/Shoot9"; libname := Shootlib, libname; with(Shoot);
with(plots):
N1 := 1.0; N2 := 2.0; N3 := .5; Bt := 6; Re_m := N1*Bt; gamma1 := 1;
FNS := {f(eta), fp(eta), fpp(eta), g(eta), gp(eta), m(eta), mp(eta), n(eta), np(eta), fppp(eta)};
ODE := {diff(f(eta), eta) = fp(eta),
        diff(fp(eta), eta) = fpp(eta),
        diff(fpp(eta), eta) = fppp(eta),
        diff(g(eta), eta) = gp(eta),
        diff(gp(eta), eta) = N1*(2.*g(eta)+(eta-2.*f(eta))*gp(eta)+2.*g(eta)*fp(eta)+2.*N2*N3*(m(eta)*np(eta)-n(eta)*mp(eta))),
        diff(m(eta), eta) = mp(eta),
        diff(mp(eta), eta) = Re_m*(m(eta)+(eta-2.*f(eta))*mp(eta)+2.*m(eta)*fp(eta)),
        diff(n(eta), eta) = np(eta),
        diff(np(eta), eta) = Re_m*(2.*n(eta)+(eta-2.*f(eta))*np(eta)+2.*N2/N3*m(eta)*gp(eta)),
        diff(fppp(eta), eta) = N1*(3.*fpp(eta)+(eta-2.*f(eta))*fppp(eta)-2.*N2*N2*m(eta)*(diff(mp(eta), eta)))
       }:
   
blt := 1.0;
IC := { f(0) = 0,
        fp(0) = 0,
        fpp(0) = alpha1,
        g(0) = 1,
        gp(0) = beta1,
        m(0) = 0,
        mp(0) = beta2,
        n(0) = 0,
        np(0) = beta3,
        fppp(0) = alpha2
      };
BC := { f(blt) = .5,
        fp(blt) = 0,
        g(blt) = 0,
        m(blt) = 1,
        n(blt) = 1};
infolevel[shoot] := 1;
 

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

 

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)?

how to find the first 6 non-zero terms of MacLaurin series of the function erf(x)

how to Prove that the circumference of a circle of radius r is 2πr
on maple ?????

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.

First 1055 1056 1057 1058 1059 1060 1061 Last Page 1057 of 2434