% OCW Universidad de Cantabria
% Cálculo Simbólico y Numérico en Ecuaciones Diferenciales
% © Cursos 2001-2014,  M. Eugenia Pérez, UC
%
function y =soluschroedinger(a0,a1,N)
a=zeros(N,1);
a(1)=a1;
a(2)=-a0*3/2;
a(3)=-a(1)/2;
a(4)=(a0-3*a(2))/12;
for i=3:N-2
    a(i+2)=(a(i-2)-3*a(i))/((i+2)*(i+1));
end
syms x
h=a0;
for i=1:N
    h=h+a(i)*x^i;
end
nterminosdesarrollo=h
ezplot(h,0,4)


