Unanswered Questions

This page lists MaplePrimes questions that have not yet received an answer

Hey, this is not the I've had this encounter. I want to open this saved document but when I open it and Maple starts up it just hits me with "A problem was encountered while opening the workbook. Database is not opened". How can I get to open it properly and see my math notes?

How can this be prevented?

When I try to upload the file in this message it says "Cant open a null file"

Any help?

 

Jacob

 

plz code and sole that this integral!!!!!!!!!!

!!!

Hi, so I have few problems here. I need to;

Create a MxN Matrix/Lattice, where N and M can be any positive integer, that contains a random selection of -1/1s at each entry.

Need to sum every entry, then multiply by -1 to find “H”.

Need to multiply each neighbour to find its bond energy, so if it’s the same you get 1 else -1, but only its direct neighbours once, so if it was a 2x2 matrix there would be only 4 values and then sum them.

I don’t seem to be able to set up the code so that it does it for any size matrix, as I only know how to write it out basic for a 2x2. Also, not so important, but I wanted to know if could create a loop that would find every iteration of possible setups i.e. for a 2x1 you can have 1.1, -1.1, 1.-1, -1.-1. And then give the solutions outlined earlier for each of the the possibilities [There being 2MxN ]

Cheers in advance.

In response to Markiyan Hirnyks reference to Mathematica's FindDistribution command http://www.mapleprimes.com/questions/219668-Determine-The-Distribution-Function.  It would be nice to have something similar coded in Maple. 

solve Derivative in maple

see in photo

 

What do people find a good screen size is? 

Personally I find two monitors would be ideal for using Maple and especially Maplesim.  I'm going to be adding a second monitor, probably a 20".  Maybe 24" is better?  I think 20" would be enough.  Just wondering how many monitors people out there work with and their monitor sizes.

Dear all

I have n, p two integers greater than one

I would like to minimize the following real number by a fixed positive constant if its possible

(n/(n+p))^(n/p)-(n/(n+p))^((n+p)/p)

Many thanks

Dear all,

I am creating an animation, and I was wondering if I can add a multiplier to my equation in a specific range.  So starting from z:=0.2 add a multiplier (z+1). The code I have so far is added. Does anyone know a code for this?

Kind regards

restart; 
with(plots); 
a := -1/2; b := 1/2; c := -2; d := 2; n := 20; 
g := proc (x) options operator, arrow; value(Int(sigma(t), t = 0 .. x)) end proc; 
sigma := proc (z) options operator, arrow; 2*sqrt(2*h^2-4*z^2)*z/h^2 end proc; 
h := i/n; 
for i to n do 
an2[i] := plot(sigma(z), z = -(1/2)*h .. (1/2)*h, view = [a .. b, c .. d], color = AQUAMARINE); 
an3[i] := plot(2*g(x), x = 0 .. (1/2)*h, view = [a .. b, c .. d], color = RED) 
end do; 
p := plots[display]([seq(an2[i], i = 1 .. n)], insequence = true); 
q := plots[display]([seq(an3[i], i = 1 .. n)], insequence = true); display(p, q)

 


equation 1 : xi+1=xi− (f·gy−fy·g)/(fx ·gy −fy ·gx)
equation 2: yi+1=yi− (fx·g−f·gx)/(fx·gy-fy·gx)

My quesiton are, deriving equations (1) and (2) above and constructing a single Maple function called newt2d that implements both of these recurrence relation.

I apolgize in advance if I don't write my question correctly.  This is my first time posting a question. 

> restart;

with(plots);

pr := .72; p := 0; n := [.5, 1, 1.5]; s := 0; a := .2; b := 0; L := [red, blue, green]; l := 0; k := 1;

for j to nops(n) do R1 := 2*n[j]/(1+n[j]); R2 := 2*p/(1+n);

