% y' = f(x)
f = @(x) exp(-x .^ 2) ;
g = @(x,y) f(x) ; 
x0_xE = [ 0 1 ] ; % estremi d'integrazione
y0 = 0 ;          % valore iniziale
[ x , y ] = ode45( g , x0_xE , y0 ) ;
y_verifica = F_integrale( f , 0 , x ) ;

% y' = a(x)y
a = @(x) -2*x ;
g = @(x,y) a(x) * y ;
x0_xE = [ 0 3 ] ; % estremi d'integrazione
y0 = 1 ;          % valore iniziale
[ x , y ] = ode45( g , x0_xE , y0 ) ;
y_verifica = exp(-x .^ 2) ;

% my"+cy'+ky=f(x)
f = @(x) 5000 * onda_quadra(x, 2) ;
%w = 2*pi*10 ;
%f = @(x) 10000 * sin(w*x) ;
m = 1 ;     % massa
c = 10 ;    % dissipazione
k = 10000 ; % costante elastica
G = @(x,y) ( f(x) - c*y(2) - k*y(1) ) / m ;
g = @(x,y) [ y(2) ; G(x,y) ] ; 
x0_xE = [ 0 4 ] ; % estremi d'integrazione 
y0    = [ 0 0 ] ; % valori iniziali 
[ x , y ] = ode45( g , x0_xE , y0 ) ;
plot( x , y(:,1) ) ;
%B = 10000 / abs( k - m*w^2 + 1i*w*c ) ;