Demuestre cómo resolver la solución numérica de la ecuación diferencial estocástica dx=senxdt+xdw en MATLAB, donde w es el proceso de Wiener.

dxdt=@(t,x)sin(x(1))+x(1)* randn;

x0 = 1;%El valor inicial aquí es 1 y se puede personalizar.

tspan =[0 10]; % intervalo de solución

[t, val] = ode45(dxdt, tspan, x0); val)

Cuadrícula abierta

set(gca, 'xtick', 0:1:10, 'ytick', 1:.2:3.5)