MaplePrimes Questions

Hey, I have a problem loading one of my maple files as the headline states.

I have seen a sugestion on how to fix it:

http://www.mapleprimes.com/questions/129377-There-Were-Problems-During-The-Loading

but I do not know how to use it proberly, and can therefore not fix my problem.

What do I do?

Elektriske_felter.mw


 

I want to generate table of 2 , how can i do this?

Dear Community,

Would someone have a good and easy to understand/implement description of the Den Iseger algorithm for the numerical inversion of Laplace transform? Even better if someone would have a Maple script to do it, that would be superb.

Tx in advance,

best regards

Andras

hi...

how i can save results in text file?

for example in attached code maple i want save results for different ''x'' and fy(x), wherer x is between 0 to 1 i.e. (0,fy(0)) , (0.1,fy(0.1)) ,(0.2,fy(0.2)),....

in text file two column should be create.one of them for ''x'' and the onther for ,fy(x)

thanks...

Numerical_.mw

restart; E := 0.169e12; mu := 0.658e11; hl2 := 4; D1 := 1; n := 3; `αn` := 0.; beta := 16.474184; xi := 1.5; lambda := .1; dsys5 := {D1*(diff(y(x), x, x, x, x)) = `αn`/(1-y(x))^n-beta*(lambda*cosh(xi*(1-y(x)))-1/2*(lambda^2+1))/sinh(xi*(1-y(x)))^2, y(0) = 0, (D(y))(0) = 0, ((D@@2)(y))(1) = 0, ((D@@3)(y))(1) = 0}; dsol5 := dsolve(dsys5, 'maxmesh' = 900, numeric, output = listprocedure); fy := eval(y(x), dsol5); fy(1)

HFloat(0.3189919226404612)

(1)

NULL

NULL

NULL

NULL



Download Numerical_.mw

 

Hi

may every one help to me for dsolve this differentia1l equation?

error:

Error, (in dsolve/numeric/bvp) singularity encountered

Turbulent2-kw.mw

dsol1 := dsolve({diff(theta(eta), eta, eta)-3*Omega(eta)*(F(eta)*(diff(theta(eta), eta))-theta(eta)*(diff(F(eta), eta)))/(2*K(eta))+((diff(K(eta), eta))/K(eta)-(diff(Omega(eta), eta))/Omega(eta))*(diff(theta(eta), eta)) = 0, diff(F(eta), eta, eta, eta)+Omega(eta)*(3*F(eta)*(diff(F(eta), eta, eta))-(diff(F(eta), eta))^2)/(2*K(eta))+((diff(K(eta), eta))/K(eta)-(diff(Omega(eta), eta))/Omega(eta))*(diff(F(eta), eta, eta))+Omega(eta)/K(eta) = 0, diff(K(eta), eta, eta)+Omega(eta)*(1.5*F(eta)*(diff(K(eta), eta))-K(eta)*(diff(F(eta), eta)))/K(eta)+((diff(K(eta), eta))/K(eta)-(diff(Omega(eta), eta))/Omega(eta))*(diff(K(eta), eta))+(diff(F(eta), eta, eta))^2-Omega(eta)^2 = 0, diff(Omega(eta), eta, eta)+Omega(eta)*(3*F(eta)*(diff(Omega(eta), eta))+Omega(eta)*(diff(F(eta), eta)))/(2*K(eta))+((diff(K(eta), eta))/K(eta)-(diff(Omega(eta), eta))/Omega(eta))*(diff(Omega(eta), eta))+Omega(eta)*(diff(F(eta), eta, eta))^2/K(eta)-Omega(eta)^3/K(eta) = 0, F(0) = 0, K(0) = 0, Omega(0) = 0., theta(0) = 1, theta(1) = 0, (D(F))(0) = 0, (D(K))(1) = 0, (D(Omega))(1) = 0, ((D@@2)(F))(1) = 0}, numeric, method = bvp[middefer], output = listprocedure, initmesh = 512)

Error, (in dsolve/numeric/bvp) singularity encountered

 

NULL

plots[odeplot](dsol1, [(D(F))(eta), eta])

Error, (in plots/odeplot) input is not a valid dsolve/numeric solution

 

``



Download Turbulent2-kw.mw

 

 

 

Hello:

I want to make a package or module. Not sure. I cant get my head around what to do. I have loads of formulas I would like to be able to use without going back through work sheets. At present just want to keep it simple.

Have loaded a sample work sheet to use

