Unanswered Questions

This page lists MaplePrimes questions that have not yet received an answer

Hi I have encounted this error:

Error, (in diffalg/DAparse_polynomial) unknown indeterminate or coefficient outside base field, -Cp+Cf/(1+Kw*(-sigma*`ΔPi`+`ΔP`)*(1-(1/2)*Fp/Ff)/Ks)

The text highlighted in bold is an equation that I typed in.

 

Can I know what it means?

1,1,2,2,2,3,3,2,2,1,2,2,2,3,3,1,1

i would like to research the ordering mechanism of this kind of series Before sorting.

How I can simplify result? For factor or using rule.

Thanks

 

can maplesim simulate a robot arm and hand catch a throwing ball and picking objects?

can maplesim simulate a balancing a robot legs?

will these simulation require complex formula to draw the robots arm and hand and maintain the joints stick together during motions?

can maplesim output the specification to autodesk for 3D print parts?

how much of a maplesim? 

i can not see the price in maple official site from Hong Kong area.

HI, 

Yesterday I replied  denitsastaicova's question while providing an attached file. 

Surprisingly, I wasn't able to save the worksheet after having executed it. The only way I found to save it was to copy-paste it's content in a new worksheet and save this latter without any execution.
Here is the file
events.mw

For info I'm working with  Standard Worksheet Interface, Maple 2015.2, Mac OS X, December 21 2015 Build ID 1097895
on Mac OS Mojave 10.14.3
If I open this file and execute it (and even if I execute only the commands from the initial restart to the three dsolve), I get this pop-up error message when trying to save it:

Does my worksheet contain some undesired character or is it a bug ?

I have this very ugly ODE system which I solve with "dsolve, numeric, complex=true". It works ok but I'd really like to be able to use events. The problem is that the solution of the system has very small (~10^(-Digits)) imaginary part and as we know "events" works only for real ODEs.

My question is does anyone know a way to force dsolve to solve the ODE as a real system? In Mathematica, there is the function "chop" which does the trick. In Maple, I tried with fnormal and also by forcing some terms to be real but that makes it hit a singularity at 0.01 which doesn't exist when it integrates it as a complex system. In Mathematica, it integrates ok with "chop" without any random singularities. 

I'm currently trying separating the system into real and imaginary parts, but considering how ugly the system is, it may not be so simple (and so far, it doesn't work). So something to trick dsolve would come handy.

Thanks in advance.

How can plot two  and bipolar figures as attached and extract their data.

Thank You.

I find many tools about differential 

i had some differential ideal

how to use these tools to research what application can be applied with differential ideals I discovered?

If A is an m × n matrix and B is a p × q matrix, then the Kronecker product C = A ⊗ B is the mp × nq block matrix. Assume I have matrix C and want to find matrix A and B. This problem has known solution called "Nearest Kronecker Product". So I just need a function like this: A,B:=NearestKroneckerProduct(C)  which minimizes ||C - A ⊗ B||F  where denotes frobenius norm. 

Here is an Image that shows a mathematica code:

I took it from the link below:

https://mathematica.stackexchange.com/questions/91651/nearest-kronecker-product

Also there is a matlab code:

https://gist.github.com/tholden/58dd9a8991daa750ae36a633fe7060a4/revisions

can you convert mathematica code to maple?

If I do:

df:=DataFrame(Matrix(3,4,[seq(1..12)]), rows=[a,b,c],columns=[A,B,C,D]);Tabulate(df, width=100)

 

The font that Maple uses for the Tablulate is much larger than the font used to display the Dataframe. How does one choose the font size that Tabluate() uses? 

Peter

We put A = {1,2,3,4} and B is the set of lower cubes
to 100 of positive integers.Determine A ∪ B  A ∩ B , A ∩ {6,8} and A \ B 

Hi,

i am trying to solve a differential equation numercially but since it is second degree,I can not achieve the proper answer.What should I do?

thank you
diff1.mw
 

restart

h := 1-.8*x+.5*(x^2-x)

1-1.3*x+.5*x^2

(1)

z := 3^(3/2)*(((2+k*De^2*(diff(p(x), x))^2*h^2)^3*h^2+108*k*De^2)/k)^(1/2)

3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2)

(2)

c := (1/6)*((-54*De+z)*k^2*h^2)^(1/3)/(k*De*h)-(1/2)*h*(2+k*De^2*(diff(p(x), x))^2*h^2)/(De*((-54*De+z)*k^2*h^2)^(1/3))-(1/2)*h*(diff(p(x), x))

(1/6)*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3)/(k*De*(1-1.3*x+.5*x^2))-(1/2)*(1-1.3*x+.5*x^2)*(2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)/(De*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3))-(1/2)*(1-1.3*x+.5*x^2)*(diff(p(x), x))

(3)

