I have a 2-D differential equation which I solve with ParametricNDSolve
with initial time given as a parameters. I solve the equation using
twoeqnpara = ParametricNDSolve[{x'[t]==-3x[t] - y[t],y'[t] == x[t],x[0]==a,
y[0]==b},{x,y},{t,0,100},{a,b}].
When I try to plot with
pp = ParametricPlot[Evaluate@Table[{x[a,b][t],y[a,b][t]}/.twoeqnpara,
{a,-2,2,0.5},{b,-2,2,0.5}],{t,0,40},PlotRange->All]
I don't get any plot. But when I do for single parameter
pp = ParametricPlot[Evaluate@Table[{x[a,a][t],y[a,a][t]}/.twoeqnpara,{a,-2,2,0.4}],
{t,0,40},PlotRange->All,PlotLegends->Range[-2,2,0.4]]
It works and give a plot.
What I think the issue is when I use the command
Table[ x+y, {x,-2,2,0.5},{y,-2,2,0.5}]
, what the command is taking a single value of x and all the values of y. What I need is all the pairs possible from the 2 arrays so that can be used as a initial condition for plotting.