MaplePrimes Questions

Search Questions:

Latest Questions Latest Questions Feed

I need to solve an ode of the type ay''+by'+cy=f(x) using cubic b spline.

can any one help me with the code or algorithm. Thanks

N := 4;
print(`output redirected...`); # input placeholder
4
y := sum(A[2*n].cos(2.*n.x), n = 0 .. N);

eq1 := diff(y, `$`(x, 2))+(a+2*q*cos(2*x))*y

eq2 := map(combine, eq1, trig)

for i from 0 to 4 do eq4[i] := coeff(eq2, cos(2*n*x)) end do

From these I want to extract the co-ffficients of cos(0x),cos(2x),cos(4x)..

and form a simultaneous linear equation containg A0,A2,A4

The solution is 

aA0+qA2=0

2q*A0+(a-4)*A2+q*A4=0

Can anybody tell me how to do it

so I'm trying this:

restart;

sigma := 0.143e-18;

l_0 := 1.87;

l0 := 1.87;

roll := rand(0 .. 25.0);

f_gauss := proc (x) options operator, arrow; exp(-(1/2)*x^2/`σ_x`^2)/sqrt(2*Pi*`σ_x`^2) end proc;

f_norm := proc (dx) options operator, arrow; int(f_gauss(x), x = -(1/2)*dx .. (1/2)*dx) end proc;

sol_gauss := proc (mix) options operator, arrow; evalf(eval(-ln((int(f_gauss(x)*exp(-2*sigma*N2O*sqrt((1/4)*l_0^2-x^2)), x = -(1/2)*dx .. (1/2)*dx))/f_norm(dx))/(sigma*N2O), [N2O = 0.25e20*mix/100])) end proc;

for ii to 10 do

a := roll();

eval(sol_gauss(a), [dx = l_0, `σ_x` = l0])

end do

Hello people in mapleprimes,

I want to ask you about how to make a function of function which makes a logarithmic derivative of a function.

For example, x^3 is mapped as h(x^3)=3h(x), h(x+y)=(x/(x+y))*h(x)+(y/(x+y))*h(y),

h(g(x)*k(x))=h(g(x))+h(k(x)).

I hope someone give me a hing to create h.

 

Best wishes.

taro

 

 

 

 


environment
- Windows 10 LTSB 64-bit
- Intel Haswell/Skylake

sequence:
1. install Windows unattended
2. install scheduletask for software installation (runs with systemaccount on windows starts and available network)
3. task runs a cmdfile(1) that map a networkshare and runs another cmdfile(2) on the networkshare
4. cmdfile(2) runs softwareinstallation one by one, every install is a seperate cmdfile
5. cmfile(maple.cmd) starts maple installation with the following cmdline
-->
START "Maple 2016" /WAIT "N:\Math\Maple\Maple2016Windowsx86Installer.exe" --mode unattended --desktopshortcut 0 --shortcutname "Maple 2016" --configureMatlab 0 --matlabpath "" --licenseType network --serverName "licserver.company.domain" --portnumber 27000
<--

Note: licenseserver and port are configured correctly, tested with manual installation of Maple

effects:
- installer starts and install some files (found unter C:\ProgramFiles\Maple 2016)
- installer persists in memory 
- no startmenuentry for maple found
- no installler logfile in maple directory found
- install not finished

- try it with --optionfile ... -> hangs
- try copying all files to local harddisk and run it from there -> hangs
- try it running with a local administrator (using psexec) -> hangs

- if im logged in and run the installer cmdline manually with adminconsole or with systemaccount (psexec -i -s -d cmd) the installation runs perfect, but this is not what i need to install hundreds of machines

Dear Community,

I have a complicated function with two independent variables, namely QD( tD,  rD ). I would like to explore its behavior in a semilog plot in the range of tD = 1.0 .. 1000.0 with changing rDs, where I would like to change the values of rD with a slider of the Explore command. I've developed the attached code, yet it gives me an error message, "unexpected option... " Could you pls. have a look at it what do I do wrong?

Tx for the kind help in advance

best regards

Andras

 Explore_VEHQD.mw

Combining graphs...

September 22 2016 montse 5

Hi

Can anyone tell me how to combine several graphs into one? Fo instance I have one graph with x going from 0-1,another with x going from 1-20 and yet another with x going from 20-21. I can plot them individually but how do I put all three together?

I'm sorry I'm not doing a very good job explaining, but I hope it is understandable.

 

 

 

Hello

Thanks to the help of several list members I managed to translate and update some old procedures.  These procedures were part of a packcage and they all had help files.   Here is an example

 

# E. Mendes - 25/04/94

`help/text/sampling`:=TEXT(``,
`FUNCTION: sampling - finds the discretized model`,
` `,
`CALLING SEQUENCE:`,
` sampling(f)`,
` sampling(f,k)`,
` sampling(f,k,vars)`,
`PARAMETERS:`,
` f - state space`,
` k - order of approximation`,
` vars - variables`,
``,
`SYNOPSIS:`,
`- sampling(f,k) returns a discrete approximation for a continuous`,
` system`,``,
` dx(t)/dt=f(x(t),u(t)) `,
``,
`EXAMPLE: `,``,
`> with(linalg):`,
`> f:=vector([x2,x3,-x1+x2^2+2*x3]):`,
`> sampling(f,1)`,
``,
` 2`,
` [- x1, x2 + x1 ]`,``,
`SEE ALSO: fixpoind, fixpoinc`):

 

