Greeting

can u help me how to solve this problem using the galerkin finite element method with graphs.  I attachment below.

IBC

Hi everyone, I've just learnt some basic tools and I need your help. 

Hello all,

I would like to know how to make computations on Quantum Field Theory questions.

For example: to write Feynman diagrams for some interactions, obtain the S matrix of scattering amplitude, then deduce the diffrential cross section.

If anyone is working on these issues in the Maple community, I would be very interested.

Best regards

Jean-Luc Paillet

Hello all

when i tried to use the library "with(VLA)" i always have this problem :

 with(VLA);
Error, invalid input: with expects its 1st argument, pname, to be of type {`module`, package}, but received LA
 

but i have sow in many videos that its working and they applied it in same way.

any help please > >  

Hello people in mapleprimes

I asked a question for a similar thing before.
https://www.mapleprimes.com/questions/222751-PDEtoolsdeclare

The question I had then is not solved yet in my mind, so I think  I will ask you
again. It is about differentiation of a composite function.

with(PDEtools):
declare(f(x,y,z),g(t,y));
a:=f(x,y,z);
b:=x=g(t,y);
c:=eval(a,b);
d:=diff(c,y);
OFF;
e:=d;

h:=convert(e,diff);

For both e and d, the notation of D__1 and D__2 appear.
And, even if, with convert, I changed them to another form,
what I can get is more complicated form of h.

Is it inevitable to use D__1 and D__2 brabra in the differentiation of
composite function?

I am writing this question thinking that if there is a way of go around it,
 I want to know.

And, though this is the second question, even if I use latex( ) command,
natually, codes of Tex I can get is that of the output shown on the screen, including D__1 and D__2, even though I want it to be shown with f__x and f__y or that including prime. Then, I have to edit that tex file, changing D__1 to f__x etc. which I think is very complicated modification of the text file.
I wish there is a way to modify it in Maple so as it to be direct TeX codes which  do not require me to modify.
Is it impossible to do such a thing?

Thanks in advance.

diff_2.mw

How to find the range of poset?

i would like to search the range of range is equal to the range of domain after solve system.

it may not exactly equal but around 80% of elements in range belong to the range of domain

if it is not factorisation domain,

how to generate ideal with function in maple?

Is there a way in Maple to numerically integrate (with some symbolic variables retained), involving vectors?

Alternatively, is there a way to simply numerically integrate the following expression involving vectors, without any remaining symbolic variables remaining?

Integrate over vectors p1, p3
(with vector x remaining as free variable after integration)
( e^{-i p1 . x} (p1² + p3²) )  / ( (p3 - p1)² (p3²*a² +1)² )

If absolutely required for numerical integration, i.e. no other way to get maple to perform a semi-numerical integration, then vector x, and the scalar variable "a" can also be specified a value, but x should remain a vector. 

Of course, if the integration above can be done analytically, or even partly analytically (e.g. if the vectors are expressed in spherical polar co-ordinates, and some of the variables like Sin \theta etc. can be integrated over), that would be very useful as well. 

Hi,

I am trying to complete this assignment but I am having trouble on this question. It is asking me to use a graphical approach to estimate a 𝛿 > 0 such that for all x...

0 < | x - c | < 𝛿 implies that | f(x) - L | < ε

It says I can do this by plotting my f(x), y₁, and y₂ over the interval [c - 𝛿, c + 𝛿] with a y-range of [L - ε, L + ε]. y₁, and y₂ are banding lines defined by ( L - ε ) and ( L + ε ), respectively. The value of ε is shown below. I dont know where to start but this is what I have so far...

f(x) := ( x ( 1 - cos(x) ) ) / ( x - sin(x) );

The limit ( L ) of this function as x approaches 0 ( c = 0 ) is...

L := ( limit ( f(x), x = 0) );

The limit is 3 ( L = 3 ).

 ε = r = Student Id / 5000000   = 0621748 / 5000000   = 155337 / 1250000

 ε := 155337 / 1250000;

y₁ := ( L - ε );

y₂ := ( L + ε );

Originally, it had asked me to graph the function with the banding lines (y₁, y₂) together within the interval [-0.2, 0.2], which I have done. I just do not know how to find a value for 𝛿 in this case. Please help.

Thanks.

This doesn't work directly:

is(RootOf(_Z^3-_Z-1, index = 1)+RootOf(_Z^3-_Z-1, index = 2) = 0);
                                      FAIL

