## Find limit function for Picard Approximations...

Hello, I need to find a limit function of recurent progression: y0 = a; yn = a + int(x * y(n-1), x) for n=infinity. Thank you very much for any solution P HK

## Drawing straight lines in Sketch...

Is it possible to use Sketch and draw a straight line? I want to use document mode to retype my class notes so I want to include graphs and illustartions. I am taking PHY-2048 where we are studing vectors, dot product, cross product etc.

## Converting Maple to Fortran...

Hi! I have a Maple worksheet comprising of TEXT + MATH, mainly for loops and if statements, using matrices. Can someone please let me know how I can turn this sheet into Fortran code (very detailed instructions necessary!) Thanks for your help! :)

## Primitive Element in GF(2^32)...

I need to generate pimitive polynomials of degree 17 over GF(2^32) for use in an LFSR working over 32-bit words. Does anyone know how it can be done?

## Maple 9 Remembering Function assignments?...

I'm trying to solve a partial differential equation with two boundary conditions below.  The general solution contains arbitrary functions of the non-differentiated variable.  These functions are solved for and assigned but do not appear in the final solution return.  Can anybody help me with this?

> restart;

> l:=lambda;

## coloring lines in contourplot...

Hi, I want to draw the contourplot of this two family of lines x*y and x*y^2. with(plots): contourplot( {x*y, x*y^2}, x=0..10, y=0..20); I would like that the first function (x*y) is colored in red ad the secod (x*y^2) in blue. How could I do it? I tried with "colouring= [red, blue]" but the coloring is not what I want. Thanks in advance for the help. JJ (By the way..I wish you the best to all of you during 2006!!)

## Loops and Conditional Execution Exercise...

This is the simple code that I came up with (after some frustration) to solve the problem of separating the sequence of positive integers from 1 to 20 into one sequence of primes and one sequence of nonprimes, and eventually displaying those sequences. I'm just curious to see what other ways this same question could have been solved. >

## Plot to PS from Maple 10 for Mac?...

Hi. I'm a very new Maple user, but I don't think that this qualifies as a Newbie question. Has anybody had problems when trying to write directly to a PS file from within the Maple GUI on OS 10.4.3 for Mac? Consider the following: currentdir(kernelopts(homedir)): plotsetup(ps,plotoutput="testplot.ps"): plot(sin(x),x=0..2*Pi); Running this code from the command-line works as expected, but running this from the GUI interface usually returns: Error, Error in device driver: plot terminated The strange thing is that about 10% of the time, the above code actually will work from the GUI. I cannot produce failure or success consistently in the GUI. I have tried this on both a G4 and a G5. Am I the only user who has had this problem?

## repartition function in Maple...

I try to deduce distribution fonction from a parametrized set of points with Maple. At this stage I have a dataset define through a relationship like y=f(x), and I want to obtain the distribution of y points, given that x in [-A,A]. The theoretic formula is : F[Y](z) = P(y<z). i.e. F[Y](z) = Int(delta[f(x)< z],x=-A..A); and I use the piecewise Maple function to implement it ( Int(piecewise(f(x)< z,1,0),x=-A..A) ), but for Maple : Int(piecewise(f(x)< z,1,0),x=-A..A) = piecewise(Int(f(x)< z, x=-A..A),1,0)) which is totally different !

## Stand-alone applications...

I was able to create a stand-alone application from a M-file with MatLab 6.0 R12. Is it possible to do the same with Maple 9.5? I mean, can I take a Maple worksheet and convert it into a executable file which can be run on a PC that does not have Maple 9.5 installed? I would like to know since it is easier for me to create GUI applications with Maple 9.5 than with MatLab.

## strange numerical integration...

Using the following integrand

.71428491807021e-1*(1/2-1/2*erf(.50253255206611e-2*xi*2^(1/2)))/Pi^(1/2)*
exp(-1/2*(.14285698361405*xi-5.6995353453900)^2)*2^(1/2)*(1/2+1/2*
erf(1.1171516020037*(5.4091420603773+.44756682898115*xi)*2^(1/2)))

I get -0 for Int( %, xi= -99.496042185589..infinity): evalf(%,14);

which certainly is false. But if I cut off at a reasonable upper bound
(say exp(...) <= 1E-16) I get what I expect (up to rounding errors).

I consider that as a bug and wonder whether it is in the NAG library
or through Maple calling it - any explanation?

## Plotting Modes of Elliptical Wave Equation...

I was wondering if some one had the knowledge as to why the plot3d command will not work for the dirichlet elliptical wave equation, at fixed time t, yet maple is able to evaluate this function at specific values for theta and r. (the coord system for the 3d plot should work with "ellcylindrical" but does not)

the modes of the function look like (for Cosine-Elliptic)

MathieuCE(j,q,theta)*MathieuCE(j,q,I*r)*cos(lambda*t)

j=0, 1, 2...n where q is a 'zero' value of MathieuCE(j,q,I*1) 0<theta<2*Pi, 0<r<1, t>=0

the error stated is "Plotting error, MESH must be a list of lists of vertices, or an hfarray"

## Splitting a polynomial into even and odd terms...

I would like to split a polynomial into even and odd terms. Has this capability been provided in a package? PolynomialTools seems the obvious choice, but doesn't do this. Here's one approach
```SplitPolynomialEvenOdd := proc(poly::polynom(anything,v), v)
description "return the even and odd parts of a polynomial in v";
local p;
p := collect(poly,v);
if p::`+` then
return selectremove(t -> degree(t,v)::even, p);
elif degree(p,v)::even then
return (p,0);
else
return (0,p);
end if;
end proc:

SplitPolynomialEvenOdd(x^2 + 3*x + 1, x);
```