Table of Contents

←Programs

DM41X Progs

DM41x RAW

HP Maths Benchmark bm.raw
HP Trig Benchmark bmt.raw
8-Queens Benchmark nqu.raw
Savage Benchmark svg.raw
Memo - as below - uses ascii file name “SITES” instead of “MEMO” memo.raw
TR - from HP41cx User guide vol. 2 tr.raw
Sigma - from HP41cx User guide vol.2 sigma.raw

Using Solve and Integrate together

  1. Create a MYSOLVE wrapper that handles putting the correct SOLVE program name in ALPHA etc
  2. Create the program that is to be solved and which has as it's function an integral
  3. Create the program that defines the integral

MYSOLVE

LBL ¬MYSOLVE
¬MOD                    // the label of the program to be solved is put in ALPHA
SOLVE                   // we XEQ SOLVE (the Advantage Pac solver)
RTN

Solve : MOD

The program to be solved 

LBL ¬MOD
¬BESSEL           // We put the name of the program with the integral in ALPHA
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
INTEG             // xeq INTEG from Advantage Pac
ABS               // our current guess, make it positive
X<=Y?             // Y has the Integration's estimate of accuracy - we use this to limit our search
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?              // if we still don't have zero we put back the last guess and SOLVE continues
LastX
RTN

Integral : BESSEL

LBL ¬BESSEL
SIN
RCL 00
*
COS
PI
/
RTN

Running it

After a few minutes of Goose-flying we get the answer 2.4048

Useful Memo prog

MEMO

By Gary Goodman from HP document collection. Dated 24/3/85

65c22 Timer values

for VIA Tone generator

Needs Ladybug module

The required count is calculated from Φ2 and the audio frequency required.

Count (decimal)

{ Φ2 / ( 2 x F ) } - 2

for 1000Hz and a 1MHz Φ2 :

Count = { 1,000,000 / ( 2 X 1000 ) } - 2 = 498

498d → $01F2

Run as XEQALPHATIALPHA

Prompt F for required audio tone frequency.

Enter freq

R/S

Result is the 2-byte value for T1CL & T1CH

e.g.

XEQALPHATIALPHA 

1000

R/S

Result: 

$001F2

The Counter should be set to 

T1CH = $01
T1CL = $F2

Entered as follows (make sure you're in DECS mode in Ladybug while entering the program - otherwise the 1000000 won't be interpreted correctly). When viewing/editing the program later it will look different! 

LBL ¬TI
DECS
WSIZE 20
¬ F
PROMPT
STI 01
1000000
ENTERI
LDI 01
2
x
/
2
-
HEXS
RTN

BFil

Butterworth HPF/LPF Filter designer

Examples

5MHz 5-pole LPF

XEQALPHABFILALPHA 

MHZ ?

5 R/S 

N

Order/poles = N

Positive for HPF

Negative for LPF

5 pole LPF = -5

5CHSR/S 

COEFFICIENTS
0.618

(at this point you can RCL 01 etc. up to RCL 05 to see the calculated Butterworth Coefficients)

R/S 

LPF
393.5 -12

The calculated component values are now in Reg 01 - Reg 05

1MHz 7 pole HPF

XEQALPHABFILALPHA 

MHZ ?

1 R/S 

N

Order/poles = N

Positive for HPF

Negative for LPF

7 pole HPF = 7

7R/S 

COEFFICIENTS
2.247

(at this point you can RCL 01 etc. up to RCL 07 to see the calculated Butterworth Coefficients)

R/S 

HPF
17.88 -06

The calculated component values are now in Reg 01 - Reg 07

Page Info

Page created Thu May 26 18:35:00 2022 by John Pumford-Green

Page last updated: 14/02/26 11:19 GMT

— John Pumford-Green 24/04/23 16:07

Further Information

,