Hello Everyone;

Hope you are fine. I I have make a code but there is error and dsolve comand is working. I also need to plot. Kindly help me. I am waiting kind response. Code is attached.

Download Question1.mw

Thanks

I would like to collect in derivatives of a function(s) in a worksheet where I have the Physics package loaded but It does not want to do it. When I bring the expression to a blank worksheet it does everything just fine. Am I missing something small? Or is there additional commands/options I need to input? 

Thoughts?

CollectDiff.mw

The highly oscillatory integrand is: 

(sin(x + sqrt(x)) + x*BesselJ(0, x^2))/(1 + x): # int(%, x = 0 .. infinity, numeric); # Note that the upper limit of `x` is NOT finite.

Unfortunately, Maple still cannot return a result in a few minutes. For instance, 

Download oscillatoryInt.mws

I tried to specify another method (as described in evalf/int), but it seems that this doesn't work. Any ideas?

For example, if I want to replace all "x" in t:={x1,x2,x3...x100} to get t:={y1,y2,y3...y100}, how can maple do this, please? Thanks a lot! 

The answers from the book dont make any sense after using their solution. Its not producing the expected results. 

This is the translation by the way: "A battery that has a voltage U, is in series with a resistor R, and a condenser with capacity C. The current strength on a given time "t" in this circuit is given by i=i(t)=U/R*e^-(t/R*C). Determine R and C based on the data given in table 5.6, the current is U=100V."

Greetings,

The Function

Download Mapleprimes_Question_Book_2_Paragraph_5.12_Question_5.mw

Good day! 

I have been wrestling with what appears to be a simple problem, but with no success and so, I thought I would reach out for support. 

Consider a two-dimensional strip of any given width, x, and depth, y.

I wish to maximize the quantity of strips that can fit inside a larger two-dimensional surface of any given width, X, and depth, Y. Overhanging is not permitted and all smaller strips must lie within the interior of the larger surface. 

I'm looking for a solution method that can calculate the best configuration necessary to fit the maximum number of smaller strips in the interior of the larger rectangle. Can someone suggest a way (if possible) to solve for this?

Thanks for reading.

diff(f(x), x) = a-a/((1+a*x)*(1/(1+a*x))^a)-a*f(x) where a is a constant and an initial condition f(0)= b is known Thank you.

Hello,

How do I isolate the variable ng in the following expression:

rho := `ρg`*ng + `ρp`*np + `ρw`*nw;
B := np/Kp + nw/Kw + ng/Kg;

C := (B/rho)^(1/2);

As you can see, I am running the isolate command, but it is not including the variable C, or it returns the wrong result.

isolate(((np/Kp + nw/Kw + ng/Kg)/(`ρg`*ng + `ρp`*np + `ρw`*nw))^(1/2), ng);
                               ( Kp nw + Kw np) Kg
                   ng = -     -----------------

                                          Kw Kp       

isolate(C=((np/Kp + nw/Kw + ng/Kg)/(`ρg`*ng + `ρp`*np + `ρw`*nw))^(1/2), ng);
                             0 = 0

Thanks.                

I think some form of simplify() would do but I am not sure how.

See the following script for more details:

Download rearrangingterms.mw

It must be something similar to this https://www.mapleprimes.com/maplesoftblog/201455-Rearranging-The-expression-Of-Equations, but my case is slighty different.

I was not able to apply ICs (initial conditions) using the advanced parameter settings of the flexible beam component in the attached file. My intention is to have a beam moving at t=0s (velocity of end frame <> 0). The computed initial values for the elastic coordinates (see output console for the disabled component below encircled in red) indicate that the flexible beam has no kinetic energy at t=0s.

A workaround is to attach a rigid body to the end frame and apply the ICs to the body. (Inertia and mass of the rigid body can be set to zero.)

How to apply initial conditions to the disabled component so that it moves as the workaround?

Imposing_IC_v.msim

In certain tasks, I need to find all accurate positive (not just nonnegative) roots that exist of some multivariate polynomials like: 

