Ronan

1022 Reputation

14 Badges

13 years, 160 days
East Grinstead, United Kingdom

MaplePrimes Activity


These are replies submitted by Ronan

@dharr Thank you.This may interest you. I was following a series of math videos form N. J. Wildberger on YouTube, on the origins of Boolean Algebra.  It rapidly leads into an efficient  canonical form of Boole's algebra. This is the link to the first in the series .A brief history of logic: Aristotle and deduction | Math Foundations 251 | NJ Wildberger - YouTube

It is great to hear that one's question is of use to another.

@Carl Love My sloppy code didn't help. I knew I should swap jm and h before I first posted.

 im := 1/2*im; istep := 1/2*istep; jm := 2*jm; h := 1/2*h; 

 I then assumed because V1 transformed back to V0 when the transform was applied to it that the code must be correct. Well done to @dharr for catching that.

Also I hadn't checked the n=3 case against the original matrix which would have caught the problem.

@Carl Love Compare the timings

Tonight's

n = 25:
memory used=4.00GiB, alloc change=256.00MiB, cpu time=63.14s, real time=57.23s, gc time=13.16s
V0 = RTABLE(18446746400855522110, (100110111010010011100011010100

  10hellip... 33554400 row vector entries not shown), 

  Vector[row], Handle)
V1 = RTABLE(18446746400855518854, (111000110010101001111000111101

  10hellip... 33554400 row vector entries not shown), 

  Vector[row], Handle)

Previous Nights

n = 25:
memory used=0.50GiB, alloc change=0.50GiB, cpu time=140.00ms, real time=154.00ms, gc time=0ns
V0 = RTABLE(18446746173536082150, (100110111010010011100011010100

  10hellip... 33554400 row vector entries not shown), 

  Vector[row], Handle)
V1 = RTABLE(18446746173536074798, (100110111010010011100011010100

  10hellip... 33554400 row vector entries not shown), 

  Vector[row], Handle)

 

@Carl Love That code is amazing. I had no idea something could be improved so vastly. Thank you for the clarifiaction on the 1D input.

@Carl Love Ok. I'm using Maple 2019.2. Is a different syntax required because 1D still did not like the +=. In the loops you have "to n do" Should I change these to "for i to n do" etc? It getting late here 2am. I will give it another shot tomorrow evening.

LAM_BM_Transform-1D.mw

@dharr That is very fast.

Can a compiled procedure be saved to disk? And then used in apackage?

 

@Carl Love  I tried the code but it does not transform the vector, I could be doing something wrong. I corrected "istart+=h" because was getting an error message about unterminated loop, to "istart:=istart+h". I ran this for n=3 so the vector only has 8 elements.

why do you have to make a copy of the vector inside the procedure? 

LAM:-Copy(2, V, 1..-1, R, 1..nv);

I tried copy and pasting the whole code but would probably only introduce more errors tidying it up.

LAM_BM_Transform.mw

 

 

Hello,

 Nice work. At some point I will have to study up on the Draghilev method. I gather it applies to dynamic systems.

Here is an SG platform I post a few years ago. It only deals with a static version of the reverse kinematics. It would be interesting to see the method applied to it.

https://www.mapleprimes.com/posts/208301-Stewart-Gough-Reverse-Kinematics

 


Hope this helps

restart

``

T := sqrt(6)*sqrt((n-1)*(-15*k^2*n+15*k^2-30*k*n+30*k-15*n+15+sqrt(1920*n*q*sigma^6-1920*q*sigma^6+225*k^4*n^2+5760*Q*n*sigma^3-450*k^4*n+900*k^3*n^2-5760*Q*sigma^3+225*k^4-1800*k^3*n+1350*k^2*n^2+900*k^3-2700*k^2*n+900*k*n^2+1350*k^2-1800*k*n+225*n^2+900*k-450*n+225)))/(12*sigma^2*(n-1)), -sqrt(6)*sqrt((n-1)*(-15*k^2*n+15*k^2-30*k*n+30*k-15*n+15+sqrt(1920*n*q*sigma^6-1920*q*sigma^6+225*k^4*n^2+5760*Q*n*sigma^3-450*k^4*n+900*k^3*n^2-5760*Q*sigma^3+225*k^4-1800*k^3*n+1350*k^2*n^2+900*k^3-2700*k^2*n+900*k*n^2+1350*k^2-1800*k*n+225*n^2+900*k-450*n+225)))/(12*sigma^2*(n-1)), sqrt(-(6*(n-1))*(15*k^2*n-15*k^2+30*k*n+sqrt(1920*n*q*sigma^6-1920*q*sigma^6+225*k^4*n^2+5760*Q*n*sigma^3-450*k^4*n+900*k^3*n^2-5760*Q*sigma^3+225*k^4-1800*k^3*n+1350*k^2*n^2+900*k^3-2700*k^2*n+900*k*n^2+1350*k^2-1800*k*n+225*n^2+900*k-450*n+225)-30*k+15*n-15))/(12*sigma^2*(n-1)), -sqrt(-(6*(n-1))*(15*k^2*n-15*k^2+30*k*n+sqrt(1920*n*q*sigma^6-1920*q*sigma^6+225*k^4*n^2+5760*Q*n*sigma^3-450*k^4*n+900*k^3*n^2-5760*Q*sigma^3+225*k^4-1800*k^3*n+1350*k^2*n^2+900*k^3-2700*k^2*n+900*k*n^2+1350*k^2-1800*k*n+225*n^2+900*k-450*n+225)-30*k+15*n-15))/(12*sigma^2*(n-1))

