I want to see the stream line plot for y' = -sin(x).
The solution of this ODE is coded as:
x = var('x')
y = function('y')(x)
de = diff(y,x) == -sin(x)
h = desolve(de,y)
solution = h.function(x)
show(solution)
However, I am not able to get the right expression in the code for the stream line plot using the functions already defined above.
I can make it happen like this:
x,y=var('x','t')
streamline_plot(cos(x),(t,-10,10),(x,-5,5))
But it is not elegant, nor does it take advantage of the variables h
, solution
or y
already defined.
What would be the more natural expression for the general (and a particular solution) plots?