Advanced Search

Ronan

1506 Reputation

16 Badges

14 years, 224 days
East Grinstead, United Kingdom
Contact Ronan

MaplePrimes Activity


These are answers submitted by Ronan

No offical way, but...

Ronan 1506
July 23 2021
1 0

Maple doesn't seem to have that formatting at present. You could try and use a single cell table scaled to the correct width.

Not really a plesent solution I admit.

 

Table_Approx_A4.mw

R(...) is a function...

Ronan 1506
July 17 2021
5 0

You need a space or a multiplication * between the R and bracket. Then it solves.

solve(R*(sigma + mu)/nu = (N*b + R*sigma)/(nu + mu), R)

R(...) is a function.

Hope this helps

Sum a finite amount of terms as an appro...

Ronan 1506
December 13 2020
1 1

Would this be acceptable as the sum converges
My internet is faulty at present so I can't get the document to display.

restart;


fd := j -> 256/3*j^5*(j - 1)^(2*j - 4)/(j + 1)^(2*j + 4);
S := x -> sum(fd(n)*ln(1 - 1/n^2), n = 2 .. x);
plot(S(x), x = 2 .. 20);
fd(10);
                             "(->)"
evalf(S(10));
for x from 2 by 5 to 200 do
    x, evalf(S(x));
end do;

 

Download MP_sum.mw

 

 

DirectSearch:-SolveEquations...

Ronan 1506
November 29 2020
0 0

Maybe I miss understand what you are after.
 

> 

restart; with(DirectSearch)
> 

E := (-7.29511346879067*10^9*b-12693.49752)*((D@@2)(phi))(x)+(5.36829162574305*b+0.3487242314e-3)*((D@@4)(phi))(x)+(2.02316788462733*10^17*b+4.107030812*10^12)*phi(x)+1.74615339888993*10^17*b

(-7295113469.*b-12693.49752)*((D@@2)(phi))(x)+(5.36829162574305*b+0.3487242314e-3)*((D@@4)(phi))(x)+(0.2023167885e18*b+0.4107030812e13)*phi(x)+0.1746153399e18*b

 (1)
> 

dsolve(E, phi(x))
> 

ivp := {E, phi(-0.406e-3) = 1, phi(0.406e-3) = 1, (D(phi))(-0.406e-3) = 0, (D(phi))(0.406e-3) = 0}

{(-7295113469.*b-12693.49752)*((D@@2)(phi))(x)+(5.36829162574305*b+0.3487242314e-3)*((D@@4)(phi))(x)+(0.2023167885e18*b+0.4107030812e13)*phi(x)+0.1746153399e18*b, phi(-0.406e-3) = 1, phi(0.406e-3) = 1, (D(phi))(-0.406e-3) = 0, (D(phi))(0.406e-3) = 0}

 (2)
> 

N := dsolve(ivp, phi(x))
> 

A := simplify(rhs(dsolve(ivp, phi(x))))
> 

W := evalf(subs(phi(x) = simplify(rhs(dsolve(ivp, phi(x)))), (2.74505742817409*10^8+3.18054074382591*10^8*phi(x))^2-(1.63833629390402*(2.74505742817409*10^8+3.18054074382591*10^8*phi(x)))*(diff(phi(x), x, x))+3*(6.98727355657655*10^7+32282.4885498329*(diff(phi(x), x)))^2+2.68414581287153*(diff(phi(x), x, x))^2 = 0.9e17))
> 

 

> 

with(DirectSearch)

[BoundedObjective, CompromiseProgramming, DataFit, ExponentialWeightedSum, GlobalOptima, GlobalSearch, Minimax, ModifiedTchebycheff, Search, SolveEquations, WeightedProduct, WeightedSum]

 (3)
> 

Digits := 30

30

 (4)
> 

SolveEquations(eval(W,{x=-0.0000039}))

Warning, complex or non-numeric value encountered; trying to find a feasible point

  

[0.147107922240261059617024e-2, Vector(1, {(1) = 0.383546505968e-1}), [b = 0.127348584174732672649668035297e-4], 661]

 (5)
> 

SolveEquations(eval(W, {x = -0.3e-5}))

Warning, complex or non-numeric value encountered; trying to find a feasible point

  

[0.419387893880979948510784e-2, Vector[column](%id = 18446745433166527174), [b = 0.127348584175363152130255238047e-4], 791]

 (6)
> 

SolveEquations(eval(W, {x = -0.2e-5}))

Warning, complex or non-numeric value encountered; trying to find a feasible point

  

[0.483612833457601000000e-5, Vector[column](%id = 18446745433826848886), [b = 0.127348584176092433330136395418e-4], 510]

 (7)
> 

SolveEquations(eval(W, {x = -0.1e-5}))

Warning, complex or non-numeric value encountered; trying to find a feasible point

  

[0.230450895918287510431716e-2, Vector[column](%id = 18446745433229447406), [b = 0.127348584176852498189365168024e-4], 666]

 (8)
> 