(1/12)*6^(1/2)*((n-1)*(-15*k^2*n+15*k^2-30*k*n+30*k-15*n+15+(1920*n*q*sigma^6-1920*q*sigma^6+225*k^4*n^2+5760*Q*n*sigma^3-450*k^4*n+900*k^3*n^2-5760*Q*sigma^3+225*k^4-1800*k^3*n+1350*k^2*n^2+900*k^3-2700*k^2*n+900*k*n^2+1350*k^2-1800*k*n+225*n^2+900*k-450*n+225)^(1/2)))^(1/2)/(sigma^2*(n-1)), -(1/12)*6^(1/2)*((n-1)*(-15*k^2*n+15*k^2-30*k*n+30*k-15*n+15+(1920*n*q*sigma^6-1920*q*sigma^6+225*k^4*n^2+5760*Q*n*sigma^3-450*k^4*n+900*k^3*n^2-5760*Q*sigma^3+225*k^4-1800*k^3*n+1350*k^2*n^2+900*k^3-2700*k^2*n+900*k*n^2+1350*k^2-1800*k*n+225*n^2+900*k-450*n+225)^(1/2)))^(1/2)/(sigma^2*(n-1)), (1/12)*(-6*(n-1)*(15*k^2*n-15*k^2+30*k*n+(1920*n*q*sigma^6-1920*q*sigma^6+225*k^4*n^2+5760*Q*n*sigma^3-450*k^4*n+900*k^3*n^2-5760*Q*sigma^3+225*k^4-1800*k^3*n+1350*k^2*n^2+900*k^3-2700*k^2*n+900*k*n^2+1350*k^2-1800*k*n+225*n^2+900*k-450*n+225)^(1/2)-30*k+15*n-15))^(1/2)/(sigma^2*(n-1)), -(1/12)*(-6*(n-1)*(15*k^2*n-15*k^2+30*k*n+(1920*n*q*sigma^6-1920*q*sigma^6+225*k^4*n^2+5760*Q*n*sigma^3-450*k^4*n+900*k^3*n^2-5760*Q*sigma^3+225*k^4-1800*k^3*n+1350*k^2*n^2+900*k^3-2700*k^2*n+900*k*n^2+1350*k^2-1800*k*n+225*n^2+900*k-450*n+225)^(1/2)-30*k+15*n-15))^(1/2)/(sigma^2*(n-1))

(1)

``

solve(T[1], Q)

-(1/3)*q*sigma^3

(2)

``

solve(T[1], q)

-3*Q/sigma^3

(3)

solve(T[1], sigma)

(-3*Q*q^2)^(1/3)/q, -(1/2)*(-3*Q*q^2)^(1/3)/q-((1/2)*I)*3^(1/2)*(-3*Q*q^2)^(1/3)/q, -(1/2)*(-3*Q*q^2)^(1/3)/q+((1/2)*I)*3^(1/2)*(-3*Q*q^2)^(1/3)/q

(4)

solve(T[1], k)

``

solve(T[1], n)*`~nosolution`

``


 

Download MP_Help_C.mw

@AHSAN 

Hello, I only palyed with MapleSim a few years ago, so can't answer your specific question. Try this link.

https://www.youtube.com/c/maplesoft/playlists  There are several series on MapleSim.

Hope this helps.

@janhardo with the abs in the Maple equation the answer came out negative. I dont know how abs reacts in the fsolve command.

I notice your arctan functions have Pi/90. Have you set GE to working degrees? Radians are better for for general problem solving. That is throwing your equation for solving in Maple.

@janhardo Thanks for that. I checked it out. All I had to do was remove the absolute value from the Area equation after pasting it into Maple.

restart

`` ``

A2 := (-R^2+4*r^2)^(1/2)*R*(-1/2)+arctan((-R^2+4*r^2)^(1/2)/R)*R^2-arctan((-R^2+4*r^2)^(1/2)*R/(R^2-2*r^2))*r^2

-(1/2)*(-R^2+4*r^2)^(1/2)*R+arctan((-R^2+4*r^2)^(1/2)/R)*R^2-arctan((-R^2+4*r^2)^(1/2)*R/(R^2-2*r^2))*r^2

(1)

``

r := 1

1

(2)

A2

-(1/2)*(-R^2+4)^(1/2)*R+arctan((-R^2+4)^(1/2)/R)*R^2-arctan((-R^2+4)^(1/2)*R/(R^2-2))

(3)

fsolve(A2 = (1/2)*Pi*r^2, R)

1.158728473

(4)

``

``


 

Download GE_Circle_puzzle.mw

 

Hello,

A bit of a side line question here. I got a copy of GE the other day as it is on special offer .How did you make the overlapping area of the circles?

First 10 11 12 13 14 15 16 Last Page 12 of 24