MaplePrimes Questions

I'm begginer and I don't know, whether this equation can be solved in Maple:

equation

The only variable is N.

i=1,2,3,...,n

n=the nearest smaller natural number from the second expression.

The equation sides probably won't be equal precisely, also I'm looking for some N when equation side's values are close enough.

Can you help me with programming this into Maple?

 

restart

int(sin(x), x)

-cos(x)

(1)

int(x^2, x)

x^2*_X

(2)

``

 

Download 12.mw

Given the sequence defined by the recursive relation a[n+1] = r*a[n](1-a[n])
You need to use the procedure iterate.
Throughout this problem you should choose initial values in the interval 0<a0<1.
(a) Let r=3/2. Calculate a moderate number of terms in the sequence (between 10 and 20). Does the sequence appear to be converging? If so to what value? Does the limit depend upon your choice of initial value? Plot the terms you have calculated
(b) Let r=2.8. Calculate a moderate number of terms in the sequence (between 10 and 20). Does the sequence appear to be converging? If so to what value? Does the limit
depend upon your choice of initial value? Plot the terms you have calculated How does this sequence differ from that in part (a).
(c) Let r=3.2. Calculate a moderate number of terms in the sequence (between 10 and 20). Show that the sequence does not appear to converging. Plot the terms you have calculated and describe how the sequence behaves in this case.
(d) Consider intermediate values between 2.8 and 3.2 to determine more precisely where the transition in behaviour takes place. Provide a few plots (no more than 4) showing the values you have investigated.
(e) Consider the values of r in the range 3.43<r<3.46. Determine as accurately as you can the value of r for which the period of oscillation doubles.
(f) As r increase further period doubling occurs. Try to find the when the sequence appears to oscillate between 8 values.
(g) Let r =3.65 and calculate a considerable number of terms (at least a few hundred) and plot your values.
(h) For r=3.65 choose a0=0.3 and then a0=0.301. Find and plot some terms in the sequence for each initial value. Determine how long the terms in the two sequences remain close together and when they begin to depart significantly from each other.

Hello,

In this procedure I get the error, 

Error, `:=` unexpected

I know what it means I just can't seem to resolve it. Any help would be greatly appreciated! Thank you in advance for looking at this code!

Kind regards,

Gambia Man

HamilMat := proc (K::integer) local ni, mi, nj, mj, N, Hamil, Eigenvec, i, j, res; option remember; global Vij, U, L; N := K^2; ni := Vector(N); mi := Vector(N); nj := Vector[row](N); mj := Vector[row](N); for i to N do for j to K do res := (i+K-j)/K; if type(res, integer) = true then ni[i] := j; nj[i] := j; mi[i] := res; mj[i] := res end if end do end do; Hamil := Matrix(N, shape = symmetric); for i to N do for j from i to N do if i <> j then Hamil(i, j) := Vij(ni[i], mi[i], nj[j], mj[j]) elif i = j then Hamil(i, j) := Vij(ni[i], mi[i], nj[j], mj[j])+(1/2)*(ni[i]^2+mi[i]^2)*Pi^2/L^2 end if end do end do; return Eigenvec := Eigenvectors(Hamil, output = ['values', 'vectors']), Hamil end proc

Error, `:=` unexpected

 

 

Download small_error.mw

There are no users with 0 reputation.  It appears all users with 0 reputation and negative reputation have been erased.  One user I can not find who is or now was a legitimate user is John Mcloone an employee at Mathematica who made a post here.  I can only think during the recent spam attack that all users with 0 or negative reputation were removed.  Some of those users had legitimate questions.  Where did those users, John Mcloone and their posts go? 

Hi all,

 

I tried to find the real solution of the unlinear integral equation:

 

exp(-h^2/T)*(Int(exp(-x^2/T)*BesselI(0, h*x/T)*x, x = 0 .. 1))/T

 

but I got the warning and an complex solution:

 

 solve(subs(T = 1, eq)-.99 = 0, h)

Warning, solutions may have been lost

-1.232350242*10^(-32)-1.130417828*I

 