I have browsed the help documentation but I must confess I don't feel comfortable to modify the help database.  Do I have to wrap all the functions up as a package and then write the help files?   I am lost here.  Any help will be most welcome.

Many thanks

 

Ed

 

Hello,

I would like to silmplify a trigonometric equation with some squares. I'm sure this equation can be reduced but i didn't manage to simplify it with Maple.

Here is the equation (named condition in the maple file) that I would like to simplify:

simplification_condition_de_compatibilité.mw

May you help me to simplify this trigonometric equation?

I would like simplification if time the pattern cos()^2+sin()^2 appears.

Thanks a lot for your help

Pick any point P and let Q = (1, 1, 1). [Your choice for P can be anything other than the origin or

Find an equation for the line l1 that passes through P and the origin. Plot the line segment formed by P and the origin in Maple 

How would I plot this?


Please, I need assistance with this problem.

Here is the problem I am trying to solve:

restart:
with(plots):
with(LinearAlgebra):
with(PDEtools):
with(Student):

myPDE1 := D11*diff(w(x,y), x$4) + 2*(D12+2*D66)*diff(w(x,y), y$4) + D22*diff(w(x,y), x$2, y$2) - G*diff(w(x,y), x,y)= 0;

pdsolve(myPDE1);

pdsolve(myPDE1, build);

"Boundary conditions";
"(Note:the domain for the problem is a rectangle)";
bc1 := w(0,y) = 0; # @ x=0 edge;
bc2 := w(a,y) = 0;  # @ x=a edge;
bc3 := w(x,0) = 0; # @ y=0 edge;
bc4 := w(x,b) = 0; # @ y=b edge;
bcx1 := -D11*D[2](w)(0,y) - D12*D[2](w)(0,y) = 0; # @ x=0 edge;
bcx2 := -D11*D[2](w)(a,y) - D12*D[2](w)(a,y) = 0; # @ x=a edge;
bcy1 := -D12*D[2](w)(x,0) - D22*D[2](w)(x,0) = 0; # @ y=0 edge;
bcy2 := -D12*D[2](w)(x,b) - D22*D[2](w)(x,b) = 0; # @ y=b edge;

sol := [myPDE1, bc1, bc2, bc3, bc4, bcx1, bcx2, bcy1, bcy2];

pdsolve(sol);

"Note:
and D11, D12, D22, D66 and G are constant.
The intention is to find the critical value for G"

I need help with how I can handle the boundary conditions for the problem. Thanks a million.

Hi,

I want to join two animations, but displayed subsequently, not simultaniously, As far as I got, I have:

plot2 := plots[animate](plot, [tau[1], z = -(1/2)*J .. (1/2)*J, color = red, legend = shear*stress], h = 1 .. 1+1/sqrt(3));
plot3 := plots[animate](plot, [tau[2], z = -(1/2)*h .. (1/2)*h, color = red, legend = shear*stress], h = 0 .. 1);
display([plot2, plot3], insequence = true);

my functions are difined as

tau[1] := h-1+(4/3)*(1/sqrt(3)-h+1)*sqrt(3)*(sqrt(3)*((1/4)*J^2-z^2)/(J^2*b))
tau[2] := (4/3)*sqrt(3)*((1/4)*h^2-z^2)/(h^2*b)

Both plots work individually, however, they have successive ranges. I wondered if I could make 1 plot, with an if statement wether to plot tau[1] or tau[2] depending on where the animation is in the range h. Does anyone know about a solution?

I need  some examples s.t. the computation of their lexicographic Groebner basis is heavy?

Thank you so much.

Already searched and browsed multiple different threads and still cannot find a solution.

Apologizing the noob nature of this question.

 

 

In this code below, Why is the factor command not working?

f := a^2+x^2-2*ax;

a^2+x^2-2*ax

(1)

factor(f);

a^2+x^2-2*ax

(2)

expand((x-a)*(x-a));

a^2-2*a*x+x^2

(3)

``

NULL

NULL


Download factor_polynomial_2_multivariable.mwfactor_polynomial_2_multivariable.mw

 

Hello people in maple primes

I have a question, which is about the matrix shown in http://www.mapleprimes.com/questions/217852-HOW-I-Convert-Root-Of-In-To-Another-Common-Form

Why can't C below be shown with beta?

A := Matrix(3, 3, [[-a, a, 0], [0, 0, -sqrt(l*b*c*(j+k))/(j+k)], [2*j*sqrt(l*b*c*(j+k))/((j+k)*l), 2*k*sqrt(l*b*c*(j+k))/((j+k)*l), -c]]);
B:=subs(l*b*c*(j+k)=alpha,A);
C:=subs(j*alpha^(1/2) = beta,B);
e:=subs(alpha^(1/2) = gamma,B);

Best wishes.

taro

 

 

1 2 3 4 5 6 7 Last Page 1 of 1358