nsd := 16*a*b*c*(9 + a^2 + b^2 + c^2)*(b*c + a*(b + c) + 3*(a + b + c)) - (3 + a + b + c)^2*(a*b + 3*c)*(3*b + a*c)*(3*a + b*c): # assume((a, b, c) >~ 0);

According to fsolve/details, for one general equation, the fsolve command only computes "a single real root", so it is inadequate to tackle this question. But if I use the solve command with floating-point values, the computation cannot finish in ten minutes instead! (Maybe a longer time will suffice, yet this is rather unacceptable.) 

Download solve_numerically.mws

Is there any workaround to obtain those (finitely many) positive solutions completely?

Dear all:

    I am using maple2022 to learn Fourier Series package - Maple Application Center (maplesoft.com) . the author provide the package here:

cgi.math.muni.cz/kriz/xsrot/fourierseries/  I use the lastest version: Package FourierTrigSeries version 0.41 

there is no .mla so I do not know how to use this package.(the package is a little old). 

please enlight me if you know how to install this old package.

but the author provide the code inside fourierseries-structs.mws 

here is what I do:

1. first delete some lines which is not useful(which I shared in the attachment as FourierTrigSeries.mw)

and remain everything inside module and save as .mw (maybe the orginal .mws is a old version of maple files)

2. export as .FourierTrigSeries.mpl (also attached)

3. use the examples.mws (the code also from author)

I got error as:

the package can be recognized, but the first line wil give some error which is out of my knowledge field.

everything mentioned above is here, including the original author's files.

fouriertrigseries_0_41en_mapleprime.zip

Could you please have a look. your idea is valuable to me.

PS: this package is not available The Fourier Series package for Maple - Maple Application Center (maplesoft.com) whenI try to learn  Teaching Fourier Series with Maple II 

if anyone keep it before, please share it.

thanks for your help.

rockyicer 

restart;  
with(geometry):  
with(plots):  
_EnvHorizontalName = 'x':  _EnvVerticalName = 'y':
point(A, -1, 9):                                                                                                       
point(B, -5, 0):
point(C, 6, 0):
triangle(ABC,[A,B,C]):
midpoint(M1,A,C): midpoint(M2,B,C):midpoint(M3,A,B):
rotation(J, C, Pi/2, 'counterclockwise', M1):triangle(AJC,[A,J,C]):
rotation(Ii, C, Pi/2, 'counterclockwise', M2):triangle(BIC,[B,Ii,C]):
rotation(K, A, Pi/2, 'counterclockwise', M3):triangle(AKB,[A,K,B]):
midpoint(O1,K,J): coordinates(O1):
midpoint(O2,A,Ii): coordinates(O2):  
poly:=[coordinates(A),coordinates(J),coordinates(Ii),coordinates(K)]:   

display(draw([A(color = black, symbol = solidcircle, symbolsize = 12), 
B(color = black, symbol = solidcircle, symbolsize = 12), 
C(color = black, symbol = solidcircle, symbolsize = 12), 
J(color = black, symbol = solidcircle, symbolsize = 12), 
polygonplot(poly,color = "DarkGreen", transparency = 0.5),
ABC(color = red ),
BIC(color = green),
AKB(color = grey),
AJC(color =blue)]),
textplot([[coordinates(A)[], "A"],[coordinates(J)[], "J"],[coordinates(Ii)[], "I"],   
[coordinates(B)[], "B"], [coordinates(K)[], "K"], 
[coordinates(C)[], "C"]], 
align = [above, right]),  axes = none);
Error, (in geometry:-draw) the option must be of type equation or name. I don't see how to correct this error/
 

Hi, I have the equation below and I'm trying to use "collect " to get the coefficient of dudy^2 and dvdx^2. I tried

aa1 := subs({gamma[1] = alpha[3]-alpha[2], gamma[2] = alpha[6]-alpha[5]}, collect(aa, {dudy^2, dvdx^2}))

But is not collecting the coefficient of dvdx^2. Please can someone help?

