MaplePrimes Questions

I want to calculate this limit

limit((b^(d+1)-1)/((b-1)*b^d), d = infinity)

but it always return limit((b^(d+1)-1)/((b-1)*b^d), d = infinity)

 

Is it becaus it cant recognize what b symbol is ?

Another easy question from a newbie:

Say I have a procedure which, along the its execution, needs to pass a parameter to a different procedure (which, in this case, would be a subroutine) to be elaborated and returned to go on.

This process may happen a few times during the execution and the passing parameter, which may have the same name but different value, needs to be worked on by the subrooutine.

My question is: how do I invoke a procedure from a parent procedure passing one or more parameters and have them returned?

Thank you

Hi,

I am trying to use implicit plot. The plot is OK but I want to put labels such as what Latex produces :

$\frac{\Omega}{\omega_n}$ for x axis and $a_0 \mathrm{(m)}$ for y axis

How can I apply this in my maple code as below:

plot1:=implicitplot(a3, Omega_r=1.5..2.5, a=0.00000001..0.1, labeldirections=[horizontal, horizontal], axes=boxed, labels=["W/w_n",typeset("a_0 (m)")], labelfont=[SYMBOL]):

What I have put as bold does not work for me, it is making everything in Greek :) . I want combination of Greek and math.

 

Thasnks,

Bahareh

 

 

 

 

I am trying to design a question, where studenst have to find the equation of a linear function given its graph. For this to work, I have to be able to draw gridlines, which for some reason is not possible in Maple TA with the ordinary plot function. I know it should be possible using an applet; given the problems this is going to create with all the different browsers, however, this is not a viable solution. I have found another solution which is unreasonably complex but should work. The problem is that sometimes it doesn't work. When I press Refresh algorithm preview in the Question Designer I get "Broken Maple plot. Verify your plot statement" roughly every third time I refresh the preview. Any help would be greatly appreciated.

 

The algorithm is here:

$a1=range(-4,4,1);
$a2=range(1,4,1);
condition:not(eq($a1,$a2));
$px=range(0,5,1);
$y0=range(0,5,1);
$x0=if(gt($a1,0),$px,-$px);
$f=$a1/$a2*(x-$x0)+$y0;
$G=plotmaple("
 plots[display](
  [plots[coordplot](cartesian, [-10 .. 10, -10 .. 10], grid = [11, 11], color = grey),
   plot($f, x = -10 .. 10, y = -10 .. 10, color = [green, blue, red], thickness = 2, tickmarks = [[1], [1]], labels = [``, ``]),
   plots[textplot]([0, 10, 'y'], align = {above, right}),
   plots[textplot]([10, 0, 'x'], align = {above, right})
  ],
  labels = [``, ``], axes = normal, view = [-10 .. 10, -10 .. 10]),
 plotoptions='width=350,height=350'");
$g=gcd($a1,$a2);
$at=$a1/$g;
$an=$a2/$g;
$a=if(eq($an,1),$at,$at/$an);
$ans=$a*x+($y0-$a*$x0);

Hi!

This question is related to http://www.mapleprimes.com/questions/204419-Derivatives-Of-Splines-Are-Not-Defined and http://www.mapleprimes.com/questions/42114-Problem-With-Spline-Integrating , however I have not been able to apply the solutions given there to my problem.

I have a set of points given by

and certain function value points given by

where e1 is the function I am approximating.


Using

I come up with my piecewise function.

When I do diff(e4,x), however, the points at the nodes show "float(undefined) x=0.2..."(the node).

As it turns out, the value of the derivative on the left of the node is not equal to that of the right side by a factor of 10^(-7), in other words, numerically unimportant but high enough for maple to realise it is not the same number. How could I tell Maple that I am happy choosing, for example, the value given by the function on right side of the node?

I welcome any suggestions.

Many thanks in advance.

Every time I try to write a procedure I get stuck.

This time is no different:

restart;

global a:=0.081819221, PI:=3.1415926535897932384626433832795;
Ecce:=proc(lt)
lat:=lt*(PI/180);
b:=((1-a*sin(lat))/(1+a*sin(lat)))^(a/2));
t:=ln((tan(PI/4)+lat/2));
d:=3437.7468*t*b;
return d;
end proc;

I digit the following to get a result
Ecce(45.2112);