sol1 := dsolve([diff(diff(diff(f(eta), eta), eta), eta)+f(eta)*(diff(diff(f(eta), eta), eta))+R1*(1-(diff(f(eta), eta))^2) = 0, diff(diff(theta(eta), eta), eta)+pr*k*f(eta)*(diff(theta(eta), eta))+R2*pr*k*(diff(f(eta), eta))*theta(eta)+(2*(a*(diff(f(eta), eta))+b*theta(eta)))/(1+n[j]) = 0, f(0) = 1, (D(f))(0) = b*((D@@2)(f))(0), (D(f))(1.8) = 0, theta(0) = 1+s*(D(theta))(0), theta(1.8) = 1], numeric, method = bvp);

fplt[j] := plots[odeplot](sol1, [eta, diff(diff(f(eta), eta), eta)], color = L[j], axes = boxed); tplt[j] := plots[odeplot](sol1, [[eta, theta(eta)]], color = L[j], axes = normal) end do; plots:-display([seq(fplt[j], j = 1 .. nops(n))]);

plots:-display([seq(tplt[j], j = 1 .. nops(n))]);

 

staganation_point1.mw
 

can we chage the axis sir ?? like  f'' vs eta to f'' vs lambda.

``

restart

l := 1:

1

 

1.5

 

.5

 

[blue, green, red, yellow]

(1)

``

for j to nops(p) do R1 := 2*n/(n+1); R2 := 2*p[j]/(n+1); R3 := 2/(n+1); sol1 := dsolve([diff(diff(diff(f(eta), eta), eta), eta)+f(eta)*(diff(diff(f(eta), eta), eta))+R1*(1-(diff(f(eta), eta))^2)-M*(diff(f(eta), eta)) = 0, diff(diff(theta(eta), eta), eta)+pr*f(eta)*(diff(theta(eta), eta))-R2*pr*(diff(f(eta), eta))*theta(eta)+R3*(A*(diff(f(eta), eta))+B*theta(eta)) = 0, f(0) = 1, (D(f))(0) = L+b*((D@@2)(f))(0), (D(f))(7) = 1, theta(0) = 1+s*(D(theta))(0), theta(7) = 0], numeric, method = bvp); plots[odeplot](sol1, [eta, ((D@@2)(f))(eta)], color = red); fplt[j] := plots[odeplot](sol1, [eta, f(eta)], color = K[j], axes = boxed); tplt[j] := plots[odeplot](sol1, [[eta, theta(eta)]], color = K[j], axes = normal); fplt[j] := plots[odeplot](sol1, [eta, diff(f(eta), eta)], color = K[j], axes = boxed) end do:

 

 

plots:-display([seq(fplt[j], j = 1 .. nops(n))]);

 

sol1(0)

sol1(0)

(2)

sol1(.1)

[eta = .1, f(eta) = 1.05958091104306206, diff(f(eta), eta) = .643210624614908300, diff(diff(f(eta), eta), eta) = .881482678165403044, theta(eta) = .623284688471349546, diff(theta(eta), eta) = -.578039450700496560]

(3)

sol1(.2)

[eta = .2, f(eta) = 1.12800452943200891, diff(f(eta), eta) = .722346769554029544, diff(diff(f(eta), eta), eta) = .706526135439307756, theta(eta) = .568123251856343492, diff(theta(eta), eta) = -.525530979400813946]

(4)

sol1(.3)

[eta = .3, f(eta) = 1.20351830506746449, diff(f(eta), eta) = .785511903074783246, diff(diff(f(eta), eta), eta) = .561442941644520022, theta(eta) = .518103974464032668, diff(theta(eta), eta) = -.475257424178228970]

(5)

sol1(.4)

[eta = .4, f(eta) = 1.28466826824405134, diff(f(eta), eta) = .835505660630676662, diff(diff(f(eta), eta), eta) = .442470716586289281, theta(eta) = .472985640642506311, diff(theta(eta), eta) = -.427567049032814172]

(6)

sol1(.5)

[eta = .5, f(eta) = 1.37026161183094430, diff(f(eta), eta) = .874752886901313142, diff(diff(f(eta), eta), eta) = .345911467377074400, theta(eta) = .432494259338694842, diff(theta(eta), eta) = -.382764248064397461]

(7)

sol1(.6)

