Hello,

please explain how to write a code to calculate and output the actual area using integration for y=X^3 over range (0,2) using left-hand rule and 200 subdivisions?

Thank you 

I currently have a function quadsum(n) that determines the [x,y] solutions of the above equation for an integer n. :

quadsum:= proc(n::nonnegint)
local
k:= 0, mylist:= table(),
x:= isqrt(iquo(n,2)), y:= x, x2:= x^2, y2:= y^2;
if 2*x2 <> n then x:= x+1; x2:= x2+2*x-1; y:= x; y2:= x2; end if;
while x2 <= n do
y:= isqrt(n-x2); y2:= y^2;
if x2+y2 = n then k:= k+1; mylist[k]:= [x,y] end if;
x:= x+1; x2:= x2+2*x-1;
end do;
convert(mylist, list)
end proc:

How would I alter this so that I get [x,y] for n= (5^a).(13^b).(17^c)(29^d) for non-negative integers a,b,c,d?

Hi everyone,

Consider the following worksheet (Maple 13):

Euler.mw

 (or as a Google Drive link)

https://drive.google.com/file/d/0Bzr3EyK8arkOUkZPcEE5b0NySVU/view?usp=sharing

Choose execute entire worksheet on both Maple 13 and Maple 18.

Maple 13 is blindingly fast at 40 seconds, entire wroksheet upon completion

Maple 18 crawls and hangs forever.

Can anyone recognise what the problem with Maple 18 is?

Running on a two core 1.7GHz AMD Athlon.

Many thanks

--

Yiannis

Hi

I can not my software,I get this error"the specified procedure could not be found"

what should I do?

Given a  vector v (parallel to the principle axis) and a center point c,

and generating vectors u(beta) orthogonal to v for 0 < beta < 2pi,

and real positive parameters a and b 

I want to plot the hyperboloid  c + a*v*sec(alpha) + b*u(beta)*tan(alpha)

for -pi/2 < alpha < pi/2.

I am able to generate something of a plot, but can't control the size and appearance.  

  Hi there,

  I want to use maple 13 for calculating mean value theorem for differentiable function  f:=piecewise(-2≤x≤0, -x² ,0≤x≤2,x²) 

on the interval [-2,2]. But an error occured, that is,"function must be continuous".Any help will be appreciated.

REGARDS

 Yegan

please is there any one can help me to find a solution of a sytem of 3 non linear equations each with 3 variable and with more than 30 unknown coefficients

this is the system

solve({EEE_x(x, y, z) = 0, EEE_y(x, y, z) = 0, EEE_z(x, y, z) = 0}, {x, y, z})

where x,y,r are the unknowns

and the three equations are simply the partial derivative with respect to x,y and z repectively

EEE_x(x,y,z):=(&DifferentialD;)/(&DifferentialD; x) EE(x,y,z)

EEE_y(x,y,z):=(&DifferentialD;)/(&DifferentialD; y) EE(x,y,z)

EEE_z(x,y,z):=(&DifferentialD;)/(&DifferentialD;z)EE(x,y,z)

the main equation is EE where (it has 3 variables and more than 30 qunknowns coefficients

(x, y, z) ->

1
----------------------------------------------------------------
2
/ 2 2 2\
hh \ii + jj x + ll z + mm y + 100. y + nn y z + oo x + pp z /

/ 2 2 2 2 3 2
\p z y + q z y + l z x + g z x + o z y + n z x + m y x

2 2 2 3 2 2 2
+ j y x + k y x + i z y + d z y + f z x + h z y

2 2 4 3 2 3
+ e y x + u z y x + v z y x + a + b x + c x + r x + s z

2 2 4 3 4 \
+ t z + bb z + cc y + dd y + ee y + ff y + gg z + aa x/

1.  a procedure quadsumstats whose input is an integer n. This procedure should return a list of length 

n whose kth  entry is the number of solutions to
x^2 + y^2 = k 
for
1 <= k and k <= n

I am sort of confused as to how to construct that list of length n and how to obtain integer solutions to the equation in maple.

2.

a procedure firstCount(k) that finds the first integer
n
with
k
representations as
"x^2+y^2= n." What does it mean for an integer to have k representations?

I am trying to use maple to plot a poincare section for the following Hamiltonain:

H:=(1/8)*(p1^2+16*p2^2-4*p1*p2*cos(q1-q2))/(3+sin(q1-q2)^2) - cos(q2)-8cos(q1)

I've been using Maples built in command as follows:

poincare(H, t=-5000..5000, ics, stepsize = 0.1, iterations = 1, scene = [q2,p2]);

for a given set of initial conditions ics. My problem is however I need to restrict the plotting value to p1>0, as otherwise i seem to get two overlapping maps as seen below:



How exactly can i do this?
Thanks
Connor

This is a question I have also submitted to the technical support, I am worried that it is a bit too technical for them, however:)

