Showing posts with label HP 48G. Show all posts
Showing posts with label HP 48G. Show all posts

Saturday, July 31, 2021

HP 48G: Truth Plots

HP 48G: Truth Plots





Calculators Featured:  HP 48G, HP 48G+, HP 48GX


Truth Plots


On the HP 48G, a truth plot tests an inequality or equality for each point of a Cartesian grid.   If the statement is true, the pixel is plotted.  If not, the pixel remains unchanged.  Since each pixel is being tested, a truth plot does take some time to plot.  


Example: The inequality x + y ≤ 2, plots a pixel every time inequality is true.


To access a truth plot on the HP 48G:


1.  Press [ right shift ], [ 8 ] (PLOT).  

2.  The expression is stored in EQ.  Default:  X is the independent variable and Y is the dependent variable.  

3.  Make any necessary adjustments.  If you have to reset your parameters, press [ left shift ], [ 8 ] (PLOT), (PPAR), (RESET).


Truth plots are not unique to the HP 48G series: the current HP Prime can easily and quickly execute truth plots.  The truth plots are renamed "Advanced Graphing App".


Here are examples of truth plots, with the HP 48GX is one the left, the HP Prime is on the right.  The window is set up as x = [-5, 5] and y = [-5,5 ]. 


x^2 + y^2 ≥ 5




(x^3 - y^2) mod 5 ≥ 3




sin(x*y) < 3*cos(x*y)

Radians angle mode is used


Eddie


All original content copyright, © 2011-2021.  Edward Shore.   Unauthorized use and/or unauthorized distribution for commercial purposes without express and written permission from the author is strictly prohibited.  This blog entry may be distributed for noncommercial purposes, provided that full credit is given to the author. 


Saturday, July 24, 2021

HP 48G: Multiple Equation Solver and The Equation Library

HP 48G: Multiple Equation Solver and The Equation Library 





Calculators Featured:  HP 48G, HP 48G+, HP 48GX


The Equation Solver


Any program in the HP 48G family, we can set up the calculator to solve a set of equations.  The commands involved:


STEQ:  Store a list of equations in the variable EQ.  The variable EQ is used in the solver.


MINIT:  This command initializes the solver.  The system variable MPAR is created.


MITM:  This optional command allows the solver to have a custom title and arrange the variables in a custom order.  Arguments:  2:  title string, 1:  the list of variables.  The list can have a blank string, " ".  All the variables must be included.  


MSOLVR:  Starts the solver.


In the solver:


( var ):  store

[ left shift ] ( var ):  solve

[ right shift ] ( var ):  recall


An RPL Template:


<< ... 

{ list of equations }  STEQ MINIT

"title" { list of variables } MITM

MSOLVR

... >>


HP 48G Program:  HCURVE


<<  DEG

'T=R*TAN(Δ/2)'

'L=R*Δ*π/180'

'C=2*R*SIN(Δ/2)'

3 →LIST STEQ MINIT

"Horizontal Curve"

{ Δ R T L C } MITM

MSOLVR >>


Variables:


Δ = central angle (in degrees)

R = central radius

T = tangent length

L = arc length of the horizontal curve

C = chord length


One of Δ or R must be known.  When Δ and R are both entered or solved for, L, T, and C can be solved for.


Example 1:


Inputs:

Δ: 45

R: 100


Outputs:

T = 41.4213562373

L = 78.5398163398

C = 76.536686473


Example 2:


Inputs:

L:  1706

Δ: 54


Outputs:

R = 181.012221943

T = 2.2303337494

C = 164.35658198


The Equation Library


The equation library can be accessed by pressing [ right shift ] [ 3 ].  A list of various subjects are listed.  Options:


[ F1 ] (SI):  SI units are used

[ F2 ] (ENGL):  English units are used

[ F3 ] (UNIT): turn units on and off


Once a topic is selected, another menu s presented where all the equations are contained.  Depending on the topic, one or more equations are selected.  When ready, select (SOLV).  


Example 1:  Triangle


Plane Geometry -> Triangle -> Equation 2 of 6


A=(b*h)/2

P=b+√(v^2+h^2)+√((b-v)^2+h^2)





Example:  (default SI units are used)


