Differences
This shows you the differences between two versions of the page.
| Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
| public:calculator:progs:dm41x_progs [14/02/26 11:10 GMT] – [Useful Memo prog] john | public:calculator:progs:dm41x_progs [14/02/26 11:19 GMT] (current) – [Using Solve and Integrate together] john | ||
|---|---|---|---|
| Line 30: | Line 30: | ||
| LBL ¬MYSOLVE | LBL ¬MYSOLVE | ||
| ¬MOD // the label of the program to be solved is put in ALPHA | ¬MOD // the label of the program to be solved is put in ALPHA | ||
| - | SOLVE // we XEQ SOLVE (the Advantage Pac solver | + | SOLVE // we XEQ SOLVE (the Advantage Pac solver) |
| RTN | RTN | ||
| </ | </ | ||
| - | === MOD === | + | === Solve : MOD === |
| The program to be solved | The program to be solved | ||
| Line 40: | Line 40: | ||
| < | < | ||
| LBL ¬MOD | LBL ¬MOD | ||
| - | ¬BESSEL | + | ¬BESSEL |
| - | STO 00 | + | STO 00 // We store our current guess for X in R 00 to be accessible by the integral finction |
| + | 0 // put upper and lower limits for the integration in Y and X | ||
| PI | PI | ||
| - | INTEG | + | INTEG // xeq INTEG from Advantage Pac |
| - | ABS | + | ABS // our current guess, make it positive |
| - | X<=Y? | + | X< |
| - | ClX | + | ClX // if our guess has returned a value from the integral that's less than the inherent accuracy then that's GOOD ENOUGH, we set X to zero and SOLVE completes |
| - | X≠0? | + | X≠0? |
| LastX | LastX | ||
| RTN | RTN | ||
| Line 54: | Line 54: | ||
| </ | </ | ||
| - | === The integral | + | === Integral : BESSEL === |
| < | < | ||
| Line 67: | Line 67: | ||
| </ | </ | ||
| + | |||
| + | |||
| + | === Running it === | ||
| + | |||
| + | * '' | ||
| + | * put initial guesses in '' | ||
| + | * < | ||
| + | * try '' | ||
| + | * Run '' | ||
| + | |||
| + | After a few minutes of Goose-flying we get the answer '' | ||
| + | |||
| + | |||