by **Cuchulainn** » Mon Aug 20, 2012 6:57 pm

Hi JL,

1.- in the first test case, you give as initial condition the value in t=1: y(1)=1. I find problems when dealing with it, since I can't apply the algorigthms... how am I supposed to deal with it?

DD In principles it's just a different start point. Can you take t' = t - 1.

2.- While with test case number 2 I find that the exact solution is prety similar to those offered by the numerical methods applied (i.e predictor-corrector and runge-kutta orders 2 and 4), in the test case 1, when I choose a large T (like T=3) the exact solution goes way quickly...

DD The exact solution goes to the value 2 when T gets big.

How many mesh points are you taking?

Adaptive meshing would use eq. (6.43) to find the correct meshe size.