ex37ach2.m

these codes are for solving OED with matlab

% Derives series solution for ex37ch2
echo on
syms x a b c n y res
y = 1 + a*x^2 + b*x^4 + c*x^6
res = x*diff(y,2) + 2*diff(y,1) + x*y^n;
res = taylor(res,8)
solve(6*a+1)
y = subs(y,a,ans)
res = subs(res,a,ans)
solve(20*b-1/6*n,b)
y = subs(y,b,ans)
res = subs(res,b,ans)
solve(42*c+1/120*n^2+1/72*n*(n-1),c)
y = subs(y,c,ans)
yp = diff(y)
echo off