Is there a way to have Maple compute the limit of a recursive sequence? For example, the sequence below

a_n = (a_(n-1) +2)/(a_(n-1) + 3), a_1 =2

is bounded and montonically decreasing, thus has a limit, but Maple is not able to compute the limit. Any ideas?

Hi,

This is for a code that analyzes composite laminates with many plys of varying orientation stacked on top of each other, and this is the first semester I've used Maple.

I'm solving for the roots of an equation which will always have a positive and a negative root, and I do this for each ply (could be 100 plies in a laminate).  I want to collect the positive root for each ply, put them in a list or array, and then pick out the smallest value of all the plies ...

I'm trying to figure out the forces on this arm crane when alpha angle is changed from 0 degrees, 30 degrees, 45 degrees, 60 degrees and 90 degrees. What I'm confused about is how to input the angle in maple. Do I add 45 degrees and alpha angle together?

Apparently the load (L) is not unknown but I'm not sure how its not unknown when I am not given any number for it. How do you find out the load (L)?!

Here's the picture:
http://postimage.org/image/x4s13q3b7/

Dear friends,

I am reporting with a brief comment concerning the integral int(1/(1+x^a), x=0..infinity) with a>=2 a real number. This was evaluated here.

Now Maple 15 (X86 64 LINUX) will quite happily compute this in its most simple form involving the sine when a is not a positive integer or a rational number. If it is, however, a beta function term results,...

Hi,

I am new to Maple and I am trying to write a CFD code. I want my code to be as modular as possible, therefore, I want to use subroutines  in my code. (Because I will need to do some steps several times )

Here is one of the main subroutines that I want to call for it in various times with different parameters. I am familiar with matlab so let be show what I mean first in matlab notation.   Lets say my subroutine is called MAT 

I am very new to Maple. I learned how to graph polar coordinates. The next thing I want to learn how is to highlight the area (not to calculate, but to show) between two polar graphs. Can you tell me what is the code for it? This is my set up:

with(plots);

polarplot([2+cos(x), 2], x = 0 .. 2*Pi, color = [red, blue]);

Thanks.

Hi I am new to using Maple,

Mybe this is a very basic question.

I am trying to find solutions to the following nonlinear system of equations using fsolve.

This is the right hand side of a system of ODE.

nonlinearsystem.mwnonlinearsystem.mw

I do get two solutions for the given parameter values. My question are:

As discussed in http://www.mapleprimes.com/questions/144586-Op-Command-Isnt-Working-As-Expected , the indices change some time at run time so previously working command is not working for me.

I am trying to use 

eq3 := applyop( u->Parts( u, phi[i](y) ), { 1,10,13,20 }, eq2 );

which clearly use the 'op' to get operands from eq2, but since the indices change at every run or on some run, the code is useless. The problem is...

Hi all,

I am new to Maple and I have (possibly a simple) question. I'm solving 2 coupled differential equation, it all goes well, until I want to plot them.

I have the output in the following form:

sol = {w1(x) = [ -a*x^4+b*x^3-c*x]}, {w2(x) = [-d*x^4+e*x^3-f*x^2+gx-h]}

My goal is now to plot the solution in 1 figure with w1(x) between x=0..100000 and w2(x) between x=100000..150000

I have searched on the web for a way to get the formula's of w1(x...

Hello

I am trying to plot solution of ode0 together with the maximum and minimum values but I am having difficulty since the first plot is a solution and second is values. I should have a plot with two line one represent the solution of ode0 and second (the max an min). Any advise or suggestion?

This is the code:

> restart;
with(DEtools); with(plots); Nsols := 5; Ntstep := 10;
 k := 0; A := 0.37e-1; B := 0.2e-6;
ode0 := diff(U(t), t) = -(A+B*U(t))*U(t);

I don't understand why 'expand(x+y)^18' the smart popup work, but with 'expand(x+y)^19'  (or more then 19) i have no smart popups!

Any idea?

Thank you all

Hi 

I am a Matlab user and just started using maple because I need more precision for my calculations. I failed to make Matlab interface work, therefore, I am writing my codes from the begining. And I am having a problem with a nested for loop.

The thing that I am trying to make is (in Matlab)

for j=1:N;
    for i=mod(j-1,2):2:j;
D(i+1,j+1)=2*j;
    end
end

In maple I tried something like

Hi Maple Community!!

Using dsolve, Mathieu 15 gives me the following solutions:

Given that I have specific values for b, e and g. How can I plot x(t...

I have a matrix of order 14 with whose entries being variables about u_1 to u_12. I want to get its determinant, but it return no results. It explains that it is too big for maple to deal with. So I wonder how to deal with such kind of this problem?

Thank you very much for kind attention!

I am trying to solve a system of 4 odes by using the following program:

with(plots):

Pr:=6.2:Ha:=10^(-12):Nt:=0.1:Nb:=0.1:Nr:=0.1:Gr:=0.1:Ec:=0.1:Rb:=0.1:Lb:=1:Pe:=1:tw:=0.2:Le:=5:

fcns:={f(y),T(y),S(y),G(y)}:

sys:=diff(f(y),y$3)+f(y)*diff(f(y),y$2)-diff(f(y),y)*diff(f(y),y)-Ha*diff(f(y),y)+Gr*(T(y)-Nr*S(y)-Rb*G(y))=0,

diff(T(y),y$2)+Pr*f(y)*diff(T(y),y)+Nb*diff(S(y),y)*diff(T(y),y)+Nt*diff(T(y),y$2)*diff(T(y),y$2)+Pr*Ec*(diff(f(y...