Inputs:

P: 100 cm

b: 40 cm

h: 35 cm


Outputs:

V = 20 cm

A = 700 cm^2


Example 2:  Longitude Waves


Waves -> Longitudinal Waves -> Equation 1 for 4


s = sm cos( k * x - ω * t )

v = λ * f

k = 2 * π / λ

ω = 2 * π * f


Example:  (default SI units are used)


Inputs:

λ: 0.2264 cm

f: 1400 Hz

t: 15 s

sm: 1.5 cm


Outputs: 

v = 314.44 m/s

k = 0.279750013677 r/cm

ω = 8976.45943005 r/s


Eddie


All original content copyright, © 2011-2021.  Edward Shore.   Unauthorized use and/or unauthorized distribution for commercial purposes without express and written permission from the author is strictly prohibited.  This blog entry may be distributed for noncommercial purposes, provided that full credit is given to the author. 


Saturday, July 17, 2021

HP 48G: Business Applications

 HP 48G:  Business Applications


Calculators Featured:  HP 48G, HP 48G+, HP 48GX

Percent Addition 

HP 48G Program: %PLUS

<< OVER SWAP % + >>

Stack:
2:  amount
1:  percent

Example 1:   150 + 17%

150 [ ENTER ] 17 (%PLUS) 
Result:  175.5

To subtract, make the percent negative.

Example 2:  170 - 18%

180 [ ENTER ] 18 [ +/- ] (%PLUS)
Result:  139.4

Interest Rate Conversion to APR

The program →APR converts a nominal interest rate (NOM) to the effective interest rate (APR/EFF) given the rate and the compound periods per year.

HP 48G Program:  →APR

<< "→APR"
{ { "I%YR:" } { "PYR:" } }
{ } { } { }
IF INFORM 
THEN OBJ→ DROP
SWAP 1 % OVER / 1 +
SWAP ^ 1 -  100 * 
"APR" →TAG END >>

Example   

nominal rate: 7.5%, periods per year: 6

Result:  
APR: 7.738318055 (%)

Source:
Thomsett, Michael C.  The Real Estate Investor's Pocket Calculator 2nd Ed.  American Management Association: New York.  2017

Present Value of a Growing Payment

The program PVGROW calculates the present value of a growing annuity.   Each payment grows by the rate g% while the discount rate is r%.

PV_grow = PMT_base / (R% - G%) * (1 - (1 +G%) / (1 + R%) ^ n )


HP 48G Program:  PVGROW

<< "PV GROWING ANNUITY"
{ { "1st Pmt:" } { "I/YR:" } { "Growth Rate:" } { "N: "} }
{ } { } { } 
IF INFORM
THEN OBJ→ DROP 4
ROLLD 1 % SWAP 1 %
SWAP - LASTARG 1 + 
SWAP 1 + / 4 ROLL 6
1 SWAP - 3 ROLLD /
* "PV" →TAG END >>

Example

1ST PMT:  1000.00
I/YR:  3.5 (%)
GROWTH RATE:  4 (%)
N: 36

Result:  
PV:  37890.784118

Source:
"Present Value of a Growing Annuity" financeformulas.net  https://financeformulas.net/Present_Value_of_Growing_Annuity.html   Retrieved December 13, 2020.  


Interest Only Loan

The program INTONLY calculates payments from:

Interest-only payment
Payment (P&I) when the principal starts being due

HP 48G Program:  INTONLY

<< "INTEREST ONLY"
{ { "LOAN:" } { "I/YR:" } { "INT YRS:" } { "MAIN YRS:" } }
{ } { } { }
IF INFORM 
THEN OBJ→ DROP 12
* 'N' STO ROT DUP
'PV' STO ROT DUP
'I%YR' STO 1 % * 12
/ NEG 'IMPT' DUP2 
STO →TAG SWAP DROP
0 'FV' STO 12 'PYR'
STO 'PMT' TVMROOT
'PMT' →TAG END >>

Example

LOAN:  185000.00
I/YR:   4.55 (%)
INT YRS:  5
MAIN YRS: 30

Results:

IMPT: -701.4583333333   ($701.46 each month for the first 5 years)
PMT: -942.871969103   ($942.87 each month for the rest of the 30 years)

