Ronan

Ask a Question

Create a Post

14 Badges

Member for: 13 years, 43 days

Location: East Grinstead, United Kingdom

Contact Ronan

MaplePrimes Activity

These are questions asked by Ronan

Problem with MutableSet...

Asked: Ronan 1331 Product: Maple 2018

March 03 2019

2 4

.I am testing out the MutableSet command, There seems to be a problem with the number of emements that can be operated on. When N =100000 things are ok. When N=1000000 (or there abouts) the process doesn't terminate last line. Uorks ok for sets. Using Maple 2018.2 on Win10 64bit.

(1)

$A := {seq(i^2, i = 1 .. N)}$

(2)

(3)

Download Mutable_set_issue.mw

efficient operation on sets...

Asked: Ronan 1331 Product: Maple 2018

February 26 2019

3 3

I have 2 questions here

1st I want to extend this set Ug3 := {1, 19, 23, 29} by repeatively subtracting 6 fron its elements. first pass should give

Ug3 := (Ug3-~ 6) union Ug3  Ug3 := {-5, 1, 13, 17, 19, 23, 29}

Now I dont want to add any elenments less than 6 so -5 need to go

Second pass produces

Ug3 := (Ug3-~ 6) union Ug3                Ug3 := {-11, -5, 1, 7, 11, 13, 17, 19, 23, 29}

Obiviously no elements less than 6 to be added on second pass either.

At this point there is no further elements that could be added to the set.

I am looking for an efficient method. I can use For loops and if then checks. The sets grow in length quickly. Repetition lengths are 6, 30 210, 2310,... I use sets because the numbers get ordered automaticially. Are sets efficient or should I use some form of table?

How do I add and multiply combinations o...

Asked: Ronan 1331 Product: Maple 2018

February 05 2019

2 1

This should be simple but I am missing something. I want to sum the product of combinations of numbers.

I know I can assign values to a0...a[n-1]. But that isn't nice.

(n,k):= (5,4);
add(mul~(combinat:-choose([a(||)[](0..n-1)], k))) ; =a0*a1*a2*a3+a0*a1*a2*a4+a0*a1*a3*a4+a0*a2*a3*a4+a1*a2*a3*a4 which is correct.
add(mul~(combinat:-choose([(0..n-1)], k))) ; #should be 24 but I get 0

How to save procedures to a library?...

Asked: Ronan 1331 Product: Maple 2018

October 25 2018

1 2

I have asked this before but am still confused. I have a half-dozen procedures I want to save. Don't want them in a module/package. I am using windows 10.

I just can't get the syntax correct on this.

Obviously after saving restart and load to test.

libname;
       "C:\Program Files\Maple 2018\lib", 

         "C:\Users\Ronan\maple\toolbox\CodeBuilder\lib", 

         "C:\Users\Ronan\maple\toolbox\OEIS\lib", 

         "C:\Users\Ronan\maple\toolbox\personal\lib", 

         "C:\Users\Ronan\maple\toolbox\UTF8\lib"
libdir := "C:/Users/Ronan/maple/toolbox/personal/lib";
     libdir := "C:/Users/Ronan/maple/toolbox/personal/lib"
NULL;

LibraryTools:-Save(Pedal, cat(kernelopts(homedir), "/maple/toolbox/personal/lib/Pedal.mpl"));
Error, (in LibraryTools:-Save) could not open `C:\Users\Ronan/maple/toolbox/personal/lib/Pedal.mpl\Pedal.m` for writing

Is the for..do loop in my procedure caus...

Asked: Ronan 1331 Product: Maple 2018

October 22 2018

3 9

I trying to get a proc to work. In the 1st half of the document, I derive four sections for a pedal curve and plot them fine. In the second I try doing it with a procedure called Pedal. Am having problems. I think it might be something to do with the for.. do loop.

C := (x^2+y^2+12*x+9)^2-4*(2*x+3)^3

(1)

[(-x^2-12*x-9+2*(8*x^3+36*x^2+54*x+27)^(1/2))^(1/2), -(-x^2-12*x-9+2*(8*x^3+36*x^2+54*x+27)^(1/2))^(1/2), (-x^2-12*x-9-2*(8*x^3+36*x^2+54*x+27)^(1/2))^(1/2), -(-x^2-12*x-9-2*(8*x^3+36*x^2+54*x+27)^(1/2))^(1/2)]

(2)

X^4+4*X^3*r+2*X^2*Y^2+4*X^2*Y*s+6*X^2*r^2+2*X^2*s^2+4*X*Y^2*r+8*X*Y*r*s+4*X*r^3+4*X*r*s^2+Y^4+4*Y^3*s+2*Y^2*r^2+6*Y^2*s^2+4*Y*r^2*s+4*Y*s^3+r^4+2*r^2*s^2+s^4-8*X^3-24*X^2*r+24*X*Y^2+48*X*Y*s-24*X*r^2+24*X*s^2+24*Y^2*r+48*Y*r*s-8*r^3+24*r*s^2+18*X^2+36*X*r+18*Y^2+36*Y*s+18*r^2+18*s^2-27

(3)

(4)

(5)

xp1 := 1; yp1 := 2; lprp1 := -coeff(ltg1, y)*x+coeff(ltg1, x)*y+K1; for i to nops(sol1) do s := eval(sol1[i], x = r); print("s ", s); K1 := solve(eval(lprp1, [x = xp1, y = yp1]), K1); print("K1 ", K1); sol2 := solve([ltg1, lprp1], [x, y]); Pedal || i := [rhs(sol2[1, 1]), rhs(sol2[1, 2])]; p9 || i := plot([op(Pedal || i), r = -20 .. 20], colour = "MediumSeaGreen"); unassign('K1', 'sol2') end do

"K1 ", -4*(-r^2-12*r-9-2*(8*r^3+36*r^2+54*r+27)^(1/2))^(3/2)-4*r^2*(-r^2-12*r-9-2*(8*r^3+36*r^2+54*r+27)^(1/2))^(1/2)+16*(8*r^3+36*r^2+54*r+27)^(1/2)*r-48*r*(-r^2-12*r-9-2*(8*r^3+36*r^2+54*r+27)^(1/2))^(1/2)+192*r^2+96*(8*r^3+36*r^2+54*r+27)^(1/2)-36*(-r^2-12*r-9-2*(8*r^3+36*r^2+54*r+27)^(1/2))^(1/2)+576*r+432

(6)

(7)

This I can't get to work. All 4 sections of the pedal curve are the same and should not be.

Parse:-ConvertTo1D, "`%1` is not a module or member", RonanRoutines

"K||i", -8*r^3+(-4*(-r^2-12*r-9-2*(8*r^3+36*r^2+54*r+27)^(1/2))^(1/2)+48)*r^2-8*(-(-r^2-12*r-9-2*(8*r^3+36*r^2+54*r+27)^(1/2))^(1/2)-3)^2*r-4*(-r^2-12*r-9-2*(8*r^3+36*r^2+54*r+27)^(1/2))^(3/2)+48*r^2+576*r+432+96*(8*r^3+36*r^2+54*r+27)^(1/2)-36*(-r^2-12*r-9-2*(8*r^3+36*r^2+54*r+27)^(1/2))^(1/2)