Gonzalo Garcia

140 Reputation

6 Badges

14 years, 55 days

MaplePrimes Activity

These are replies submitted by

@Carl Love Thanks.

The article, in this case, is available for free download.


Yes, as you say, I see that some variables are not used...Maybe a "debug" that I forgot to delete? Please, find attached the "corrected" file. Even for the simple function (t-0.5)^2, I see that the approximation is not "very fine"....

I reiterate my gratitude for your useful help.


@Carl Love @mmcdara @vv If it's not a nuisance, someone might check the above attached file, please? I "suspect" that I am doing something wrong ...


I have tried to implement the algorithm descrbied in https://www.sciencedirect.com/science/article/pii/S0898122105002397


Many thanks.

@mmcdara  Thanks for your comments.

I agree with you. Maybe the best is to implemennt a "more actual" algorithm, such as this by J.M. Calcin




I have trieed to encode it into Maple, but I am not sure tat that the code be correct



@vv Thnaks for your comment.  

@Axel Vogt Thanks for yuor suggestion, I know the "DirectSearch" package, but I think that it works (very fine) for determinstic optimization, and not for stochastic optimization as the indicated algorithm.

@vv Many thanks. This works!!

@Carl Love Thanks you! Your code seems work fine!

@vv Can you attach the Maple 2015 file, please?


@vv Sorry, the code returns me the error

Error, invalid proc termination

What do you think about this code?  GKLS.mw


I reiterate my gratitude for your help.

@Carl Love  Many thanks for your reply,

I have coded this GKLS.mw  to compute the centers (I think that the calculus of the radii is easy).

However, CreaParam(2, 15, .6)


Return me 2 centers instead of 15. Can you help me please?


@vv  Thanks for your answer, it sounds good!


Sorry, but I do not understand yuor code at all. The radius R1 has to be equal to 0.6, and the center X1 not necessarily be the (0,0).

Also, your code return me an error: 


Error, invalid sequence

(for r in R do   do v:=Xrnd(r) until check(v,r); S(++NS) := [v, r]; od)


Many thanks for your answer, but it is seems that PHCPack is not free: http://homepages.math.uic.edu/~jan/PHCpack/node11.html#SECTION000110000000000000000





@acer thanks for your comments.

@tomleslie thanks for your reply! Really, it is easier.

@dharr thanks for your comments!

1 2 3 4 5 Page 2 of 5