Can anyone help me to find a real solution for this issue (if possible)...?

I would like to thank you in advance.

 

Hi. I have a commutation table, that I've got as output of Greg Reid algorithm.

Can I use to define Lie algebra for components DifferentialGeometry & Lie Algebras?

Worksheet is attached...

FCT.mw

 

In the following problem though b and c are same (except the way denominator 2 is hanfled), command ' a-b ' readily answers zero, but a-c not so. Why? Only on condition of assumption real it gives zero!

a := (1/2)*(kappa*omega^2+omega^3)*(Y+(1/2)*(-sqrt(N)*omega^(3/2)*sin(theta[2])*cos(varphi[2])*lambda__b+sqrt(N)*omega^(3/2)*sin(theta[1])*cos(varphi[1])*lambda__a)/(kappa*omega^2+omega^3))^2/omega:

b := (1/2)*(kappa*omega^2+omega^3)*(Y+(1/2)*(-sqrt(N)*omega^(3/2)*sin(theta[2])*cos(varphi[2])*lambda__b+sqrt(N)*omega^(3/2)*sin(theta[1])*cos(varphi[1])*lambda__a)/(kappa*omega^2+omega^3))^2/omega:

a-b;

0

(1)

c := (1/2)*(kappa*omega^2+omega^3)*(Y+(-sqrt(N)*omega^(3/2)*sin(theta[2])*cos(varphi[2])*lambda__b+sqrt(N)*omega^(3/2)*sin(theta[1])*cos(varphi[1])*lambda__a)/(2*(kappa*omega^2+omega^3)))^2/omega:

a-c;

(1/2)*(kappa*omega^2+omega^3)*(Y+(1/2)*(-N^(1/2)*omega^(3/2)*sin(theta[2])*cos(varphi[2])*lambda__b+N^(1/2)*omega^(3/2)*sin(theta[1])*cos(varphi[1])*lambda__a)/(kappa*omega^2+omega^3))^2/omega-(1/2)*(kappa*omega^2+omega^3)*(Y+(-N^(1/2)*omega^(3/2)*sin(theta[2])*cos(varphi[2])*lambda__b+N^(1/2)*omega^(3/2)*sin(theta[1])*cos(varphi[1])*lambda__a)/(2*kappa*omega^2+2*omega^3))^2/omega

(2)

"(->)"

0

(3)

``

Why the answer is not given as zero?

``

``

 

Download what_is_the_difference_between_b_and_c.mw

What difference therms b and c make for Maple? Are they not same?

Ramakrishnan V

rukmini_ramki@hotmail.com

On the Physics Research and development page there are only 3 versions of the Physics package available.  What happened to the earlier final updates for Maple 16, 15 etc...

What is the latest one that can work on Maple12?  Was there one for M12?  What is the earliest versions of Maple the Research & Development Physics packages can work on?

are these quations in wave form for harmonic ocillator?

i can not plot them

plot(exp(I*2*t) - exp(-I*2*t), t=0..10);

plot(cos(t) - I*sin(t), t=0..10);

 

why same solution but their plot graph are different?

plot(cos(t) + sin(t), t=0..10);
complexplot(exp(I*t), t=0..10);

 

Hello all, first time posting here.  Anyways, I am trying to make a velocity vector plot from known u and v velocity components at known x and y locations i.e. at (0.1, 0.1), u=-0.5 m/s and v=0.43 m/s.  My calculations have been performed with a FORTRAN code, and the results output in matrix format.  What I am left with are 4 vectors: x location, y location, u velocity component, and v velocity component.  I know that MATLAB has a specific function, called "quiver," for creating a velocity vector plot from known data points (through the use of vectors containing those points).  I have trouble believing that there is not an equivalent function, however convoluted it may be, in Maple to create these plots.  I have tried to manipulate the fieldplot, vectorplot, and arrow functions to no avail.

Any ideas on how to create a velocity vector plot from known data points?  I am using Maple 18.

Here is a link to an example of what I am looking for:

https://people.sc.fsu.edu/~jburkardt/py_src/spiral_data/spiral_matplotlib.png

