Dear all,

I would like to ask you to help me with the following problem that I got error.

Thanks in advance

Is it possible to animate plot like this parametically:
animate(plot, [[cos(t), sin(t), t = 0 .. A]], A = 0 .. 2*Pi, scaling = constrained, frames = 50)

however given the x and y components as solutions of an implicit equation.

I know I could run RootOf. But it seems that there is a part missing due to a branch cut :-/

I had problems plotting the solutions of my PDE with the derivative as the axis. Anyone knows what went wrong?

p1 := pds:-plot(diff(f(x,y),y), y = 0 .. 10, x = 0, numpoints = 100)

I seem to be getting a plot of y=0. Any help would be greatly appreciated:))

fyp2.mw

I have a system of pde as follow,

PDE := [(x*y+1)*(diff(f(x, y), y, y, y))+(x+(3/4)*f(x, y))*(diff(f(x, y), y, y))-(1/2)*(diff(f(x, y), y))^2+T(x, y) = (1/4)*x*(diff(f(x, y), y))*(diff(diff(f(x, y), y), x))-(1/4)*x*(diff(f(x, y), x))*(diff(f(x, y), y, y)), (x*y+1)*(diff(T(x, y), y, y))/(.733)+(x/(.733)+(3/4)*f(x, y))*(diff(T(x, y), y)) = (1/4)*x*(diff(f(x, y), y))*(diff(T(x, y), x))-(1/4)*x*(diff(f(x, y), x))*(diff(T(x, y), y))]

sys_ode := diff(g(y), y, y, y)+(3/4)*g(y)*(diff(g(y), y, y))-(1/2)*(diff(g(y), y))^2+h(y) = 0, (diff(h(y), y, y))/(.733)+(3/4)*g(y)*(diff(h(y), y)) = 0

ics := g(0) = 0, h(0) = 1, (D(g))(10) = 0, g(10) = 0, h(10) = 0

sol2 := dsolve([sys_ode, ics], numeric)

BC := {T(0, y) = h(y), T(x, 0) = 1, T(x, 10) = 0, f(0, y) = g(y), f(x, 0) = 0, f(x, 10) = 0, (D[2](f))(x, 0) = 0}

pds := pdsolve(PDE, BC, numeric)

module() ... end module

pds:-plot(T, y = 0 .. 10, x = 0);

Error, (in pdsolve/numeric/plot) unable to compute solution for x<HFloat(0.0):
solution becomes undefined, problem may be ill posed or method may be ill suited to solution

When I try to use the solution of the Ode as the boundary condition for PDE, by subbing g(y) and h(y) into BC. The plot returns me the error. Anyone knows the reason behind this and how to solve? Any help would be really greatly appreciated. Thanks

hi.

i have problem with this eq.

i need your help.

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

Am here again, pls help me check out this adm code, cant get a result.

below is the attached file

 adm_2_method.mw

Consider the following code snippet:

with(DifferentialGeometry):
DGsetup([x],M);
RemoveFrame(M);

By itself it is, of course, pretty meaningless, but that is not the point. My problem is that after Maple has executed it, it will annoyingly continue using M> for each new execution group that is inserted. Having removed the frame with the command RemoveFrame, I would have expected that not to be the case. How can the frame be truly removed?

The keyboard shortcut sequence for executing an entire worksheet is Alt+e+e+w. Occasionally, I hit Alt+e+w+w by accidence. In that case, Maple sometimes, not always!?, gets completely stuck and can only be shut down through the job list window (Windows Ctrl+Alt+Del). Can anybody else replicate that? And if so, has this erratically odd behaviour been fixed in later versions of Maple?

S:=...+a^4*b*c*x^2+...

So I know it is possible to collect for example any letter collect(S,x)

But I want to collect in particular for an entire product in this case

e.g. collect(S,a^4*b*c)

This obviously doesnt work but is it possible to get it working somehow?

There are 4 points which I want to plot as solid spheres in 3d. Here is the code:

Rplot := pointplot3d(`<,>`(r, Transpose(Vector([0, 0, 0]))), style = point, symbolsize = 40, color = [red, grey, grey, green])

plotsetup(ps, plotoutput = `E:\\.../C.eps`, plotoptions = `width=2000, height=2000, noborder`); print(plots:-display(Rplot, axes = boxed, labels = [x, y, z], symbol = solidsphere));

plotsetup(default, plotoptions = `width=2000, height=2000, noborder`); display(Rplot, axes = boxed, labels = [x, y, z], symbol = solidsphere)

In Maple itself everything looks fine, but once I try to export it to ".ps" I see this.

The problem lies in "symbol = solidsphere". When I change it to "diamond", for example, export is successful:

But I need solidsphere. What should I do to fix this error (except manual editing in Adobe Illustrator)?

This may be a silly question, but does there exist some simple way of (Taylor) expanding an expression of 'small' functions in terms of these functions.

A simple example: Assume that diff(f(x),x) and g(x) are two functions both with range, say, in [-a,+a], where a << 1, and consider the following expression:

sqrt(1 + diff(f(x),x)) * (2 + g(x));

Its expansion to first order in terms of diff(f(x),x) and g(x) should be 2 + diff(f(x),x) + g(x). My problem is that mtaylor does not accept functions as variables to expand on, and I would prefer not to have to substitute back and forth with some 'placeholders'.

Can we calculate the following equations in Maple?

Substituting equations (21) and (22) into (17), and then obtain equation (23). How to do that? I have done this, but the results are complex and large. They are not in a sum form, but in an expansion form. The reference and the maple file are attached.

Hope for your help.

Best wishes,

Kang

Dynamic_buckling_of_thin_isotropic_plates_subjected_to_in-plane_impact.pdf

gg.mw

Can Maple 17.0 work on Windows 10 Home Edition? or it can be only Windows 10 Pro Edition

Trying to plot the function 7*x^2+22*xy+7*y^2+14*xz*sqrt(3)+14*yz*sqrt(3)-5*z^2 = 180, I tried using implicitplot3d to plot it, with ranges I'm quite certain should contain the surface but for whatever reason all I get is a blank plot with no graphics. Here's my input:

Loading plots

with(plots);
[animate, animate3d, animatecurve, arrow, changecoords,

complexplot, complexplot3d, conformal, conformal3d,

contourplot, contourplot3d, coordplot, coordplot3d,

densityplot, display, dualaxisplot, fieldplot, fieldplot3d,

gradplot, gradplot3d, implicitplot, implicitplot3d, inequal,

interactive, interactiveparams, intersectplot, listcontplot,

listcontplot3d, listdensityplot, listplot, listplot3d,

loglogplot, logplot, matrixplot, multiple, odeplot, pareto,

plotcompare, pointplot, pointplot3d, polarplot, polygonplot,

polygonplot3d, polyhedra_supported, polyhedraplot, rootlocus,

semilogplot, setcolors, setoptions, setoptions3d, spacecurve,

sparsematrixplot, surfdata, textplot, textplot3d, tubeplot]
implicitplot3d(7*x^2+22*xy+7*y^2+14*xz*sqrt(3)+14*yz*sqrt(3)-5*z^2 = 180, x = -50 .. 50, y = -50 .. 50, z = -50 .. 50, axes = normal);

Not sure what I'm doing wrong exactly.