Question: Need advice for manipulating large amount of Data

HI all,

 

Here is an excerpt of a procedure I am working on:

>for i from 0 to N while T[i] <= rb do

>......

>Calculations

>.....


>X[1+i] := evalf(F(T[i],X[i]),15);
>T[1+i] := evalf(T[i]+h, 15);
>end do:

where rb is the rignt boundery

For example, if N equal 1000, knowing that X(T[0])=X[0] and F is an arrow procedure that depends on T[i] and X[i],  I will have 1001 elements that I could put in a list like:

[  [t[0],X[0]], [t[1],X[1]], ....]

The objective is to plots[pointplot](list).  But I would like to know if there is a way to plot the numeric calculations AND having the possibilities on finding the value of X for a particular T that is not necessarily an exact match of a T[i] in the list because of the increment h.

I have try to make a procedure (A) in a procedure (B) so I could use A to calculate a particular X(T) and B to plot the list but without any succes.  Moreover, the calculation take times because of F that implies a lot of calculation to be done for each T[i] and X[i].

Any advice?

Mario

Please Wait...