clear all;
T = 1; % the duration of time
h = 0.01; % the step size
N = floor(T/h); % the number of total steps
x = zeros(length(N+1));
x(1) = 5; % assign the initial data

for j=1:1:N
    x(j+1) = x(j) + h * sin(x(j)^2);    
end

t = 0:h:N*h;
plot(t,x);
    