Source:
Kapfidze, Tendayi  "What Is an Interest-Only Mortgage and How Does It Work?"  Edited by Deborah Kearns  LendingTree.  LendingTree, LLC, Charlotte, NC:  2021.   https://www.lendingtree.com/home/mortgage/interest-only-mortgages/  Retrieved March 9, 2021


Eddie

All original content copyright, © 2011-2021.  Edward Shore.   Unauthorized use and/or unauthorized distribution for commercial purposes without express and written permission from the author is strictly prohibited.  This blog entry may be distributed for noncommercial purposes, provided that full credit is given to the author. 

Saturday, July 10, 2021

HP 48G: Scientific Applications

HP 48G:   Scientific Applications



Calculators Featured:  HP 48G, HP 48G+, HP 48GX


Theoretic Freezing Level


The program FRZLVL calculates the theoretical freezing level in feet above sea level given the outside temperature at a given altitude.  The results are displayed in feet (the feet unit is attached).  


HP 48G Program: FRZLVL


<< "FREEZING LEVEL"

{ { "ALT:" "DEFAULT FT" }

{ "TEMP:" "DEFAULT °C" } }

{ }  { }  { }  

IF INFORM

THEN OBJ→ DROP DUP TYPE

IF 13 == THEN '1_°C' CONVERT

ELSE '1_°C' →UNIT 

END SWAP DUP TYPE 

IF 13 == THEN '1_ft' CONVERT

ELSE '1_ft' →UNIT

END SWAP '1000_ft/°C' 

* DUP2 2 / +

"DRY" →TAG 3 ROLLD

1.5 / + "WET" →TAG END >>


Instructions


1.  Run FRZLVL

2.  ALT:  Enter altitude.  You can attach units if you want, the default is altitude is in feet.

3.  TEMP:  Enter outside air temperature.  You can attach units if you want, the default unit is degrees Celsius (°C).

4.  Results: freezing dry level in feet and freezing wet level in feet.  


Example


Altitude:  1,750 ft

Temperature:  -2°C  (28.4°F)


Results:

DRY:  750_ft

WET:  416.666666667_ft


Source:

Hewlett Packard  "Predicting Freezing Levels" HP 65 Aviation Pac 1.  Hewlett Packard.  1974


Lens Calculations 





Variables


Please note the inequality restrictions.  


R = radius of curvature

D = diameter of the lens, D ≤ 2*R

S = sag S ≤ R

θ =  tangent angle between lens axis and reflection 


HP 48G Program: LENSANG


<< 'R=(D^2+4*S^2)/(8*S)' STEQ

{ { "R/D" << 'D' STO 'R' STO EQ  'S' 0 ROOT 'S' →TAG >> } 

{ "R/S" << 'S' STO 'R' STO EQ 'D' 0 ROOT 'D' →TAG >> }

{ "D/S" << 'S' STO 'D' STO EQ 'R' 0 ROOT 'R' →TAG >> } 

{ "→θ" << 'D' RCL 'R' RCL 'S' RCL - 2 * / ATAN 'θ' DUP2 STO →TAG >> } }

TMENU >> 


Instructions


1.  Run LENSANG.  A temporary custom menu appears.  

2.  Enter two known amounts and press the appropriate key to store the variables:


If R and D are known:  R [ ENTER ] D [ F1 ] 

If R and S are known:  R [ ENTER ] S [ F2 ]

If D and S are known:  D [ ENTER ] S [ F3 ]


The third variable between R, D, and S is calculated and displayed as a result.


3.  Press [ F4 ] to calculate the angle, θ.  


Results are stored in the variables R, S, D, and θ.


Examples


Known:  R = 5.8, D = 7.6

LENSANG  5.8 [ ENTER ] 7.6  [ F1 ] (R/D)

Result:  S: 1.41821953996

[ F4 ] (θ).  Result:  θ:  40.9327245742


Known:  D = 11, S = 9

LENSANG 11 [ ENTER ] 9 [ F3 ] (D/S)

Result:  R: 6.18055555556

[ F4 ] (θ).  Result:  θ:  -62.8591312297


Source:


Tuchscherer, L.D.  "Lens Calculations - SAG, ANGLE, MIN/MAX" HP 67 Optics  Hewlett Packard.  1978.


Period of a Pendulum


The program PENDULUM calculates the period of a pendulum given two parameters:


ANGLE: the angle that the pendulum swings out.  The angle must be entered in radians.


LENGTH:  The length of the pendulum in meters. 


Unlike most calculations, the program uses an Legendre elliptic integral of the first kind to increase accuracy.  


T = 4 * √(L / g) * ∫( dx / √(1 - k^2 * sin^ x), x, 0, π/2) where k = sin(α/2)


g = Earth's gravitational constant = 9.80665 m/s^2


This program does not use units.  


HP 48G Program: PENDULUM


<< RAD 

"PENDULUM"

{ { "ANGLE:" "IN RADIANS" }

{ "LENGTH:" "IN METRES" } }

{ } { } { } 

IF INFORM

THEN OBJ→ DROP

9.80665 / √ 4 *

SWAP 2 / SIN SQ 0 θ π 2 /

→NUM 3 ROLL 'X' 

SIN SQ * 1 SWAP - √

INV 'X' ∫ →NUM *  

"PERIOD" →TAG END >>


Example


ANGLE:  π/6   (enter this as 'π/6')

LENGTH: 1.5  (m)

Result:  

PERIOD:  2.50011881685 (s)


Source:

Steers, Hugh.  "The Pit & Pendulum"  Datafile.  ISSN 1352-8254.  Handheld and Portable Computer Club.   V29 N1.  January - March 2010


Eddie


All original content copyright, © 2011-2021.  Edward Shore.   Unauthorized use and/or unauthorized distribution for commercial purposes without express and written permission from the author is strictly prohibited.  This blog entry may be distributed for noncommercial purposes, provided that full credit is given to the author. 


Sunday, June 27, 2021

Review: Victor 6500 Loan Calculator, HP 48G Retro Month

Review:  Victor 6500 Loan Calculator, HP 48G Retro Month


Victor 6500 Review






Quick Facts

Company:  Victor
Year Introduced:  about 2010, still in production
Battery:  Solar and 1 AA battery
Digits: 12
Memory Registers:  2: Independent (M) and Grand Total (GM)
Cost: About $20.00 to $25.00 US


It is rare that a desktop calculator features advanced financial or scientific functions, and that is what attracted me to the Victor 6500.  The keys respond quickly which allows for fast typing.   The display has room for 12 large digits which makes reading the display easy.

Features

*  Fix Display Switch: Float, 0, 2, 4, Add Mode
*  Round Setting Switch:  cut (always round results down), 5/4 (regular rounding rules are followed), up (always round up)
*  Time Value of Money - will be discussed next
*  %, Δ%, √
*  Tax Rate:  setting, TAX+, TAX-
*  Double Zero Key

Memory Registers

There are two memory registers:

Independent Memory (M).  This the traditional memory register that comes with most desktop calculators.  Included are M+ and M-.  Recall and clear are combined in one key RCM.

Grand Total (GT).  The grand total memory is activated every time the equals key is pressed.   

Time Value of Money

The Victor 6500 have four finance keys:

Loan = PV
Pmt = Payment
Years = number of years = n * PYR
Int = annual interest rate

The cash flow convention (inflow/outflow) is not followed, all results are positive.  The number of payments per year is set through the P/Y setting.  

What is missing is the FV (future value) key, which leaves out some of the finance problems involving calculating balance or balloon payment.   However, we can still do a lot of calculations on the Victor 6500 fairly easier.  Here are several example calculations.  In these examples, the Victor 6500 is set to 2 Decimal Places with 5/4 rounding mode.

Example 1:  Calculating Payment of a Loan
Loan:  345,000.00
Years: 30
Interest Rate:  3.3%

[ ON/C ] [ ON/C ]
345000 [ Loan ]
30 [ Years ] 
3.3 [ Int ]
[ Pmt ]  returns 1,510.95

Example 2:  Calculating Payment of a Loan with a Down Payment
Loan:  260,000.00 with 15% down
Years: 27
Interest Rate:  3.2%