[eta = .6, f(eta) = 1.36678221814533528, diff(f(eta), eta) = .771028661281065508, diff(diff(f(eta), eta), eta) = .407805382194403932, theta(eta) = .876413930517023876, diff(theta(eta), eta) = -.197648778495384870]

(8)

sol1(2)

[eta = 2., f(eta) = 2.66120522956795602, diff(f(eta), eta) = .991532161353848585, diff(diff(f(eta), eta), eta) = 0.251405465681268682e-1, theta(eta) = .635967939441598018, diff(theta(eta), eta) = -.144641270049362308]

(9)

``

``

``

 

``


 

Download staganation_point1.mw

 

 

 

 

 

Hello every one,

I'm using Maple to create a compiled program of my model Maplesim. I had a error with "getcompiledProc" command, which is :

Error, (in GetCompiledProc) non-numeric initial condition for `Main.RAB.value`(t): Float(undefined)

the "Main.RAB.value(t)" represent a probe used to visualize an output of my system.

Does anybody had a such error ?
Thank you for your responses.

Larbi ANIBA

I am writing a package to define a few Records and collect some quantities and formulae to do some arithmetic with.

At some point I started running into a severe problem: Each time I want to use one of the exports form the module, the Maple kernel dies with a "Stack limit exceeded" message. The issue is: This seems to happen only on Maple 2016. I had some issues in Maple 2015 but that got sorted by rewriting some stuff. I never had any issue like that with Maple 15. I went as far as commenting out most of my .mapleinit code to make sure that did not cause issues. The code works as intended in Maple 15 and in Maple 2015.

I am posting this here and am asking others to see if they can reproduce this. Needless to say I will SCR this. But I actually need to get this to work.

The error can be triggered by simply entering

Booster;

after loading the package. Needless to say, it should just echo the name "Booster", to get more one would need to enter eval(Booster) (due to the last-name evaluation rules).

Thanks,

M.D.

BooMBA.mw

 

Hello everybody,

 

I want to solve this pde.the desire solution is V(r,z). three boundar conditions are written that two of them are related to rhe radial and one is related to the longitudinal coordinate.

I attached the solution for you. but this solution is derived by Matlab. Now, I just want to resolve it by Maple, but I couldn't reach it. Please let me know the correct way asap.

 

Thanks a lot.PDE.mw
 

NULL

NULL

restart 

Q := diff(V(r, z), r, r)+(diff(V(r, z), r))/r-V(r, z)/r^2+diff(V(r, z), z, z)+C/r = 0

diff(diff(V(r, z), r), r)+(diff(V(r, z), r))/r-V(r, z)/r^2+diff(diff(V(r, z), z), z)+C/r = 0

(1)

NULL

NULL

pdsolve(Q)

PDESolStruc(V(r, z) = _F1(r)*_F2(z)-(1/2)*(_C1/r+_C2*r+_C3*r*ln(r))*C/_C3, [{diff(diff(_F1(r), r), r) = _F1(r)*_c[1]+(-(diff(_F1(r), r))*r+_F1(r))/r^2, diff(diff(_F2(z), z), z) = -_F2(z)*_c[1]}])

(2)

NULL

 

NULL

NULL

``

``


 

Download PDE.mwPDE.mw

 

 

Howdy!!

Just thought I'd put a quick one in. I'm trying to solve the Schrodinger equation given as:

schro := {diff(psi(x), x, x)-((b^2-a*(2*p+3))*x^2+2*ab*x^4+a^2*x^6-energy)*psi(x) = 0};

under the constraint

energy=(2p+1)b

I'm trying to plot this for a range of values of the parameters a,b and fixed p=0,1.

I'm aiming to get a solution of the form

psi(x)=(x^p) exp((-a(x^4))/4)-(b(x^2)/2)

and plotting it.

Thanks in advance

Recently I have used a FilterForge filter to generate hyperbolic tiles from images (Lua code with Möbius transform etc. available). If an input image (1 see below) is selected a tile (2) is generated if the hyperbolic parameters (5,3,2) were chosen.  How could this be made with Maple?

 

First 176 177 178 179 180 181 182 Last Page 178 of 362