Thanks, Brandon.

I have used data and the curve fitting procedure to find a curve and I now have to put error bars on the graph to complete my task? I know that there is an errorplot function but I cant seem to find a way to do it with the curvefitting thing?

 

I have attached the whole worksheet, I hope someone can help! Thanks Detectors_Tutorial.mwDetectors_Tutorial.mw

 

Even if there was some method of exporting the graph to another medium that would also be very helpful

 

I use Maple18 on unix 64. I have the MAPLE and the LD_LIBRARY_PATH properly set up to use OpenMaple. When I link my code to OpenMaple library I get segmentation fault when I try to debug using gdb. In order to isolate the problem, I created a simple temp.cc file that contains the following:

 

#include <string>

#include <stdio.h>

 

int main()

{

int x = 10;

std::string abc;

abc = "testString";

printf("this is it %s\n", abc.c_str());

}

 

When I compile the code using g++ -g -o temp temp.cc I can run ./temp and step into the code using gdb with no problems.

 

If I link to Maple Library g++ -g -o temp temp.cc -L/opt/maple18/bin.X86_64_LINUX –lmaplec

 

I can run the code ./temp but I can not step into it using gdb. I get the error as shown below when I try to step over std::string abc.

 

//-----------------------------------------------------

 

GNU gdb (GDB) 7.8

Copyright (C) 2014 Free Software Foundation, Inc.

License GPLv3+: GNU GPL version 3 or later <http://gnu.org/licenses/gpl.html>

This is free software: you are free to change and redistribute it.

There is NO WARRANTY, to the extent permitted by law.  Type "show copying"

and "show warranty" for details.

This GDB was configured as "x86_64-unknown-linux-gnu".

Type "show configuration" for configuration details.

For bug reporting instructions, please see:

<http://www.gnu.org/software/gdb/bugs/>.

Find the GDB manual and other documentation resources online at:

<http://www.gnu.org/software/gdb/documentation/>.

For help, type "help".

Type "apropos word" to search for commands related to "word"...

Reading symbols from temp...done.

(gdb) b main

Breakpoint 1 at 0x40091f: file temp.cc, line 6.

(gdb) r

Starting program: xxxxxxxxxxxxxxxxxxxxxxxxxxx

[Thread debugging using libthread_db enabled]

Using host libthread_db library "/lib/libthread_db.so.1".

 

Breakpoint 1, main () at temp.cc:6

6         int x = 10;

(gdb) n

7       std::string abc;

(gdb) n

BFD: reopening /tmp/xf-dll/xf-1274259bc6684076beea8652f625b03bf2133.tmp: No such file or directory

 

BFD: reopening /tmp/xf-dll/xf-1274259bc6684076beea8652f625b03bf2133.tmp: No such file or directory

 

BFD: reopening /tmp/xf-dll/xf-1274259bc6684076beea8652f625b03bf2133.tmp: No such file or directory

 

Can't read data for section '.eh_frame' in file '/tmp/xf-dll/xf-1274259bc6684076beea8652f625b03bf2133.tmp'

0x00000000004007a8 in std::basic_string<char, std::char_traits<char>, std::allocator<char> >::basic_string()@plt

    ()

(gdb) n

Single stepping until exit from function _ZNSsC1Ev@plt,

which has no line number information.

Segmentation fault

 with small step rather than two points?

I had writen my questions into this code hw2_finished.mw And I think there's maybe something wrong with A=[0..130]

I think you can understand what I want by reading this code but if you are confused please let me know.

Dear all

I have query related to generating commutator table from given Lie algebra spanning symmetries of PDE, I have ten vectors that I need to commute with each other to constitute commutator table. I have calculated all the commutations manually using "symmetrycommutator" command.

My queries is, can I convert this operation into short loop type program which automatically calculate cummutation of vectors from 1 to 10.

The Maple file is attaches with this query for reference.

Commutator_for_Infinitesimal_Generators_for_variable_coefficient_Zakharove_Kuznetsov_Equation.mw

First 1199 1200 1201 1202 1203 1204 1205 Last Page 1201 of 2428