aa := alpha[1]*(cos(theta(x, y, t))^4*dudx^2+2*cos(theta(x, y, t))^3*dudx*sin(theta(x, y, t))*dvdx+2*cos(theta(x, y, t))^3*dudx*sin(theta(x, y, t))*dudy+2*cos(theta(x, y, t))^2*dudx*sin(theta(x, y, t))^2*dvdy+cos(theta(x, y, t))^2*sin(theta(x, y, t))^2*dvdx^2+2*cos(theta(x, y, t))^2*sin(theta(x, y, t))^2*dvdx*dudy+2*cos(theta(x, y, t))*sin(theta(x, y, t))^3*dvdx*dvdy+cos(theta(x, y, t))^2*sin(theta(x, y, t))^2*dudy^2+2*cos(theta(x, y, t))*sin(theta(x, y, t))^3*dudy*dvdy+sin(theta(x, y, t))^4*dvdy^2)+(alpha[2]+alpha[3]+gamma[2])*(-dudx*cos(theta(x, y, t))*sin(theta(x, y, t))*thetadot-(1/2)*dudx*cos(theta(x, y, t))*sin(theta(x, y, t))*dudy+(1/2)*dudx*cos(theta(x, y, t))*sin(theta(x, y, t))*dvdx+(1/2)*cos(theta(x, y, t))^2*dvdx*thetadot-(1/4)*cos(theta(x, y, t))^2*dvdx^2+(1/2)*cos(theta(x, y, t))^2*dudy*thetadot+(1/4)*cos(theta(x, y, t))^2*dudy^2-(1/2)*sin(theta(x, y, t))^2*dvdx*thetadot+(1/4)*sin(theta(x, y, t))^2*dvdx^2-(1/2)*sin(theta(x, y, t))^2*dudy*thetadot-(1/4)*sin(theta(x, y, t))^2*dudy^2+dvdy*sin(theta(x, y, t))*cos(theta(x, y, t))*thetadot-(1/2)*dvdy*sin(theta(x, y, t))*cos(theta(x, y, t))*dvdx+(1/2)*dvdy*sin(theta(x, y, t))*cos(theta(x, y, t))*dudy)+alpha[4]*(dudx^2+(1/2)*dvdx^2+dvdx*dudy+(1/2)*dudy^2+dvdy^2)+(alpha[5]+alpha[6])*(cos(theta(x, y, t))^2*dudx^2+2*cos(theta(x, y, t))*dudx*((1/2)*dvdx+(1/2)*dudy)*sin(theta(x, y, t))+cos(theta(x, y, t))^2*((1/2)*dvdx+(1/2)*dudy)^2+2*cos(theta(x, y, t))*((1/2)*dvdx+(1/2)*dudy)*dvdy*sin(theta(x, y, t))+sin(theta(x, y, t))^2*((1/2)*dvdx+(1/2)*dudy)^2+sin(theta(x, y, t))^2*dvdy^2)+gamma[1]*(sin(theta(x, y, t))^2*thetadot^2+sin(theta(x, y, t))^2*dudy*thetadot-sin(theta(x, y, t))^2*dvdx*thetadot+(1/4)*sin(theta(x, y, t))^2*dudy^2-(1/2)*sin(theta(x, y, t))^2*dudy*dvdx+(1/4)*sin(theta(x, y, t))^2*dvdx^2+cos(theta(x, y, t))^2*thetadot^2-cos(theta(x, y, t))^2*dvdx*thetadot+cos(theta(x, y, t))^2*dudy*thetadot+(1/4)*cos(theta(x, y, t))^2*dvdx^2-(1/2)*cos(theta(x, y, t))^2*dudy*dvdx+(1/4)*cos(theta(x, y, t))^2*dudy^2)+xi*(cos(theta(x, y, t))^2*dudx+2*cos(theta(x, y, t))*sin(theta(x, y, t))*((1/2)*dvdx+(1/2)*dudy)+sin(theta(x, y, t))^2*dvdy);

hi every body 

please How to write these solutions in fraction form

Download problem1.mw