Processing math: 100%

First time here? Check out the FAQ!

Ask Your Question

Revision history [back]

click to hide/show revision 1
initial version

As said in my previous comment to @slelievre's answer (see above), the D[0](real_part) expression has no obvious meaning. Assuming that :

  • the author meant t(t)

  • tR

we can replace it by 1 :

sage: with assuming(t, "real"): real_part(t).diff(t)
1

We note also that this curious notation occurs only in the expression D[0](real_part)(1/4*pi). foo noting the original triplet, duly completed by its closing parenthesis, this replacement leads to :

sage: [u.subs(D[0](real_part)(1/4*pi)==1) for u in foo]
[-1/6*sqrt(6)*sqrt(3)*sqrt(2)*t + 1/2*sqrt(2)*e^(1/4*pi),
 1/2*sqrt(2)*e^(1/4*pi),
 e^(1/4*pi)]

It turns out that all the quantities involved are real :

sage: [[v[0].n() for v in u.subs(D[0](real_part)(1/4*pi)==1).coefficients(t)] for u in foo]
[[1.55088319691803, -1.00000000000000], [1.55088319691803], [2.19328005073802]]

The sought curve is therefore a straight line parallel to the x axis, determined by y=1.55 and z=2.19. Materializing this graph is lazily left as an exercise to the reader which is also incited to check his notations...

HTH,

click to hide/show revision 2
No.2 Revision

As said in my previous comment to @slelievre's answer (see above), above or below according to this site's software whims), the D[0](real_part) expression has no obvious meaning. Assuming that :

  • the author meant t(t)

  • tR

we can replace it by 1 :

sage: with assuming(t, "real"): real_part(t).diff(t)
1

We note also that this curious notation occurs only in the expression D[0](real_part)(1/4*pi). foo noting the original triplet, duly completed by its closing parenthesis, this replacement leads to :

sage: [u.subs(D[0](real_part)(1/4*pi)==1) for u in foo]
[-1/6*sqrt(6)*sqrt(3)*sqrt(2)*t + 1/2*sqrt(2)*e^(1/4*pi),
 1/2*sqrt(2)*e^(1/4*pi),
 e^(1/4*pi)]

It turns out that all the quantities involved are real :

sage: [[v[0].n() for v in u.subs(D[0](real_part)(1/4*pi)==1).coefficients(t)] for u in foo]
[[1.55088319691803, -1.00000000000000], [1.55088319691803], [2.19328005073802]]

The sought curve is therefore a straight line parallel to the x axis, determined by y=1.55 and z=2.19. Materializing this graph is lazily left as an exercise to the reader which is also incited to check his notations...

HTH,

click to hide/show revision 3
No.3 Revision

As said in my previous comment to @slelievre's answer (see above or below according to this site's software whims), the D[0](real_part) expression has no obvious meaning. Assuming that :

  • the author meant t(t)

  • tR

we can replace it by 1 :

sage: with assuming(t, "real"): real_part(t).diff(t)
1

We note also that this curious notation occurs only in the expression D[0](real_part)(1/4*pi). foo noting the original triplet, duly completed by its closing parenthesis, this replacement leads to :

sage: [u.subs(D[0](real_part)(1/4*pi)==1) for u in foo]
[-1/6*sqrt(6)*sqrt(3)*sqrt(2)*t + 1/2*sqrt(2)*e^(1/4*pi),
 1/2*sqrt(2)*e^(1/4*pi),
 e^(1/4*pi)]

It turns out that all the quantities involved are real :

sage: [[v[0].n() for v in u.subs(D[0](real_part)(1/4*pi)==1).coefficients(t)] for u in foo]
[[1.55088319691803, -1.00000000000000], [1.55088319691803], [2.19328005073802]]

The sought curve is therefore a straight line parallel to the x axis, determined by y=1.55 and z=2.19. Materializing this graph is lazily left as an exercise to the reader which is also incidentally incited to check his notations...

HTH,