% OCW Universidad de Cantabria
% Cálculo Simbólico y Numérico en Ecuaciones Diferenciales
% © Cursos 2001-2014,  M. Eugenia Pérez, UC
%
syms  x;
f=-(( (abs(x+1)/(x+1)) +1)/2) *(x^2-1)+(( (abs(x-1)/(x-1)) +1)/2)*(x^2-1)
% 
ezplot(f,[-5,5])
axis([-5 5 0 2])% 
pause
%  
w2=subs(f,x,x-2*0.25);
w1=subs(f,x,x+2*0.25);
w=(w1+w2)/2;
ezplot(w,[-5,5])
axis([-5 5 0 2])
pause
%  
ww1=subs(f,x,x+2*1/2);
ww2=subs(f,x,x-2*1/2);
ww=(ww1+ww2)/2;
ezplot(ww,[-5,5])
axis([-5 5 0 2])
pause
%  
www1=subs(f,x,x+2*2);
www2=subs(f,x,x-2*2);
www=(www1+www2)/2;
ezplot(www,[-6,6])
axis([-6 6 0 2])
pause
%  
syms x t
solucion=(subs(f,x,x+2*t)+subs(f,x,x-2*t))/2
ezsurf(solucion,[0,2],[-6,6])
  