SolveEquations(eval(W, {x = 0}))

Warning, complex or non-numeric value encountered; trying to find a feasible point

  

[.19220551392974036042704144, Vector[column](%id = 18446745433149170022), [b = 0.127348584177643908884594179237e-4], 605]

 (9)
> 

SolveEquations(eval(W, {x = 0.39e-5}))

Warning, complex or non-numeric value encountered; trying to find a feasible point

  

[0.2802257031156986788225e-4, Vector[column](%id = 18446745433231180254), [b = 0.127348584181041268501301157811e-4], 763]

 (10)
> 

``


 

Download MZ2-1.mw

solve...

Ronan 1506
September 24 2020
0 2

Use solve

T has 4 different value.

solve(T[1], Q); gives  -q*sigma^3/3

You get the same answer for T_2, T_3, T_4


                            

 

This Thread should help you....

Ronan 1506
September 10 2020
1 0

 

https://www.mapleprimes.com/posts/212577-Scalable-Subplots-Within-Figure

Possibly check save large calculations...

Ronan 1506
January 01 2020
1 1

 

Might be worth checking this setting.

Shift + F5...

Ronan 1506
September 06 2019
1 2

Shif +F5 stops executable math from running the blue box turns grey. You can use it inside test too. Also RMB to untick executible maths works. When in a line of test and you want to enter an equation click on math button enter equation then RMB untick executible math then click on text button again.
 

restart

``

P = 2*x+5

P = 2*x+5

 (1)

P = 2*x+5

``

Test section P = 2*x+5  now more text
Shift Return for this line.

``


 

Download Text_and_equations.mw

Here are a couple...

Ronan 1506
August 31 2019
1 1

https://www.mapleprimes.com/questions/201298-How-To-Animate-A-Bouncing-Ball-On-A-3d-Surface

or another

https://www.mapleprimes.com/questions/224449-How-To-Model-A-Bouncing-Ball-Through

 

Just Search for Ball.

 

 

Windows 10 File History...

Ronan 1506
June 01 2019
0 0

If you are using Windows 10 turn File History. That keeps a running backup of all files you work on. You coud set to save  once a hour. It has gotten me out of trouble a few times.

Evaluate - Remove output...

Ronan 1506
May 17 2019
1 1

 

They moved it to Evaluate Remove Output. Ctrl+Shift+R

This is one way of doing it...

Ronan 1506
April 26 2019
1 1

This works.
 

restart

with(plots)

[animate, animate3d, animatecurve, arrow, changecoords, complexplot, complexplot3d, conformal, conformal3d, contourplot, contourplot3d, coordplot, coordplot3d, densityplot, display, dualaxisplot, fieldplot, fieldplot3d, gradplot, gradplot3d, implicitplot, implicitplot3d, inequal, interactive, interactiveparams, intersectplot, listcontplot, listcontplot3d, listdensityplot, listplot, listplot3d, loglogplot, logplot, matrixplot, multiple, odeplot, pareto, plotcompare, pointplot, pointplot3d, polarplot, polygonplot, polygonplot3d, polyhedra_supported, polyhedraplot, rootlocus, semilogplot, setcolors, setoptions, setoptions3d, shadebetween, spacecurve, sparsematrixplot, surfdata, textplot, textplot3d, tubeplot]

 (1)

c1 := x*y+x+y = 7

x*y+x+y = 7

 (2)

c2 := x^2+x*y+y^2 = 13

x^2+x*y+y^2 = 13

 (3)

implicitplot([c1, c2], x = -5 .. 5, y = -5 .. 5, colour = [blue, green])

  

sol := `assuming`([solve({c1, c2}, [x, y], useassumptions)], [0 < x and x < 4, 0 < y and y < 4])

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

 (4)

``


 

Download This_works.mw

Using affine geometry...

Ronan 1506
January 19 2019
1 0


Hi,

This really is just a more expanded version of vv 5208  answer.

restart

P := [x3, y3, z3]

[x3, y3, z3]

 (1)

A := [x1, y1, z1]

[x1, y1, z1]

 (2)

B := [x2, y2, z2]

[x2, y2, z2]

 (3)

H lies on the line connecting AB, H can be defined as an affine combination of A and B. -∞ < λ < ∞. If 0 < λ < 1 then λ lies between A and B

 

H := simplify(expand((1-lambda)*A+lambda*B))

[(-x1+x2)*lambda+x1, (-y1+y2)*lambda+y1, (-z1+z2)*lambda+z1]

 (4)

``

Get distance between P and H

distPH := sqrt((H[1]-P[1])^2+(H[2]-P[2])^2+(H[3]-P[3])^2)

(((-x1+x2)*lambda+x1-x3)^2+((-y1+y2)*lambda+y1-y3)^2+((-z1+z2)*lambda+z1-z3)^2)^(1/2)

 (5)

Miniminist the distance by differentiating wrt λ. Then solve for λ

mdist := diff(distPH, lambda)