ode1 := ((1/2)*c*(diff(p(x), x))^2*h^4+(diff(p(x), x))*c^2*h^3+(1/10)*(diff(p(x), x))*h^5+c^3*h^2)*k*De^2+(1/2)*c*h^2+(1/6)*(diff(p(x), x))*h^3+h = 0

((1/2)*((1/6)*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3)/(k*De*(1-1.3*x+.5*x^2))-(1/2)*(1-1.3*x+.5*x^2)*(2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)/(De*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3))-(1/2)*(1-1.3*x+.5*x^2)*(diff(p(x), x)))*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^4+(diff(p(x), x))*((1/6)*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3)/(k*De*(1-1.3*x+.5*x^2))-(1/2)*(1-1.3*x+.5*x^2)*(2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)/(De*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3))-(1/2)*(1-1.3*x+.5*x^2)*(diff(p(x), x)))^2*(1-1.3*x+.5*x^2)^3+(1/10)*(diff(p(x), x))*(1-1.3*x+.5*x^2)^5+((1/6)*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3)/(k*De*(1-1.3*x+.5*x^2))-(1/2)*(1-1.3*x+.5*x^2)*(2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)/(De*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3))-(1/2)*(1-1.3*x+.5*x^2)*(diff(p(x), x)))^3*(1-1.3*x+.5*x^2)^2)*k*De^2+(1/2)*((1/6)*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3)/(k*De*(1-1.3*x+.5*x^2))-(1/2)*(1-1.3*x+.5*x^2)*(2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)/(De*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3))-(1/2)*(1-1.3*x+.5*x^2)*(diff(p(x), x)))*(1-1.3*x+.5*x^2)^2+(1/6)*(diff(p(x), x))*(1-1.3*x+.5*x^2)^3+1-1.3*x+.5*x^2 = 0

(4)

ivp := {ode1, p(0) = 0}

{((1/2)*((1/6)*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3)/(k*De*(1-1.3*x+.5*x^2))-(1/2)*(1-1.3*x+.5*x^2)*(2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)/(De*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3))-(1/2)*(1-1.3*x+.5*x^2)*(diff(p(x), x)))*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^4+(diff(p(x), x))*((1/6)*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3)/(k*De*(1-1.3*x+.5*x^2))-(1/2)*(1-1.3*x+.5*x^2)*(2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)/(De*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3))-(1/2)*(1-1.3*x+.5*x^2)*(diff(p(x), x)))^2*(1-1.3*x+.5*x^2)^3+(1/10)*(diff(p(x), x))*(1-1.3*x+.5*x^2)^5+((1/6)*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3)/(k*De*(1-1.3*x+.5*x^2))-(1/2)*(1-1.3*x+.5*x^2)*(2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)/(De*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3))-(1/2)*(1-1.3*x+.5*x^2)*(diff(p(x), x)))^3*(1-1.3*x+.5*x^2)^2)*k*De^2+(1/2)*((1/6)*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3)/(k*De*(1-1.3*x+.5*x^2))-(1/2)*(1-1.3*x+.5*x^2)*(2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)/(De*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3))-(1/2)*(1-1.3*x+.5*x^2)*(diff(p(x), x)))*(1-1.3*x+.5*x^2)^2+(1/6)*(diff(p(x), x))*(1-1.3*x+.5*x^2)^3+1-1.3*x+.5*x^2 = 0, p(0) = 0}

(5)

dsolve(ivp, p(x), numeric, parameters = [k, De])

Error, (in dsolve/numeric/make_proc) Could not convert to an explicit first order system due to 'RootOf'

 

``


 

Download diff1.mw

 

 


 

Download diff1.mw

Hello,

Is it possible to expand vector calculus identities directly in Maple without taking them to basis form (i.e. their component partial derivatives).

For example: 

restart

with(Physics[Vectors]):

%Curl(u*%Gradient(v)) = u*%Curl(%Gradient(v)+`&x`(%Gradient(u), %Gradient(v)))

 

Given the Left Hand Side, can Maple come up with the RHS on its own?

Many Thanks.

Please, I need to use Maple to solve Euler-Bernoulli Beam on Pasternak Foundation using Homotopy Perturbation Method.

The governing equation is 

initial conditions are 

the boundary condition is 

The governing equation represents Euler-Bernoulli beam on a generalized Pasternak viscoelastic foundation under an
arbitrary distributed dynamic load. in which E, I , ρ ,A are the parameters of the beam, representing Young’s modulus of elasticity, moment of inertia, density and area of cross section, respectively. K,C and Gp are spring stiffness, damping coefficient, and shear coefficient of the foundation. Moreover, y(x, t) and F(x, t) are defined as the vertical deflection of the beam and the generic arbitrary dynamic loads, respectively, where the loads distribute along the x-axis and t is time.

 

I will appreciate anyone who can help me with a Maple solution.

 

Thank you.

How should wirte a while loop for solve nonlinear equations by newton raphson method

First 92 93 94 95 96 97 98 Last Page 94 of 334