[ ON/C ] [ ON/C ]
260000 [ - ] 15 [ % ] [ Loan ]
27 [ Years ] 
3.2 [ Int ]
[ Pmt ]  returns 1,019.53

Example 3:  PITI (Payment with Principal Interest Taxes Insurance)
Loan: 300,000.00
Years: 30
Interest Rate: 3.25%
Annual Property Tax: 2,121.00
Annual Property Insurance: 541.00

[ ON/C ] [ ON/C ] [ MRC ] [ MRC ]
300000 [ Loan ]
30 [ Years ]
3.25 [ Int ] 
[ Pmt ] [ M+ ]

2121 [ + ] 541 [ = ] [ ÷ ] 12 [ = ] [ M+ ]
[ MRC ] (once)  
PITI:  1527.45

Verdict

I like the Victor 6500 since it offers a time value of money solver on a desktop calculator form.  I wish Victor would have included the future value variable.   The keys are fast and responsive.  The keyboard is pretty light but also feel stable.  

I am working on posting additional examples for the Victor 6500 in the near future.  

HP 48G Retro Month - July 2021


For the month of July 2021, every Saturday from July 10 to July 31 will have a post on the legendary HP 48G/48GX/48G+ calculator.  


Next post will be on July 10, 2021 (taking the Fourth of July weekend off).   Have a great day and thank you!  

Eddie 

All original content copyright, © 2011-2021.  Edward Shore.   Unauthorized use and/or unauthorized distribution for commercial purposes without express and written permission from the author is strictly prohibited.  This blog entry may be distributed for noncommercial purposes, provided that full credit is given to the author. 

Sunday, September 23, 2018

HP 48G: Stefanescu Horatiu: Intersection of the Diagonals of a Quadrilateral

HP 48G:  Stefanescu Horatiu:  Intersection of the Diagonals of a Quadrilateral

The following program, written by Stefanescu Horatiu, calculates the intersection point of the two diagonals of a quadrilateral (partulator in Romanian). 



The program takes place in a directory PATRULATERDIAGONALE.  In order to type in the program, you will need to create a directory by using the CRDIR command. 

This program also works on the HP 49G series (48gii, 49G, 49g+) and the HP 50g. 

HP 48G Program PATRULATERDIAGONALE - Stefanescu Horatiu

Note:  You are going to notice how all the subroutines flow together, as the first program calls the next program. You will need to enter and store all the programs for the entire routine to work correctly.

On a personal note, I like how the text is spelled out one letter at time.  Hoartiu accomplishes this by using a FOR-NEXT structure which includes the commands 1 I  SUB 1 DISP.  (with I as the counter variable).

The text of the program is in Romanian. 

The program is presented here with permission from the author, with full credit going  to Stefanescu Horatiu.  I thank Horatiu in allowing me to feature his work on my blog.

STARTDIAGONALE:

<< CLEAR CLLCD 1 110 FOR I
"ACEST PROGRAM CALCUL
EAZA COORDONATELE PUN
CTULUI DE INTERSECTIE
A DIAGONALELOR UNUI
PATRULATER OARECARE."
1 I SUB 1 DISP
NEXT 0 WAIT
INTRODUCERE >>

INTRODUCERE:  (User input screen)

<< CLLCD 
"COORDONATE PATRULATER"
{"Ax" "Ay" "Bx"
"By" "Cx" "Cy" "Dx"
"Dy" } { 2 4 } { }
{ } INFORM DROP
OBJ→ DROP 'DY' STO
'DX' STO 'CY' STO
'CX' STO 'BY' STO
'BX' STO 'AY' STO
'AX' STO CLLCD
EMURILE >>

EMURILE:  (the start of the calculation routines)

<< 'CY-AY' →NUM 'M1' STO 
'-AX*(CY-AY)' →NUM 'M2' STO 
'CX-AX' →NUM 'M3' STO
'-AY*(CX-AX)' →NUM 'M4' STO
ENURILE >>

ENURILE:

<< 'BY-DY' →NUM 'N1' STO
'-DX*(BY-DY)' →NUM 'N2' STO
'BX-DX' →NUM 'N3' STO
'-DY*(BX-DX)' →NUM 'N4' STO
ESURILE >>

ESURILE:

<< 'N1-N3*M1/M3' →NUM 'S1' STO
'N3*(M2-M4)/M3+N4-N2' →NUM 'S2' STO
PURILE >>

PURILE:

<< 'S2/S1' →NUM 'PX' STO
'(PX*M1+M2-M4)/M3' →NUM 'PY' STO
REZULTATE >>

REZULTATE:   (display the results)

<< CLLCD 1 59 
FOR I 
"Coordonate punct de
intersectie ale dia
gonalelor: P(x,y)"
1 I SUB 1 DISP NEXT 4 FIX
"P(" PX + →STR ";" +
→STR PY + →STR ")" 
+ →STR 5 DISP 0 WAIT 
STERG >>

STERG:  (cleanup, purge working variables)

<< { AX AY BX BY CX 
CY DX DY PX PY M1
M2 M3 M4 N1 N2 N3
N4 S1 S2 } PURGE
CLLCD CLEAR >>

Instructions:

Set your calculator the folder PATRULATERDIAGONALE.  Run STARTDIAGONALE and enter the coordinates.  The rest of the program runs automatically.  Start with the left hand corner and go clockwise. 

Note:  I typed Horatiu's program on an HP 49G.    In REZULTATE, for some reason I couldn't get the FIX 4 to take effect, even when I switched the calculator to Approximate mode.  So in order to get both PX and PY on the screen, I added a carriage return character (↲, [right shift] [ 2 ]) after the semicolon.  This creates a new line allowing both coordinates to be displayed.

Example:

Example 1:   (Horatiu's example)
A:  (ax, ay):  (2, 10)
B:  (bx, by):  (11, 12)
C:  (cx, cy):  (13, 6)
D:  (dx, dy):  (1, 4)

P: (px, py):  (6.46875, 8.375)
FIX 4:  (6.4688, 8.3750)

Example 2:  (my own example)
A:  (ax, ay):  (0, 8)
B:  (bx, by):  (5, 12)
C:  (cx, cy):  (7, 7)
D:  (dx, dy):  (2, 1)

P:  (px, py): (3.72649999999, 7.4625)
FIX 4:  (3.7265, 7.4625)

Source:

Horatiu, Stefanescu "Calculul Coordonatelor Punctului de Intersectie A Diagonalelor Unui Partulater Oarecare"  2018. 

Eddie

All original content copyright, © 2011-2018.  Edward Shore.   Unauthorized use and/or unauthorized distribution for commercial purposes without express and written permission from the author is strictly prohibited.  This blog entry may be distributed for noncommercial purposes, provided that full credit is given to the author.  Please contact the author if you have questions.

Sunday, August 7, 2016

HP 48G: When is Easter?

HP 48G:  When is Easter?

The following program will accurately calculate the date when Easter occurs.  Special thanks to Peter Ketel for furnishing this program to me and he has given me permission to share it on this blog.  Much appreciated, Peter!

This should be able to be adopted on the HP 49G and HP 50G easily.

HP 48G Program EASTER (591 bytes)
≪ → y
≪ y 19 MOD 1 + y
100 / FLOOR 1 + DUP
OVER 3 * 4 / FLOOR
12 – SWAP 8 * 5 +
25 / FLOOR 5 - → a b c d
≪ y 5 * 4 /
FLOOR 10 – c – 11 a *
20 + d + c – 30
MOD 30 + 30 MOD → e f
≪ IF f 24 == f 25 == a 11 >
AND OR
THEN f 1 + ‘f’ STO END
44 f – DUP
IF 21 < THEN 30 + END → g
≪ g 7 + e g
+ 7 MOD – DUP
IF 31 > THEN 31 - .0001 * .04
ELSE .0001 * .03 END
+ y +
≫ ≫ ≫ ≫ ≫

Input:  Enter the year and run the program

Tests and Results:

1976 returns 1976.0418  (April 18, 1976)
1985 returns 1985.0407  (April 7, 1985)
2017 returns 2017.0416  (April 16, 2017)
2050 returns 2050.041  (April 10, 2050)

This blog is property of Edward Shore, 2016.


Sunday, February 21, 2016

Blackjack: Approximate Chance of Busting (HP 48G)

Blackjack:  Approximate Chance of Busting

In standard game of Blackjack, the main objective is beat the dealer’s hand.  Usually this occurs if your hand has a better total (no more than 21) than the dealer’s hand.  At any time if the total of your hand exceeds 21, you bust, or automatically lose. 

As your score increases, the number of cards that can cause you to bust increase.

If your score is 11 or less, you cannot bust on the next card.  Even if you get an Ace, Aces are 1 or 11, and are always used to your best benefit.

If your score is 12, if the next card is any face card or ten, you bust.  There are 16 face cards and tens in a standard deck of 52 playing cards.  The face cards are Jacks, Queens, and Kings.  In Blackjack, face cards count as 10 points.  Below is a table of score and bust cards.

Your Score
Bust Cards
Number of Bust Cards per Deck of 52 Playing Cards
12
Face Cards, Tens
16
13
Face Cards, Tens, Nines
20
14
Face Cards, Tens, Nines, Eights
24
15
Face Cards, Tens, Nines, Eights, Sevens
28
16
Face Cards, Tens, Nines, Eights, Sevens, Sixes
32

If your score is 17 and above, basic Blackjack strategies dictate that you stand (draw no more cards).

Approximating the Odds of Busting

To calculate the odds of busting, we have to determine two things:  (I) the number of bust cards that remain, and (II) the number of cards remaining to be dealt.

The number of bust cards that remain are:

d * (4*t – 32) – b  

Where:

d = the number of decks in play.  Decks are assumed to standard decks of 52 playing cards, without jokers.
t = your score.  For our purposes, 12 ≤ t ≤ 16.  Outside this range, this formula does not make sense.
b = the number of bust cards used.  Unless you are counting cards, this is an estimate number. 

The number of cards remaining to be dealt:

52*d – n

Where:
d = the number of decks in play (same as above)
n = number of cards that have been played

The odds of busting on the next card is:
OB = (d * (4*t – 32) – b )/( 52*d – n)

Let’s calculate some approximate odds of busting. 

Scenario #1:  The table is just you and the dealer.  It is the first hand the session and 6 decks are used.  On the first hand both you and the dealer have been dealt a face card. 

Variables: d = 6, b = 2 (two bust cards used), n = 4 cards dealt (one card has an unknown value)

Score (t)
Chance of Busting
12
30.5195 (%)
13
38.3117 (%)
14
46.1039 (%)
15
53.8961 (%)
16
61.6883 (%)

Scenario #2:  Now the table are six players and the dealer.  Five hands have been dealt and now we are on the sixth hand.  There 6 decks being used and let’s assume that 88 cards have been used, 28 bust card used. 

Variables:  d = 6, b = 28, n = 88

Score (t)
Chance of Busting
12
30.3571 (%)
13
41.0714 (%)
14
51.7857 (%)
15
62.5000 (%)
16
73.2143 (%)

Scenario #3:  Now the table are six players and the dealer.  Same as above, let’s time assume that 40 bust cards have been used, a lot of high cards have been dealt.

Variables:  d = 6, b = 40, n = 88

Score (t)
Chance of Busting
12
25.0000 (%)
13
35.7143 (%)
14
46.4286 (%)
15
57.1429 (%)
16
67.8571 (%)

A quick observation is that the more bust cards that have been used, the lower the chances of busting later on.  Obviously, having a score of 15 or 16 is not a desirable score since the dealer stands on 17 and above.  15 and 16 are desirable if only a dealer has a weak up card (3-6) and the dealer’s chances of busting is high (assuming the hidden card is a 9, 10, or face card).
The approximations have been calculated using an HP 48GX:

HP 48G Program OBUST:

Stack:  t, d, b, n

≪ → T D B N
≪ 4 T * 32 – D * B – 52 D * N - / →NUM 100 * ≫


And no, I don’t count cards during a blackjack game.  In a real life game, I try to apply basic strategy and hope for the best. 

Until next time,

Eddie


This blog is property of Edward Shore. 2016



RPN HP 12C: Fibonacci and Lucas Sequences

  RPN HP 12C: Fibonacci and Lucas Sequences Golden Ratio, Formulas, and Sequences Let φ be the Golden Ratio: φ = (1 + √5) ÷ 2...