Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Hello!
Try to consider the following system of differential equations:

sys := diff(U(ksi), ksi) = (y(ksi)-(1.5*(1+0.01*ksi))*U(ksi)/ksi)/(1.5+0.015*ksi+0.002*ksi^2),
diff(y(ksi), ksi) = U(ksi)*(0.002+(1.5*(1+0.01*ksi))*(0.002*ksi^2)/(ksi^2*(1.5+0.015*ksi+0.002*ksi^2))-19.3^2)-y(ksi)*(0.002*ksi^2)/(ksi*(1.5+0.015*ksi+0.002*ksi^2))

with the boundary conditions: cond:= y(0.8) = 0, y(1) = 0

And Maple gives me zero solution for this system , i mean U(from 0.8  to 1) = 0 and y(from 0.8 to 1) = 0

How can i get some other solutions?


P.S. i need numerical solution => dsolve( [sys,cond],type= numeric)

Please, i need some help.

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 ?

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

 

 

 

 

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?

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)

 

 

In this work the theme of vector analysis shown from a computational point of view; this being a very important role in the engineering component; in civil and mechanical special it is why, using the scientific software Maple develops interactive solutions for long processes through MapleCloud calculations. At present the majority of professors / researchers perform static classes open source leaves; so that our students learn and memorize commands, thus generating more time learning in the area. Loading Bookseller VectorCalculus develop topics: vector algebra, differential operators, conservative fields, etc. Maplesoft making processes provide immediate calculations long operation Embedded Components displayed in line with MapleNet integrations. Today our future engineers to design solutions and will be launched in the cloud thus being a process with global qualification in the specialty. Significantly Maple is a scientific software which allows the researcher to design their own innovations and not use themes for their manufacturers.

 

III_CRF_2015.pdf

CRF_2015.mw

 

L.AraujoC.

 

 

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 recently got myself a 64 bit computer and have noticed that I cannot use the option compile=true in dsolve/numeric. Take the following simple example:

dsolve({diff(x(t),t)=x(t),x(0)=1},numeric,compile=true);
Error, (in dsolve/numeric/SC/preproc) unable to compile (rc=1), please try again, and if that fails verify your Windows compiler installation

I'm using Windows 10, but had the same problem with Windows 8.1 on the same machine.
The Compiler:-Compile examples in the help page all work.

What do I have to do to make the option compile=true work in dsolve/numeric?

You may safely assume that I don't know any technicalities about these things.

hi, I am tyying to solve this equation but there is arising an error, plz help me,
VIM.3rd_order.mw

First 1221 1222 1223 1224 1225 1226 1227 Last Page 1223 of 2228