% deut: Solves ODE for toy deuteron problem.
global x0  E;
% Domain outer boundary ...
xmax = 60.0;
% Tolerance parameters ...
abstol = 1.0e-8;
reltol = 1.0e-8;
options = odeset('AbsTol', abstol * [1 1]', 'RelTol', reltol);
% Parameters ...
x0 = 2.0; 
E = -0.1
% Integrate ODE ...
[xout yout] = ode45(@fcn_deut, [0.0 xmax], [0.0 1.0]', options);
% Make plot and save as JPEG
figure(1);
clf;
hold on;
axis square;
xlabel('x');
ylabel('u');
plot(xout, yout(:, 1));
% print('deut.jpg','-djpeg');
