What's the original ODE you are trying to solve? The solution you give
has symbols in it that aren't in your diffeq.
One change I would recommend is to use sympy.pi instead of math.pi.
math.pi is a float approximation to pi, whereas sympy.pi is pi
exactly. Alternatively, you can use w = symbols('omega') if you want
the solution to match the general text version.
Aaron Meurer
On Mon, Jun 27, 2022 at 2:14 PM Federico Manfrin
<federicomanf...@gmail.com> wrote:
>
> Hi, I'm writing a notebook to reproduce what I can see in this wikipedia page:
> https://it.wikipedia.org/wiki/Circuito_RC
>
> The problem is the CAP 4, I'm looking forward to get the result that follow
> that text:
> dove K è una costante. Dunque: (latex folleow)
> {\displaystyle v_{C}(t)=v_{C}(0)e^{-{\frac {t}{\tau }}}+K\cos(\omega t+\theta
> )}
>
> This is my code, running in a colab notebook:
> # define the independent variable ‘t’
> from sympy.abc import t
> import sympy as sp
> import math
> w = 2 * math.pi * 50
>
> # define dependent variable in symbol form
> vc = sp.symbols('vc', cls=sp.Function)
> C, R, vo = sp.symbols('C R vo')
> cos = sp.cos
>
> diffeq = vc(t).diff(t) + (1/(R*C))*vc(t) - (vo*cos(w*t)/(R*C))
> res = sp.dsolve(diffeq, ics={vc(0): 0})
>
> res
>
> The result is really confused, compared to the nice wikipedia result.
> Is there a way to get the same result? (follow the latex)
>
> {\displaystyle v_{C}(t)=v_{C}(0)e^{-{\frac {t}{\tau }}}+K\cos(\omega t+\theta
> )}
>
> Thanks so much
>
> --
> You received this message because you are subscribed to the Google Groups
> "sympy" group.
> To unsubscribe from this group and stop receiving emails from it, send an
> email to sympy+unsubscr...@googlegroups.com.
> To view this discussion on the web visit
> https://groups.google.com/d/msgid/sympy/8f5df611-4aa2-4f09-9f5b-a9d528c24333n%40googlegroups.com.
--
You received this message because you are subscribed to the Google Groups
"sympy" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to sympy+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/sympy/CAKgW%3D6JRCQdUktVUEvY3o4R-iyT7sAGBJMgmX0-aGYcPngLOTg%40mail.gmail.com.