Intend to place package in Maple2016 library folder.

 

restart

Library Folder

C:\Program Files\Maple 2016\lib

``

``

with(LinearAlgebra):

``

``

Points

Point2 := proc (x1, y1) options operator, arrow; [x1, y1] end proc

proc (x1, y1) options operator, arrow; [x1, y1] end proc

(1)

Point3 := proc (x1, y1, z1) options operator, arrow; [x1, y1, z1] end proc

proc (x1, y1, z1) options operator, arrow; [x1, y1, z1] end proc

(2)

a1 := Point3(x1, y1, z1)

[x1, y1, z1]

(3)

a2 := Point3(x2, y2, z2)

[x2, y2, z2]

(4)

a3 := Point3(x3, y3, z3)

[x3, y3, z3]

(5)

Line through 2 points

LinePts := proc (p1, p2) options operator, arrow; (p2[2]-p1[2])*x+(p1[1]-p2[1])*y-p2[2]*p1[1]+p1[2]*p2[1] end proc

proc (p1, p2) options operator, arrow; (p2[2]-p1[2])*x+(p1[1]-p2[1])*y-p2[2]*p1[1]+p1[2]*p2[1] end proc

(6)

``

Line

Line := proc (a, b, c) options operator, arrow; a*x+b*y+c end proc

proc (a, b, c) options operator, arrow; a*x+b*y+c end proc

(7)

Lines s

LnPrll := proc (l1, l2) options operator, arrow; is(coeff(l1, x)*coeff(l2, y)-coeff(l2, x)*coeff(l1, y) = 0) end proc

proc (l1, l2) options operator, arrow; is(coeff(l1, x)*coeff(l2, y)-coeff(l2, x)*coeff(l1, y) = 0) end proc

(8)

Lines t

LnPrpnd := proc (l1, l2) options operator, arrow; is(coeff(l1, x)*coeff(l2, x)+coeff(l1, y)*coeff(l2, y) = 0) end proc

proc (l1, l2) options operator, arrow; is(coeff(l1, x)*coeff(l2, x)+coeff(l1, y)*coeff(l2, y) = 0) end proc

(9)

``

``

Quadrance between 2 points (2D)

Quadrance := proc (a1, a2) options operator, arrow; (a1[1]-a2[1])^2+(a1[2]-a2[2])^2 end proc

proc (a1, a2) options operator, arrow; (a1[1]-a2[1])^2+(a1[2]-a2[2])^2 end proc

(10)

``

Quadrance between 2 points (3D)

Qd3 := proc (a1, a2) options operator, arrow; (a1[1]-a2[1])^2+(a1[2]-a2[2])^2+(a1[3]-a2[3])^2 end proc

proc (a1, a2) options operator, arrow; (a1[1]-a2[1])^2+(a1[2]-a2[2])^2+(a1[3]-a2[3])^2 end proc

(11)

Triple Quad Formula (3 colinear points)

TQF := proc (Q1, Q2, Q3) options operator, arrow; (Q1+Q2+Q3)^2 = 2*Q1^2+2*Q2^2+2*Q3^2 end proc

proc (Q1, Q2, Q3) options operator, arrow; (Q1+Q2+Q3)^2 = 2*Q1^2+2*Q2^2+2*Q3^2 end proc

(12)

Spread*PolynomialsNULL

SpreadPoly := proc (n, s, rn) with(orthopoly); s-factor(1/2-(1/2)*T(n, 1-2*rn)) end proc

proc (n, s, rn) with(orthopoly); s-factor(1/2-(1/2)*T(n, 1-2*rn)) end proc

(13)

SpreadPoly(3, s, r)

s-r*(4*r-3)^2

(14)

``

 

Download Testpackage.mw

Very new to using Maple or coding in general. Ran into a problem that I couldn't figure it out.

Any help would be greatly appreciated.

