Question: Help: Differential Equation on Decomposition Method (Procedure contains error)

Hello there, I have a problem solving the differential equation based on decomposition method. There are 3 operator which is N, R, G and u. N=4, g=-77, R=79. Below is the procedure, and yet there is an error at the last end statement. Please help asap. Thanks Regards, Sheby > Adomian:=proc(ICs,Eqs,n) local nEq,t0,m,v0,R,N,g,nOrder,u,k,nu,tempN,A,i; ICs:=[0,[0,2]]: Eqs:=[[u->-77*diff(u,t),u->79*diff(u,t)^2 + 4*u*diff(u,t),t->9]]: n:=13; nEq:=nops(ICs)-1; nOrder:=[seq(nops(ICs[m]),m=2..nEq+1)]; t0:=ICs[1]; v0:=[seq(ICs[m],m=2..nEq+1)]; R:=[seq(Eqs[m,1],m=1..nEq)]; N:=[seq(Eqs[m,2],m=1..nEq)]; g:=[seq(Eqs[m,3],m=1..nEq)]; u:=array(1..nEq,0..n); for m from 1 to nEq do u[m,0] :=g[m](t); for i from 1 to nOrder[m] do u[m,0]:=v0[m,nOrder[m]-i+1]+int(subs(t=s,u[m,0]), s=t0..t); end do; end do; for k from 0 to n-1 do nu:=seq(sum('lambda^i*(u[m,i])','i'=0..k),m=1..nEq); for m from 1 to nEq do tempN:=unapply((N[m])(nu),lambda); A[m,k]:=(D@@k)(tempN)(0)/k!; end do; for m from 1 to nEq do u[m,k+1]:=-R[m](seq(u[i,k],i=1..nEq))-A[m,k]; for i from 1 to nOrder[m] do u[m,k+1]:=int(subs(t=s,u[m,k+1]),s=t0..t); end do; end do; end do; [seq(sum('u[m,k]','k'=0..n),m=1..nEq)]; end:

Please Wait...