hi Daniyar, You mean setting Tspan as Inf ? if yes then :
i do not think its logic to make the bound as Inf , but instead make a large value , ok try this :
>>f=@(y,x) sin(x+y).*exp(-x.^2-y.^2)
>>[t,y]=ode15s(f,[1 2],[0 Inf])
>>[t,y]=ode15s(f,[1 2],[0 1e+9])
does this helps ?