But Maple can compute the bounding interval for a polynomial RootOf (convert(..., RootOf, form=interval)), and evalrC can perform operations on complex intervals:

evalrC(RootOf(_Z^3-_Z-1, index = 1)+RootOf(_Z^3-_Z-1, index = 2));
       INTERVAL(0.6623589782, 0.5622795119, 0.6623589795, 0.5622795122)

evalrC(Re(%));
Error, (in evalrC/INTERVAL) not implemented yet: 4

Re and Im are obviously positive though. Why doesn't Maple try to do that automatically, and also, how do I check if the point is inside the interval? (Without writing my own function to handle real intervals, unions of real intervals, and complex intervals.)

the infint  summation  maple cant evaluated it  when my variables  are more than one decimal number 

;

In the link below is my code which generates numerical data of a Fourier series of an impulse function (not to be confused with the DIRAC DELTA IMPULSE).  The issue is as k increases the harmonic frequencies become faster than the time step which causes resolution issues at higher k values.  Is there a way to iteratively change the time step for each iteration of k value?  However, the range for t would have to correspondingly decrease or the file size would balloon by the corresponding change in k value.

So lets say my time step for k = 1 is 0.001 and my time range is 0..1, now for k = 1000 my time step would be 0.000001, but my time range would shrink from 0..1 to 0..0.001.  I hope this is not too confusing.  Maybe experimenting with my code will aid in understanding my request.

Thanks for any assistance

MATRIX_loop.mw

Maple will multiply two matrices for me perfectly.  

However, when I change the value in one of the cells and try to re-evaluate the product, Maple starts giving me the sum of something in all of my cells.   

Are you simply not allowed to ever change the initial values you set for a matrix?  This doesn't seem to make sense to me.  

I have attached a screen shot to show an example with a random matrix.  

1.

eliminate({-a+x+y+z, -b+x*y+y*z+x*z, -c+x*y*z, -P+x^4+y^4+z^4}, {x, y, z});
    [{x = (sqrt(c^2*(a^2*b^2-2*a*b*c-4*b^3+c^2)/b^4)*b^2+a*b*c-c^2)/(2*b*c),
      y = -(sqrt(c^2*(a^2*b^2-2*a*b*c-4*b^3+c^2)/b^4)*b^2-a*b*c+c^2)/(2*b*c),
      z = c/b}, {-a^4+4*a^2*b-4*a*c-2*b^2+P}],
    [{x = -(sqrt(c^2*(a^2*b^2-2*a*b*c-4*b^3+c^2)/b^4)*b^2-a*b*c+c^2)/(2*b*c),
      y = (sqrt(c^2*(a^2*b^2-2*a*b*c-4*b^3+c^2)/b^4)*b^2+a*b*c-c^2)/(2*b*c),
      z = c/b}, {-a^4+4*a^2*b-4*a*c-2*b^2+P}]

The expression for P is fine, but how are x and y square roots and z=c/b??

2.

eliminate({y*z, x+y-z}, z);
                       [{z = 0}, {x + y}]

Is this supposed to happen? I was expecting y*(x+y):

Groebner:-Basis({y*z, x+y-z}, lexdeg([z], [x, y]));
                     [x*y + y^2, z - x - y]

Although a very general question... does the rope break or not?

A 100 kg weight is hung from the center of a 1/4" polypropelene rope making an angle of 10 degrees at both ends with the horizontal.  According to some characteristics, the breaking strength of the rope is 1125 lbs or 5kN.  The safe load (safety factor of 12) drops the weight to  93.8 lbs or 0.417kN. Find the tension in the rope?  Is it strong enough?

Fairly easy to solve in Maple - just two equations two unknowns making t1 tension on one side and t2 tension on the other side

eq1 := t1*cos(80*Pi*(1/180))+t2*cos(80*Pi*(1/180)) = 100*9.8 # ΣFy=mg

eq2 := t1*sin(80*Pi*(1/180))-t2*sin(80*Pi*(1/180)) = 0 #ΣFx=0 hence basically t1=t2

solve([eq1,eq2])

                 {t1 = 2821.797537, t2 = 2821.797537}

Therefore the tension on the rope t1=t2 is 2821 N  (looks like way beyond the safety factor) so this means the rope will not break, correct?  or is it that the total tension becomes t1+t2=5642 N putting us over the 5kN and the rope breaks?

I am new to maple. how do i segregate different term in long PDE in maple