How to plot the function f := x->sum(sin(x/(k+1))/k, k = 1 .. infinity) on the interval 0..50? The series obviously converges.

Hello,

I have 2 questions about the resolution of a ODE for the simple pendulum equations.

1) I need to add 'maxfun'=10^7 in the dsolve function. With this option, dsolve works. Can you explain me in a very simple way the principle of this option ? It's just to feel better how this option runs.

2) gamma(t) is a angular parameter.

With

tgammax:=odeplot(sol,[gamma(t),x(t)],0..10,numpoints=200,color=blue,legend="z"):
display(tgammax);

I manage to plot x(gamma) but with gamma in rad.

How can i do to x(gamma) but with gamma in degree? I should do gamma_deg=360/(2*Pi)*gamma(t) but as gamma(t) is defined numerically, this lign doesn't work.

Thank you for helpe

Here some ligns of my code :

Cinit:=D(psi)(0)=0,psi(0)=-0.93,D(x)(0)=0,x(0)=1.13,D(gamma)(0)=0,gamma(0)=0.79;sol:=dsolve({sys,Cinit},numeric,'maxfun'=10^7);
tx:=odeplot(sol,[t,x(t)],0..10,numpoints=200,color=blue,legend="x"):
display(tx); 

tgamma:=odeplot(sol,[t,gamma(t)],0..10,numpoints=200,color=blue,legend="gamma"):
display(tgamma);

tgammax:=odeplot(sol,[gamma(t),x(t)],0..10,numpoints=200,color=blue,legend="z"):
display(tgammax);

Hi, I am trying to plot these two curves:

I tried:

with(plots);  A := Array(1..2)  A[1] :=plot (0.199563349672261+0.0178636902277546 x^1.14406289706794-0.0182070811144750 x^(1.13867380551454),x=50..2050,  color=red);  A[2]  := plot(0.298910542599302+0.0117459591500434 x^1.00390277106937-0.0137065176395662 x^0.970667551759677, x = 50..2050, color = blue);  display (A);

and I got:

Error, (in plot) unexpected option: .298910542599302+0.117459591500434e-1*x^1.00390277106937-0.137065176395662e-1*x^.970667551759677

so I tried

so I am not sure how I would graph these two functions...

I am coming back to maple after being years out of it (Maple 9). I am writing a paper in class of how to use Maple with Algebra. I would like to create a document so that students can manipulate the graph of  and se how parameters a, h and k effect the function.  However, I can't seem to figure out how to create a slider linked with each parameter so students can manipulate it. 

Is this possible? Any help wold be great!

Nicholas

I've been having trouble plotting in polar coordinates. My code is included below:

The ouput prints the 10 terms of the sum but the plot does not produce anything.

I also would like to know a method that I could use to plot with a going from zero to the size of the radius.

Thank you for any help.

Is it possible to show using maple if a given line in 3d space intersects a 2d plane with given co-ordinal vertices?

Ive been given a task to write a fuction that will print true/false to if a given line will intersect a triangle given the co-ordinates of the vertices. 
I have a general idea on what needs to be done mathmaticly, but i have poor maple skills and the syntax is killing me.
I would appreciate some help in this regard.

Thanks.

I've been on this question a week now and still no conclusive answer! What I need is a function that produces the inequalities that determine a triangle given the 3 points and then using a 4th point, prints true if the 4th point satisfies 2 or 3 or the inequalities and prints false if it only satisfies 1 or none of the inequalities. I need to have this solved by tonight so any quick help would be greatly appreciated!

Consider the curve defined by f(x, y) = 3+2x+y+2x^2+2xy+3y^2 = 0.Locally on the curve we can view y as a function of x, i.e. y = y(x).Compute formulas for the first and second derivative of y with resoect to x.

2. Give the Maple command(s) to compute \frac{\partial^8 f}{\partial^5 x \partial^3 y} for f(x, y) = e^{2x+ cos(y)}.

How to show the code in file with extension .lib and edit it?

I am stuck with an IVP which is

eq1:=diff(y(x),x$2)+2/x*(diff(y(x),x))+y^M=0;

ic:=y(0)=a,D(y)(0)=0;

its quite easy to find the series solution of the ode 

dsolve({eq1, ic}, y(x), series);

y(x)=a-(1/6)*exp(M*ln(a))*x^2+(1/120)*(exp(M*ln(a)))^2*M*x^4/a;

But I am facing problem when I try to solve it numerically,

dsolve(subs(a=1,M=3,{eq1,ic}),numeric);

THanks

Hello

Am trying to convert a vector projection to matrix type output.

Attached is a sample worksheet.

Transformation_matri.mw 

Consider:

> restart:
> int(exp(I*x)/x,x=-infinity..infinity,CauchyPrincipalValue)
                                                                                                   Pi I

> int(exp(I*3*x)/x,x=-infinity..infinity,CauchyPrincipalValue)
                                                                                                   Pi I

> int(exp(I*k*x)/x,x=-infinity..infinity,CauchyPrincipalValue) assuming k::positive;
                                                                                                  2 I Pi

I would expect the last line to also give me Pi I

Mathematica gives a Pi I via:

FullSimplify[Integrate[Exp[I k x]/x, {x, -\[Infinity], \[Infinity]}, PrincipalValue->True, Assumptions->{k>0}],  Assumptions->{k>0}]

Any idea what this could be?

thanks.

Congratulations to Eric Miles, a graduate teaching assistant at Colorado State University, for winning the second individual prize of The Möbius App Challenge! Eric will be receiving an iPad prize pack for his app Stretch and Translate Sin(x).

Maplesoft launched the Möbius App Challenge earlier this year to give Maplesoft community members the opportunity to get involved in the creation of Möbius Apps. This contest seeks to challenge users to create Math Apps using Maple, submit them to The Möbius Project website, then award the most useful, creative, and effective submissions. Finalists from the individual level challenge are chosen four times a year, so make sure to enter the next round of the challenge for your chance at a Maplesoft prize pack!

For full contest details and to enter, please visit: http://www.maplesoft.com/lp/mobius/index.aspx.

Jonny Zivku
Maplesoft Product Manager, Maple T.A.