I assigned

before an algebraic calculation so I would like to get  or have the program print the 70 digits of the answer and not just 10 digits. Because when I press ENTER, I get only 10 digits.

I am trying to implement Subresultant p.r.s. algorithm for calculating greatest common divisor. The algorithm decribed in the book:

My code return the correct GCD, however the sub-resultant terms are different from the result of the built-in function. The last term a[i-1] is huge and involves fractions. I think my implementation is same as the algorithm described in the textbook.

I have attached the file. Could anybody spot anything wrong in the code? Why do fractions still appear? In my code, "lsr" is last subresultant term returned from the built-in function, the second one is my result.

Download subresultant.mw

I resolved the coefficients to a 2nd order diff eq of the form:ay''+by'+cy=f(t)

I have included the .mw file for convenience at the link at the bottom of the page.  I resolved the coefficients in 2 different ways & they do not concur.  The 1st approach used the LaPlace transform & partial fraction decomposition.  The coefficient results are given by equations # 14 & 15.  The 2nd approach used undetermined coefficients where I assumed the particular solution and then applied the initial conditions to resolve the coefficients pertaining to the homogeneous solution which are given in the results listed in equation #23.  Noted in the 1st case the coeff's are A₃ & A₄ and for the 2nd approach the coeff's are A₁ & A₂.  I have worked this numerous times & do not understand why they do not concur.  So I thought I should get some fresh eyes on the problem to find where I may have gone wrong.

Any new perspective will be greatly apprecieated.

I had trouble uploading the .mw file so I have included an alternative link to retrieve the file if the code contents is illegible or you cannot dowlad the file drectly from the weblink  Download coeffs_of_homogen_soln_discrepancy.mw.  You should be able to download from the alternative link below once you paste the link into your browser.  If you cannot & wish for me to provide the file in some other fashion respond with some specific instructions & I will attempt to get the file to you.

https://unl.box.com/s/dywe90wwpy0t4ilkuxshkivz2z26mud8

Thanks 4 any help you can provide.

Download coeffs_of_homogen_soln_discrepancy.mw

Dear all,

I'm investigating the vibration performance of timber beams. I have sample data from my test which shows the vibration of the beam. I want to determine the eigenfrequency from this data. The problem I face is that I'm not finding the probber eigenfrequency. I have two data rows; time and amplitude. I'm able to plot the amplitude with SignalPlot but not the time, therefore I have to adjust the samplerate. I have the same problem with the fourier analysis. Is it possible to include the time period as well?

Regards,

Maurits

Dear all,

I need to transforme these equation from time domain to frequency domain with fourier transforms and solve it in frequency domain but i received the flowing error

any helps

thank you !

Code :

Download Fourier_TRAns_MAPLEprime.mwFourier_TRAns_MAPLEprime.mw

Hi, I'm trying to use Maple to construct some examples of symmetry solutions for certain nonlinear PDE's.  As a warm up, however, I'm working through the commands just for the heat equation in 3d: u[t]-u[x,x]-u[y,y]-u[z,z]=0 

I've gotten Maple to produce both determining equations for the symmetry infinitesimal generators via the DeterminingPDE command.  I've also gotten the command Infinitesimals to work too.

However, when I next use PDETools Invariants command, it correctly outputs invariants for most of the generator output of Infinitesimals EXCEPT it won't output anything for the simple rotation generators yd[x]-xd[y].  It will, however, output invariants if the rotation is between an independent and the dependent coordinate.

An example:
with(PDETools)
S:=[_xi[x]=y, _xi[y]=-x, _eta[u]=0]
Invariants(S,u(x,y))

*Above returns nothing, But if you instead have _xi[x]=x and _xi[y]=y then it returns the right invariants.

Thanks in advance!

Hey, this is not the I've had this encounter. I want to open this saved document but when I open it and Maple starts up it just hits me with "A problem was encountered while opening the workbook. Database is not opened". How can I get to open it properly and see my math notes?

How can this be prevented?

When I try to upload the file in this message it says "Cant open a null file"

Any help?

Jacob

plz code and sole that this integral!!!!!!!!!!

!!!

Hi, so I have few problems here. I need to;

Create a MxN Matrix/Lattice, where N and M can be any positive integer, that contains a random selection of -1/1s at each entry.

Need to sum every entry, then multiply by -1 to find “H”.

Need to multiply each neighbour to find its bond energy, so if it’s the same you get 1 else -1, but only its direct neighbours once, so if it was a 2x2 matrix there would be only 4 values and then sum them.

I don’t seem to be able to set up the code so that it does it for any size matrix, as I only know how to write it out basic for a 2x2. Also, not so important, but I wanted to know if could create a loop that would find every iteration of possible setups i.e. for a 2x1 you can have 1.1, -1.1, 1.-1, -1.-1. And then give the solutions outlined earlier for each of the the possibilities [There being 2^MxN ]

Cheers in advance.

In response to Markiyan Hirnyks reference to Mathematica's FindDistribution command http://www.mapleprimes.com/questions/219668-Determine-The-Distribution-Function.  It would be nice to have something similar coded in Maple. 

solve Derivative in maple

see in photo

What do people find a good screen size is? 

Personally I find two monitors would be ideal for using Maple and especially Maplesim.  I'm going to be adding a second monitor, probably a 20".  Maybe 24" is better?  I think 20" would be enough.  Just wondering how many monitors people out there work with and their monitor sizes.

Dear all

I have n, p two integers greater than one

I would like to minimize the following real number by a fixed positive constant if its possible

(n/(n+p))^(n/p)-(n/(n+p))^((n+p)/p)

Many thanks

Dear all,

I am creating an animation, and I was wondering if I can add a multiplier to my equation in a specific range.  So starting from z:=0.2 add a multiplier (z+1). The code I have so far is added. Does anyone know a code for this?

Kind regards

restart; 
with(plots); 
a := -1/2; b := 1/2; c := -2; d := 2; n := 20; 
g := proc (x) options operator, arrow; value(Int(sigma(t), t = 0 .. x)) end proc; 
sigma := proc (z) options operator, arrow; 2*sqrt(2*h^2-4*z^2)*z/h^2 end proc; 
h := i/n; 
for i to n do 
an2[i] := plot(sigma(z), z = -(1/2)*h .. (1/2)*h, view = [a .. b, c .. d], color = AQUAMARINE); 
an3[i] := plot(2*g(x), x = 0 .. (1/2)*h, view = [a .. b, c .. d], color = RED) 
end do; 
p := plots[display]([seq(an2[i], i = 1 .. n)], insequence = true); 
q := plots[display]([seq(an3[i], i = 1 .. n)], insequence = true); display(p, q)

equation 1 : x_i+1=x_i− (f·g_y−f_y·g)/(f_x ·g_y −f_y ·g_x)
equation 2: y_i+1=y_i− (fx·g−f·g_x)/(f_x·g_y-f_y·g_x)

My quesiton are, deriving equations (1) and (2) above and constructing a single Maple function called newt2d that implements both of these recurrence relation.

I apolgize in advance if I don't write my question correctly.  This is my first time posting a question.