(1/2)*(2*((-x1+x2)*lambda+x1-x3)*(-x1+x2)+2*((-y1+y2)*lambda+y1-y3)*(-y1+y2)+2*((-z1+z2)*lambda+z1-z3)*(-z1+z2))/(((-x1+x2)*lambda+x1-x3)^2+((-y1+y2)*lambda+y1-y3)^2+((-z1+z2)*lambda+z1-z3)^2)^(1/2)

 (6)

lambda := solve(mdist, lambda)

(x1^2-x1*x2-x1*x3+x2*x3+y1^2-y1*y2-y1*y3+y2*y3+z1^2-z1*z2-z1*z3+z2*z3)/(x1^2-2*x1*x2+x2^2+y1^2-2*y1*y2+y2^2+z1^2-2*z1*z2+z2^2)

 (7)

Back substitute. This is the shortest possible distance  between P and H

simplify(distPH)

((x1^2*y2^2-2*x1^2*y2*y3+x1^2*y3^2+x1^2*z2^2-2*x1^2*z2*z3+x1^2*z3^2-2*x1*x2*y1*y2+2*x1*x2*y1*y3+2*x1*x2*y2*y3-2*x1*x2*y3^2-2*x1*x2*z1*z2+2*x1*x2*z1*z3+2*x1*x2*z2*z3-2*x1*x2*z3^2+2*x1*x3*y1*y2-2*x1*x3*y1*y3-2*x1*x3*y2^2+2*x1*x3*y2*y3+2*x1*x3*z1*z2-2*x1*x3*z1*z3-2*x1*x3*z2^2+2*x1*x3*z2*z3+x2^2*y1^2-2*x2^2*y1*y3+x2^2*y3^2+x2^2*z1^2-2*x2^2*z1*z3+x2^2*z3^2-2*x2*x3*y1^2+2*x2*x3*y1*y2+2*x2*x3*y1*y3-2*x2*x3*y2*y3-2*x2*x3*z1^2+2*x2*x3*z1*z2+2*x2*x3*z1*z3-2*x2*x3*z2*z3+x3^2*y1^2-2*x3^2*y1*y2+x3^2*y2^2+x3^2*z1^2-2*x3^2*z1*z2+x3^2*z2^2+y1^2*z2^2-2*y1^2*z2*z3+y1^2*z3^2-2*y1*y2*z1*z2+2*y1*y2*z1*z3+2*y1*y2*z2*z3-2*y1*y2*z3^2+2*y1*y3*z1*z2-2*y1*y3*z1*z3-2*y1*y3*z2^2+2*y1*y3*z2*z3+y2^2*z1^2-2*y2^2*z1*z3+y2^2*z3^2-2*y2*y3*z1^2+2*y2*y3*z1*z2+2*y2*y3*z1*z3-2*y2*y3*z2*z3+y3^2*z1^2-2*y3^2*z1*z2+y3^2*z2^2)/(x1^2-2*x1*x2+x2^2+y1^2-2*y1*y2+y2^2+z1^2-2*z1*z2+z2^2))^(1/2)

 (8)

``

The coordinates to H

H

[(-x1+x2)*(x1^2-x1*x2-x1*x3+x2*x3+y1^2-y1*y2-y1*y3+y2*y3+z1^2-z1*z2-z1*z3+z2*z3)/(x1^2-2*x1*x2+x2^2+y1^2-2*y1*y2+y2^2+z1^2-2*z1*z2+z2^2)+x1, (-y1+y2)*(x1^2-x1*x2-x1*x3+x2*x3+y1^2-y1*y2-y1*y3+y2*y3+z1^2-z1*z2-z1*z3+z2*z3)/(x1^2-2*x1*x2+x2^2+y1^2-2*y1*y2+y2^2+z1^2-2*z1*z2+z2^2)+y1, (-z1+z2)*(x1^2-x1*x2-x1*x3+x2*x3+y1^2-y1*y2-y1*y3+y2*y3+z1^2-z1*z2-z1*z3+z2*z3)/(x1^2-2*x1*x2+x2^2+y1^2-2*y1*y2+y2^2+z1^2-2*z1*z2+z2^2)+z1]

 (9)

``


 

Download Affine_Geometry.mw

Union...

Ronan 1506
December 29 2018
0 0

Try this, I extended your list,

L1 := [{3, 5}, {4, 5}, {4, 8, 9}, {-7, 2, 3}]:
      
L2 := {};  #empty set
                            L2 := {}
for i to nops(L1) do
L2 := L1[i] union L2
end do;
                          L2 := {3, 5}
                        L2 := {3, 4, 5}
                     L2 := {3, 4, 5, 8, 9}
                  L2 := {-7, 2, 3, 4, 5, 8, 9}
 

gridrefine...

Ronan 1506
November 29 2018
1 1

Try   implicitplot(cos(x)*cosh(y) = 1, x = -3 .. 3, y = -5 .. 5, gridrefine=2).

You can set higher values for gridrefine   3,4,5

5 6 7 8 9 Page 7 of 9