and this is what I get (in blue)
Ecce(45.2112)

Every single time. I can never have a procedure that works right away. It's getting on my nerves

I want to extract all the symbols and last trade from 

http://finance.yahoo.com/q/cp?s=%5EIXIC

I have tried:

status, data, headers := HTTP:-Get("http://finance.yahoo.com/q/cp?s=%5EIXIC");
data;

but it just gives me a bunch of jiberish! Any ideas?

The equation tan(y) = 2*tan(x) defines y implicitly as a function of x.  Well, perphas "defines" is too strong a word, since there are multiple solutions for y.  However, if I am not mistaken, there exists a unique continuous solution y(x) that goes through the origin, that is, y(0)=0, and is defined for all x.

Question 1: How do we plot the graph of y(x)?

I have a roundabout solution as follows.  Differentiate the equation tan(y(x)) = 2*tan(x) with respect to x and arrive at a first order differential equation in y(x).  Solve the differential equation with the initial condition y(0)=0.  Surprisingly, Maple obtains an explicit solution:

which we can plot:

plot(rhs(%), x=0..2*Pi);

Question 2: Is there a neat way of getting that solution with algebra only, without appealing to differential equations?

 

please i would love make a 3d plot of this Bessel function expression against the parameter p and T[0] using maple (p^(((2+deltaT[0])/(2)))*(C[1]*BesselJ((delta)/(-2)*sqrt(T[0]-4 T[g]),(-((T[g])/(T^(2)[0])*gamma*G*delta*sqrt(alpha*k))/(p))+C[2]*BesselY((delta)/(-2)*sqrt(T[0]-4 T[g]),(-((T[g])/(T^(2)[0])*gamma*G*delta*sqrt(alpha*k))/(p)) where [delta = 0.2e-2, T[g] = (1.1, .615, .48, .2962), k = (1.2, 1.3, 1.4, 1.5), G = .2, gamma = .2, alpha = 5.36, C[1] = 500, C[2] = 100,p = (.22, .23, .24, .25), T[0] = (3.666667, 2.307692, 1.714286, 1.377778)

 

we can get the continued fraction of cot(x) with the command

> convert(cot(x), confrac, x, 6)

the result is

(1+x^2/(-3+x^2/(5-(1/7)*x^2)))/x

 convert(cot(x), confrac, x, 6)

how could I convert this result to a form like

1/x-(1/10)*x+49*x/(20*x^2-210)

 

 

I have some excel data which I need to fit in the formula:

0.5a*erfc(0.5*2^0.5*(-x+m1/s1)) + (0.5-0.5a)*erfc(0.5*2^0.5*(-x+m2/s2))

In this formula the coefficient m ans s are mean and standard deviation and a is the weigth of 2 peak in a cummulative gaussian distribution.

I fitted this (I will insert the maple file) and found:

-48736.43 erfc(-0.707x + 0.707) + 48736.43 erfc(0.707x + 0.707)

Can this be correct?

And how do I find the values of m1, s1, m2, s2 and a?

data.xlsx

data_fit.mw

I want to solve numerically the nonlinear pde:

 

u_x+u_t - (u_{xt})^2 = u(x,t)

 

which method do you propose me to use with maple? (I don't mine about which boundary conditions to be used here).

 

Hello.
Have a question for use Physics[TransformCoordinates] . For example, in a Cartesian coordinate system is an arbitrary tensor. As you know, in this case the covariant, contravariant and mixed components are the same. We have a coordinate transformation such as transforming our Cartesian coordinate system in the oblique coordinate system. In this simple example, correctly calculated the covariant, contravariant and mixed components (structure: covariant and contravariant), but other mixed components (structure: contravariant and covariant) are calculated is not correct. I checked by hand.



restart

with(Physics):

Setup(mathematicalnotation = true):

ds := dx[1]^2+dx[2]^2+dx[3]^2:

Setup(coordinates = (X = [x[1], x[2], x[3]]), dimension = 3, metric = ds, spacetimeindices = lowercaselatin, quiet):

g_[]:

A[a, b] = Matrix(3, 3, {(1, 1) = 2, (1, 2) = 1, (1, 3) = 3, (2, 1) = 2, (2, 2) = 3, (2, 3) = 4, (3, 1) = 1, (3, 2) = 2, (3, 3) = 1})

A[a, b] = (Matrix(3, 3, {(1, 1) = 2, (1, 2) = 1, (1, 3) = 3, (2, 1) = 2, (2, 2) = 3, (2, 3) = 4, (3, 1) = 1, (3, 2) = 2, (3, 3) = 1}))

(1)

Define(%):

`Defined objects with tensor properties`

(2)

A[]

A[a, b] = (Matrix(3, 3, {(1, 1) = 2, (1, 2) = 1, (1, 3) = 3, (2, 1) = 2, (2, 2) = 3, (2, 3) = 4, (3, 1) = 1, (3, 2) = 2, (3, 3) = 1}))

(3)

A[`~`]

A[`~a`, `~b`] = (Matrix(3, 3, {(1, 1) = 2, (1, 2) = 1, (1, 3) = 3, (2, 1) = 2, (2, 2) = 3, (2, 3) = 4, (3, 1) = 1, (3, 2) = 2, (3, 3) = 1}))

(4)

A[`~a`,b,matrix]

A[`~a`, b] = (Matrix(3, 3, {(1, 1) = 2, (1, 2) = 1, (1, 3) = 3, (2, 1) = 2, (2, 2) = 3, (2, 3) = 4, (3, 1) = 1, (3, 2) = 2, (3, 3) = 1}))

(5)

A[`a`,~b,matrix]

A[a, `~b`] = (Matrix(3, 3, {(1, 1) = 2, (1, 2) = 1, (1, 3) = 3, (2, 1) = 2, (2, 2) = 3, (2, 3) = 4, (3, 1) = 1, (3, 2) = 2, (3, 3) = 1}))

(6)

[y[1] = x[1]-x[2], y[2] = x[2]-x[3], y[3] = x[3]]

[y[1] = x[1]-x[2], y[2] = x[2]-x[3], y[3] = x[3]]

(7)

solve((7), {x[1], x[2], x[3]})

{x[1] = y[1]+y[2]+y[3], x[2] = y[2]+y[3], x[3] = y[3]}

(8)

OK

B[a,b] = TransformCoordinates((8), A[a, b], [y[1], y[2], y[3]], [x[1], x[2], x[3]], simplifier = `@`(`simplify/size`, simplify))

B[a, b] = (Matrix(3, 3, {(1, 1) = 2, (1, 2) = 3, (1, 3) = 6, (2, 1) = 4, (2, 2) = 8, (2, 3) = 15, (3, 1) = 5, (3, 2) = 11, (3, 3) = 19}))

(9)

OK

C[a,b] = TransformCoordinates((8), A[~a,~b], [y[1], y[2], y[3]], [x[1], x[2], x[3]], simplifier = `@`(`simplify/size`, simplify))

C[a, b] = (Matrix(3, 3, {(1, 1) = 2, (1, 2) = -1, (1, 3) = -1, (2, 1) = 0, (2, 2) = -2, (2, 3) = 3, (3, 1) = -1, (3, 2) = 1, (3, 3) = 1}))

(10)

OK

D[a,b] = TransformCoordinates((8), A[a,~b], [y[1], y[2], y[3]], [x[1], x[2], x[3]], simplifier = `@`(`simplify/size`, simplify))NULL

D[a, b] = Matrix(%id = 4452149890)

(11)

Bug

E[a, b] = TransformCoordinates((8), A[~a,b], [y[1], y[2], y[3]], [x[1], x[2], x[3]], simplifier = `@`(`simplify/size`, simplify))NULLNULL

E[a, b] = Matrix(%id = 4452139458)

(12)

should be

LinearAlgebra:-Transpose(rhs((12)))

Matrix(3, 3, {(1, 1) = 0, (1, 2) = -2, (1, 3) = -3, (2, 1) = 1, (2, 2) = 2, (2, 3) = 5, (3, 1) = 1, (3, 2) = 3, (3, 3) = 4})

(13)

``



Download Transformation_tensor_components.mw

I have a function fitted and found the coeffiecients, the coeffiecient values are found with data from maple

the function is: a*(x+b/c)

my result is 5x+15

How do I get the values of a, b and c? The value of a is 5 that I know, but how does I get that maple gives me all the values of a,b  and c? The problem is that I get one value while i need the values of a, b and c.

Thanks

First 1242 1243 1244 1245 1246 1247 1248 Last Page 1244 of 2429