Question: Solving for the coefficient in a polynomial with initial conditions.

I have two polynomials, say

 

x(t) = a[0] + a[1]t + a[2]t^2

y(t) = b[0] + b[2]t + b[3]t^2

 

with the following conditions, x(0) = 2 and y(0) = 1. The polynomial are related by

 

x'(t) = y(t)

y'(t) = x(t)

 

and x(0) = 2, y(0) = 1.

 

One can show that the solution is x(t) = 2 + t + t^2 + O(t^3) and y(t) = 1 + 2t + 0.5t^2 + O(t^3).

 

I am trying to write something so that it works on a larger system. I have a set of nonlinear DEs which i want to solve this way. Here is a sample code of a much simplifed problem.

 

See the uploaded file for details.

 

edit: link is reparied. I uploaded the wrong file initially

 

edit2: Some improvements were made, but Maple still not returning.

 

 

 

2-2_tests_01.mw

Please Wait...