## 45 Reputation

3 years, 173 days

## Thank you @tomleslie for the e...

Thank you @tomleslie for the expalanation. Is there any way to avoid the jumps on pi (or -pi) in maple ?

## Dear @Kitonum ,  thank yo...

Dear @Kitonum ,
thank you for the help.

## Hey @mmcdara ! Thank you for ...

Hey @mmcdara !

Thank you for the trick !

## Hey @MapleMathMatt, Thank you very ...

Hey @MapleMathMatt,
Thank you very much for your quick reply. Yes you are right; the areas should be cumulative. I figured it out a half hour after I posted the question. That is why I forgot to comment back.  Thank you again for the help!

## I agree with you!...

Thank you @tomleslie  for your comments. I am trying to build up a code and at some point; I need to compare the discretized functions (solutions as vectors) at each iteration with the ones of the previous iteration. I tried to present my problem that way to avoid possible distractions from the rest of the code. However, I think that the example provided by @dharr  is clearer.

## it works...

Thank you @dharr for your help.

## No...

No worries @mmcdara ! my message was not that clear either,  and thank you very  much for the further details which I am really interested in.

## Thank you @Carl Love for the c...

Thank you @Carl Love for the clarification. It is more clear now than before.

## Sorry for the confusion...

@mmcdara

The range 0.5..1 for x[1] and x[2] and  -1..-0.5 for x[3] and x[4] was just an experiment for me to restrict the range of the solution in fsolve. However, as I mentioned before, I am interested in those that lie on the lines Im x_i = ±π/2.

## Clarification...

Yes, you are right. I should mention that I was talking about one solution among many… next time I will be more careful. However, as I mentioned before, I am interested in those that lie on the lines Im x_i = ±π/2.

## I guess that I found some solutions...

Thank you very much for your detailed answer. I guess that these solutions (for a=0.8) :

x[1] = -.43334870250217230+1.5707963267948966*I, x[2] = .43334870250217230+1.5707963267948966*I, x[3] = -.43334870250217230-1.5707963267948966*I, x[4] = .43334870250217230-1.5707963267948966*I

verified the criteria. Any comments will be very appreciated.
if the solutions are valid, how can I guess the solutions for N=8 ?

## @tomleslie  I am sorry! I did not ...

I am sorry! I did not insert the link. Now, I edited my questions. Thanks for letting me know!

## @Mariusz Iwaniuk  for my point of ...

for my point of view (numerically), It is better to keep  k(x) in the momentum space because you will need to deal with the convolution numerically, then you will have just to multiply both functions (k and log(1+f_1)) instead of dealing with convolution' integrals.

## @Mariusz Iwaniuk  it works, thank y...

@Mariusz Iwaniuk  it works, thank you !

## Thank you very much  @vv ...

Thank you very much  @vv  for your help. when I tried to run it, it gives me the following  error "Error, invalid sum/difference" . The cursor is pointing on

`[crt++"here"]`

 1 2 Page 1 of 2
﻿