I am debugging a C program which links against the OpenMaple API library (under Linux and with Maple 17 and 18). I am using valgrind memcheck, because I am experiencing strange behavior which could be due to writes beyond allocated blocks of memory.  

The first thing which jumps to my eye, are many errors of the types

Use of uninitialised value of size (4/8/16)

Invalid read of size (4/8/16)

Conditional jump or move depends on uninitialised value(s) 

The same errors are also printed when I use the examples that ship with Maple. For instance, I compile "simple.c" with

gcc  -Wl,--no-as-needed -lmaple -lmaplec -lrt -L /usr/lib -L $MAPLEDIR/bin.X86_64_LINUX -I $MAPLEDIR/extern/include -o simple simple.c

and run valgrind as

valgrind --tool=memcheck --error-limit=no --log-file=memcheck.log ./simple

memcheck.txt 

Some, but not all, of the errors occur in __intel_sse2_strcpy or __intel_sse2_strlen. Furthermore, according to valgrind there are definite memory leaks. which appear in the library. 

Practically this makes it hard for me to identify my potential own errors. I am a bit surprised to see so many warnings because I tend to fix my own programs until memcheck does not print these anymore (before I give it away at least). The question is: Can I consider these errors as safe to ignore? How would I distinguish real errors which may appear in my application?

Hello!

To get the phase portrait, I did this

Eq1:=diff(x(t),t)=1-d*x(t)-x(t)*v(t);

Eq2:=diff(y(t),t)=-a*y(t)+x(t)*v(t)-y(t)*w(t);

Eq3:=diff(z(t),t)=-b*z(t)+y(t)*w(t);

Eq4:=diff(v(t),t)=-p*v(t)+y(t);

Eq5:=diff(w(t),t)=-q*w(t)+c*z(t);

d:=0.012:a:=0.93:c:=40:b:=5.6:p:=5.6:q:=5.6:
ics:=x(0)=5,y(0)=1,z(0)=2,v(0)=0.5,w(0)=4;

eq1:=1-d*x-x*v;
eq2:=-a*y+x*v-y*w;
eq3:=-b*z+y*w;
p:=5.6:
eq4:=-p*v+y;
eq5:=-q*w+c*z;

solve({eq1=0,eq2=0,eq3=0,eq4=0,eq5=0},{x,y,z,v,w});

initialset:={seq(seq(seq(seq(seq([x(0)=a1,y(0)=a2,z(0)=a3,v(0)=a4,w(0)=a5],a1=0..5),a2=0..1),a3=0..2),a4=0..0.5),a5=0..4)}:

A:=DEplot([Eq1,Eq2,Eq3,Eq4,Eq5],[x(t),y(t),z(t),v(t),w(t)],t=0..140,x=5..7,y=0..2,initialset,stepsize=0.01,color=blue,linecolor=magenta,arrows=medium,axes=boxed):

Error, (in DEtools/DEplot/WhichPlot) More than two dependent variables - please indicate the desired scene.

I want phase portrait projected on x − y plane.

Any comments?

Hello,

How can I write a code to calculate the Rieman Sum for  y=x^1/2 [0..4] using 

left hand rule and 100 subdivision.

Thank you 

Hi guys, 

       After my calculations, I got the Lagrange function having three gereralized co-ordinates namely u,v and w. Hence, according to theory it should give three equations of motion. 

      I have trouble finding equations of motion from following Lagrange function(L = T -U), Can anyone guide me with this?

      Moreover, kinetic energy and strain energy equations are in the form of double integral!! 

Hello,

i need help!!

write a procedure for the taylor series sin (x) and plot it in the range (-2pi to 2 pi)

use 20 term iterations in the taylor series approximation.

Thank you very much for your help.... 

I'd like to plot the differences between terms in a sequence of vectors. Each difference term should start at the end of the last difference term, so that if I was to plot the actual term, the vector would meet at the end of the difference term. The sequence is limited in length and stored as a list.

Ex:

Suppose I have  . The difference terms would thus be . The first difference term would be plotted be plotted from <0,0>, the second starting at <1,1>, the third at <2,2>.

I can compute the difference terms, but I am not sure how to make the plot I desire. Is arrow(...) the answer, somehow?