Has anyone been able to do multivariate partial fraction decomposition in maple (here is a paper introducing the idea https://arxiv.org/pdf/1206.4740.pdf)

I often find maple generating complicated rational functions that it would be nice to visualise in other ways

Here is an example of such a function if anyone wants to have a play:

 

(a*x^3+b*x*y^2+a*x*y+b*y^2)

/(a*x^3+a*b*x*y^2+a*b*y*x^2+b*y^3)

 

Hi There:

i have big expressions with many variables of the form (say)  Q__T(t). I would like to simplify my visual parsing of them.

I was wondering if there is a way to tell maple to print particular symbols/variables in specific colors (something akin to what the Physics package does). I.e. I would imagine saying something like:

setup_print_color(Q__T(t), "red");

and then every time Q__T(t) would come up, it would be printed in my desired color.

is this easily possible? I've looked through docs and this site, but can't find any info.

thanks

 

 

Hello guys

 

I am trying to solve a non-linear 2nd order diffential equation in order to describe the turbulence phenomena for tip vortices .

The initial conditions that i have are

b(0)=1

b(8)=exp(.1842068074*ln(-(2.*(-1.+b(x)))/(2.*b(8)-7.))+9.670857391)

Moreover i know that b from 0 to 1 equals one and then it  increases propably parabololoid(from the experimntal measurements)with respect o Re(Reynolds)

I tried to solve the eqution using maple solve de interactivelly but doesnt seem to work

Do you have any suggestions?

Thank you

giorgio_3.mw

 

I attempted to use convolution in AudioTools to convolve 2 vectors with arbitrary values, but could not since the operation is expecting numeric values.  Can this be done?    Is there not a convolution operastor in the LinearAlgebra package?  See attached:

convolutionWithAu[1].mw

I have a base class with a lot of assignements. The classes that I let my students use are all child classes of that. I used to be able to hide all the assignements except those, that the students should actually work on, but now I can no longer do that. The base class is shared with several colleges and not "owned" by me, but that didn't ust to be an issue. Is there a solution to this problem?

I am new to Maple.

If i have a list (or a table) A:=[[1,2],[3,4],[5,6]

how can i split it into two lists B:=[1,3,5] and C:=[2,4,6] ? (all the 1st numbers in one list, and all 2nd numbers in another list).

Hello guys

Before starting, i must tell you that english is not my mother tongue, so sorry for my spelling errors :)

 

I have a work to do that ask me to solve the equation : e^((x))-x-5 with the Fixed Point Iteration method. (starting value=2 ; shutoff parameter = 10^(-8) and maximum number of iterations = 100).

I already try and at this time, I did that :

http://imgur.com/PGvmusk

But it doesn't seem to work. 

I am not very good with Maple...

I need some help :)

Thanks !

I'm trying to plot the lines that are based on eigenvectors of the matrix. (Functions r1(t),r2(t), q(t) are defined earlier. I can upload file if necessary)

with(LinearAlgebra):
x1:=t->-r1(t)*cos(Pi/2-q(t)/2):
y1:=t->-r1(t)*sin(Pi/2-q(t)/2):
z1:=t->0:
x2:=t->r2(t)*cos(Pi/2-q(t)/2):
y2:=t->-r2(t)*sin(Pi/2-q(t)/2):
z2:=t->0:
Ixx:=t->m*(y1(t)^2+z1(t)^2)+m*(y2(t)^2+z2(t)^2):
Iyy:=t->m*(x1(t)^2+z1(t)^2)+m*(x2(t)^2+z2(t)^2):
Izz:=t->m*(x1(t)^2+y1(t)^2)+m*(x2(t)^2+y2(t)^2):
Ixy:=t->-m*x1(t)*y1(t)-m*x2(t)*y2(t):
Ixz:=t->-m*x1(t)*z1(t)-m*x2(t)*z2(t):
Iyz:=t->-m*y1(t)*z1(t)-m*y2(t)*z2(t):

InertiaTensor:=t-><<Ixx(t),Ixy(t),Ixz(t)>|<Ixy(t),Iyy(t),Iyz(t)>|<Ixz(t),Iyz(t),Izz(t)>>:

 Here I define the matrix as a function of time. 

EV:=t->Eigenvectors(InertiaTensor(t));
V1:=t->Column(Re(EV(t)[2]),1);
V2:=t->Column(Re(EV(t)[2]),2);

Producing two eigenvectors also as functions of time.

InertiaAxis1 := t->plot(x, V1(t)[2]*x/V1(t)[1], x = -2 .. 2);
InertiaAxis2 := t->plot(x, V2(t)[2]*x/V2(t)[1], x = -2 .. 2);

InertiaAxis1(1);

And here I've got an error:

Error, (in plot) unexpected options: [-HFloat(7.378655652881484e-6)*x, x = -2 .. 2]

I tried to set "datatype=float"-option in plot function, to evaluate the coefficient by eval-function, it wouldn't help. What does this error mean exactly? 

First 1066 1067 1068 1069 1070 1071 1072 Last Page 1068 of 2428