So I am trying to solve a given ODE using calculated christoffel symbols found by maple, and in order to get the correct christoffel symbols, I need my function to be r=x-x_s, where x is not a function of t. However, I then have to solve the ODE where x is a function of t. Maple used the r value to find the christoffel symbols which has x in it, and now I want to find the origonal function of x(t), but I can't have x and x(t) in the same ODE. If I change r=x-x_s to r=x(t)-x_s, I get the wrong christoffel symbols. How can I solve my ODE?

I have a markov matrix that is 500x500. I need to take the standard deviation of the top row of this matrix after I raise it to high powers, like 10^17 by 10^12, (from a do loop.) I would also like ot make a histogram of this data too. 

Any ideas what I can do? 

I have a pair of markov matrices that are 250x250 and 500x500 in size. I am raising these matrices to high powers, like 10^17 by 10^12. (I'm using a do loop.) If I set digits:=15, Maple will use 100% of the processor that is available. If I set digits:=16+ Maple will only use 1 thread of my processor. Why? It is the same matrix same program. The only thing that changes is the digits of accuracy. This seems exactly backwards from what it should be. As a result, my AMD9590 processor will run full blast when I choose 15 digits of accuracy and gets the biggest matrix done in 1-2 hours. If i set the accuracy to 17 digits, it takes half a day or more to do the same thing!!! 2 digits of greater accuracy should not take an additional 20 hours!

Thanks 

Hi

I need a command or code to display numerical amount of parameter (i) in below textbox

e.g. display Element 2 instead of Element i when i is equal to 2

Hello,

How can i from a randomly generated 100 numbers, output the number of unique elements...

thank you very much.

hello, 

i need to use only looping to determine the larget integer in my random list...

here is how i put it together, but my result is incorrrect... Please help.

thank you

R:=rand(1..50)

seq(R(1..50),i=1..20)

L:=[%]

max:=proc(L,maxv : : evaln)

locali;

maxv:=L[i];

for i from 1 to nops(L) do

if eval(maxv)<L[i] then maxv:=L[i]end if;
end do;

end proc;

Hi all,

I want to design a procedure which can generate multiple loop automatically. While I  assign a value to it, it will create a multiple loop with same number of layers.    

For example,

P(1):

for i from 1 to n do
...
end do:

p(2):

for i from 1 to n do

for j from 1 to n do
...

end do: end do:

p(3)...

p(4)...

...

I'd appreciate any help on this topic. Thank a lot.

hello,

i am quite new into programming and i need help with writing a code that determines the largest integer in the list...

Please help..

Thank you 

I try to solve a set of differential equations with boundary conditions by dsolve command but I got this error:

--------------------------------------------------------------------------------------------------------------------------

odes := D2*(diff(P(x), x))/((3*D1*a+4*D2)*P(x))-(diff(S(x), x))/(q*S(x)-1) = 0, diff(S(x), `$`(x, 2))+(diff(S(x), x))*cotx+4*pi*(3*D1*a+4*D1)*P(x)/((q*S(x)-1)*D2) = 0

ics := P((1/2)*pi) = 1, S((1/2)*pi) = -1, (D(S))((1/2)*pi) = 0, (D(P))((1/2)*pi) = 0

sol := dsolve({ics, odes}, numeric);
Error, (in dsolve/numeric/process_input) input system must be an ODE system, got independent variables {x, (1/2)*pi}

--------------------------------------------------------------------------------------------------------------------------

Does any body knows what the problem is??

Hello

I wantt to assigne value to a variable that doesn't exist before the radio button. How can i do that?

Could you please help me?

Hey guys. I've tried to find an answer for this, but have struggled since our learning book is in danish, so the used terms may not be technically correct, so sorry for poor phrasing.

Anyway, how would you solve this problem in maple?

Find a solution for the differential equation:

d⁴y/dt⁴ - 16y = u' + u

With the effect*  u(t) = e^3t + 3e^t

I've gotten this result (By hand calculation)

y(t) = 4/65*e^3t - 6/15*e^t

Thanks for the help. It's my first post, so let me know if I should do something different next time :)

-Alex

Bonus question:

How do you calculate the transfer function in maple:

H(s) = (s + 1)/(s⁴-16)

*Don't know if 'effect' is the correct term. 

See subject. Via Google I only found links for converting from decimal to binary. But how do I write a binary number in the first place? I don't have a decimal to convert from, I have long binaries to start with! Which I don't want to convert to decimal, since I want binaries as output too.

Not sure exactly how i could achieve this but:

how do i determine the value of k for which the graphs p(x) = x^2+2x+3 and q(x) = k+5x-7x^2 enclose an area of exactly 36?

I have to do it in maple and using i guess area under the curve.

Thanks

I need to show that the least square solution x that I obtained as x=pseudoinverse(A).b is the solution of Transpose(A).A.x=Transpose(A).b with the smallest norm. I've obtained the norm for the RHS of this expression as well as the norm of x but I'm unsure of how to conclude that this is the lowest possible norm using these values

The issue I am currently having is that, while analyticity (and physics) indicates a certain function must have roots, fsolve is having trouble finding them. In fact, I have even found roots manually in a certain region myself, simply inputting into the function various values until I found them. However, fsolve does not seem to want to find these roots, and I believe it is a numerics issue: when I changed the digits around, for extremely low values of Digits, it would find a root (even though it was incorrect). Further, this problem arose elsewhere in the domain of interest for other values of Digits (in particular, for Digits:=5, fsolve failed in a region it had not failed before).

The region of interest is the "peak" of the output of poleR(M0, 0.935, mK), which should be somewhere around M0 = 0.95 or so. However, because fsolve cannot find the roots, the plot cannot be made.

Anyone have any ideas as to why fsolve cannot find the roots? I was also experiencing issues with some of these functions having multiple roots, which itself is weird as well (note that I am working over the complex plane).

Attached is the document.

pole-dragging-mapleprimes.mw

Any ideas?