TweetFollow Us on Twitter

Matrices
Volume Number:4
Issue Number:9
Column Tag:Basic School

Matrices in True Basic

By Dave Kelly, MacTutor Editorial Board

Using Cramer’s Rule and True Basic to Solve Simultaneous Equations

It is only a matter of time until each Engineering or Math student comes across the need to solve simultaneous equations. Fortunately, we live in a time when someone has already figured out the solutions and we can readily use the methods which have been discovered. In this column I hope to show how True Basic’s MAT statements can be used with Cramer’s Rule to solve simultaneous equations.

In most areas of engineering, the time comes from time to time to solve equations. Many of these solutions can be found using a pencil and paper and a little bit of logical thinking. We learn from algebra that to solve a set of equations, there needs to be one equation generated for each unknown variable. It is the responsibility of the engineer to use his expertise in his field of study to determine what the equations are. The solution of the equations are sometimes trivial, but may take time to manipulate the equations with simple algebra. Eventually the answer is found.

Computers are supposed to make life easier. Right? Well, of course! So let’s see how we can solve some equations using True Basic. First, we need a problem to solve: We want to find the current in each of three paths of a circuit. The circuit is shown in figure 1.

To determine the equations which we will use to solve the problem we will use Ohm’s law and Kirchhoff’s voltage law. For the non-engineer, Ohm’s law is a mathematical way to express the relationship between voltage, current, and resistance:

E = I * R

Figure 1.

where R is resistance, E is the voltage across the resistance, I is the current through the resistance. This simple law is basic to a beginning electronics class. Ohm’s law is named after a German physicist, George Simon Ohm, who published a pamphlet in 1827 which contained the results of his efforts to describe and relate currents and voltages mathematically. It has since been shown, however, that this result was discovered 46 years earlier in England by Henry Cavendish.

Kirchhoff’s voltage law is named after Gustav Robert Kirchhoff, a German university professor who was born about the time Ohm was doing his experimental work. Kirchhoff’s law states that the algebraic sum of the voltages around any closed path in a circuit is zero. I won’t attempt to prove this as there are many electronics books which explain how Kirchhoff’s voltage law works.

Figure 2.

Our first equation is derived by taking the loop shown in Figure 2. The path of the loop does is not important as long as the path is closed. The equation is:

-7 + 1(i1 - i2) + 3 (i3 - i2) + 1i3 = 0

which equates to:

i1 - 4i2 + 4i3 = 7.

Our second equation is derived by taking the loop where current i2 is as shown in Figure 1. The equation is:

1(i2 - i1 ) + 2 i2  + 3(i2 - i3) = 0

which equates to:

-i1 - 6i2 - 3i3 = 0.

The last equation comes from the relationship of the current source and the currents i1 and i2:

i1 - i3  =  7.

Now we have three equations and three unknown values to solve for. An easy method to use is given by Cramer’s Rule. Gabriel Cramer (1704-1752), was a Swiss mathematician, also known by his book on the theory of curves, which appeared in 1750 in Geneva. Cramer’s Rule states that if we have a system of linear equations

a11x1 + a12x2 + + a1nxn = b1,
a21x1 + a22x2 + + a2nxn = b2,
            ,
an1x1 + an2x2 + + annxn = bn

with n equations and n unknowns., the n X n coefficient matrix of the system is denoted by A. The determinant A 0, the components of the solution of the system are given by the formula,

where the matrix Bk is the same as A except that the elements aik, 1 ¾ i ¾ n, in the kth column of A have been replace by the terms bi, 1 ¾ i ¾ n, respectively.

Cramer’s Rule says that n simultaneous equations with n unknowns can be determined by taking determinants of the matrix. Each of the unknowns is calculated as follows:

and

and

and

where

The great thing about this is that with True Basic the solution is easy to calculate because of the MAT statements included in True Basic. MAT statements include MAT INPUT (for inputing an array), MAT LINE INPUT, MAT PLOT (plotting to a True Basic chart), MAT PRINT (to print a matrix, MAT READ (read an array from data statements), MAT REDIM (to reset dimensions of a matrix), and MAT WRITE (to write an array to a file). The MAT statement is an enhancement to Basic that is very welcome in the scientific/engineering community. It has been included for many years in Hewlett Packard Basic (now version 5.0).

Using Cramer’s Rule by hand the evaluation of the determinant for i3 is:

evaluation of the other determinants gives i1= 9 Amps and i2= 2.5 Amps. The evaluation of these determinants are fairly easy by hand, but with more equations and more unknowns, the it gets to be a lot of number crunching. The program at the end will do the number crunching for an n X n matrix, limited by the amount of memory you have. The number of equations is calculated in a function which uses the factorial function included in the True Basic Libraries to figure how much data is available in the DATA statements. The coefficients are entered into the DATA statements with equation coefficients in the first DATA statement and the solution coefficients in the second DATA statement. For solving for a large number of equations/unknowns you may break up the DATA into multiple lines.

Once the number of equations is determined, the program calls the sim_solve subroutine to solve the equation. Since the MAT statement does all the work, the subroutine is short and after a couple of MAT assignment statements, the matrix is ready to be crunched. True Basic also supplies us with the det function which will automatically calculate the determinant of an n X n matrix. It sure beats keeping track of a bunch of numbers by hand like I did when I went to school. Calculations are now done and our circuit has been analyzed! Now if they would only put a Mac on my desk I could use it!

{1}
! Simultaneous Equation Solver
! Dave Kelly
! ©1988 MacTutor

LIBRARY “Fnmlib”
DECLARE DEF Factrl

DEF Get_Equation_Count            !Read the number of equations
    WHEN ERROR IN
         LET Count=0
         DO
            READ value
            LET Count=Count+1
         LOOP
    USE
         RESTORE
    END WHEN
    LET n=0
    DO
       LET n=n+1
    LOOP UNTIL Factrl(n)>Count
    LET Count=n-1
    LET Get_Equation_Count=Count
END DEF

! Start Main Program
LET number_of_equations=Get_Equation_Count
DIM a(2,2),b(2,2),c(2,1),x(2)
! Set up all the equations
MAT READ a(number_of_equations,number_of_equations)
MAT READ c(number_of_equations,1)
MAT REDIM x(number_of_equations)
CALL sim_solve(a,b,c,x,number_of_equations)
MAT PRINT x                       ! Print the solutions
DATA 1,-4,4,-1,6,-3,1,0,-1        ! equation coefficients
DATA 7,0,7                        ! equation data
END

SUB sim_solve(a(,),b(,),c(,),x(),number_of_equations)
    IF det(a)=0 THEN EXIT SUB
    FOR j=1 to number_of_equations
        MAT b=a
        FOR i= 1 to number_of_equations
            LET b(i,j)=c(i,1)
        NEXT i
        LET x(j)=det(b)/det(a)
    NEXT j
END SUB

Someone told Dave Kelly that he was a programmer. They were kind, but not very accurate. His article “Hierarchical Menus & Colors Notes” [VOL. 4 NO. 6] contained more errors than accuracies. He claims that ZBasic source code must be compiled and the resources attached to make use of them. Why not use the function: “RefNum% = FN OPENRESFILE(FileName$)”?

This makes resources available during the programming phase and requires the deletion of only one command before creating the application.

The next glaring error was his slap-dash program that should have shown how to use submenus from ZBasic.[Excuse me, but I showed A way to to submenus from ZBasic, not THE only way. Actually, I like your way better now that you have shown it to us. Thanks -DK] As the following program demonstrates, submenus do not require special event handling loops, nor is it necessary to append special resources.

{2}
‘**********************************************
‘                    SubMenus from ZBasic
‘**********************************************
WINDOW OFF
COORDINATE WINDOW
WINDOW 1
TEXT 0,12,0,0
PRINT@(1,1)”MENU#”
PRINT@(1,3)”Item#”
‘----------------------Set The Menu------------------
MENU 1,0,1,”File”:MENU 1,1,1,”Quit/Q”
‘
SubMenu=150 ‘This is the number we’ll use for our SubMenu
MENU SubMenu,0,1,”I’m outa here”  ‘We’ll delete this from the menu bar 
as soon as we’ve attached it to a menu item.
MENU SubMenu,1,1,”Plain/P;Bold/B<B;Italic/I<I;Outline/O<O;Shadow/S<S”
SubMenuHndl&=FN GETMHANDLE(SubMenu)  ‘Handle to the SubMenu
‘
‘now delete the name from the menu bar
‘CALL DELETEMENU(SubMenu)  ‘(the menu itself still remains available)
‘
MENU 2,0,1,”Format”‘We’ll add the submenu here
‘
‘Create a menu item with the command key equivalent of
‘CHR$(27) and mark the item with the number of the
‘SubMenu instead of a “2” for a check mark
‘
MenuName$=”Style/”+CHR$(27)
MENU 2,1,SubMenu,MenuName$
‘
‘Use a negative insertion number(-1) and the handle to
‘our SubMenu - the “CALL INSERTMENU(SubMenuHandle,-1)”
‘should take place immediately after the root menu is created
‘
CALL INSERTMENU(SubMenuHndl&,-1)
MENU 2,2,1,”Itz Eazy With Z”  ‘Add more items and menus if you like
CALL DRAWMENUBAR ‘Redraw the bar to exclude the hidden menu
‘------------------------Events--------------------------------
ON MENU GOSUB “Handle Menu”
MENU ON ‘We’ll just track menu events
“Loop”  ‘Loop and wait
GOTO “Loop” ‘Getting dizzy?
MENU OFF
“Handle Menu”’----------------Menu Handling------------------
MenuID=MENU(0)
ItemID=MENU(1)
MENU  ‘Get results of the menu action
IF MenuID=1 AND ItemID=1 THEN END  ‘End if user selected “Quit”
PRINT@(9,1)MenuID”       “
PRINT@(9,3)ItemID”       “ ‘Else-Show menu info
RETURN  ‘Back to the “Loop”

Wake up, Dave. And try to spend a little time developing your programming skills and a little less time throwing rocks.

Regards

Chris Stasny

Well, Chris it seems you have a little chip on your shoulder about something (Do you work for Zedcor?). In my defense, there are many “programmers” out there that have had trouble getting ZBasic event processing to do everything they want it to do. It is true that some of the problems have been fixed as the users have been debugging ZBasic ever since it was released! In my opinion, there is no Basic available today for the Macintosh which is satisfactory for doing serious software development as there are glaring holes in their capabilities. I’ve had too many phone calls from disgruntled developers trying to find a solution to their problems with no solution in sight. Call it rock throwing if you like, but face it, there isn’t any such thing as LightSpeed Basic. Thank you for your letter. We encourage others to share the technology by writing for MacTutor. We are dedicated to the distribution of useful programming information without regard to race, creed or developer status. Chris, we never claim to know everything. We could use more of your good ideas. Write to MacTutor and ask for our authors kit and share some of your “expertise”.

Dave Kelly

I’m a new subscriber to MacTutor and have particularly enjoyed your column “Basic School”. I’m a BASIC programmer and frankly don’t have the time to really learn the other high level languages such as C or Pascal. It is, therefore, heartening to find a source of instruction for BASIC on the MAC.

Do you think it would be feasible to compile your articles into a separate volume? I believe such a book would sell quite well. There are many recreational and small application programmers who both love BASIC and the MAC. Keep up the great work.

Sincerely,

Julian Wan

Thank you for your kind remarks. I’ve had several people ask for a separate volume of “Basic School”. Because of the cost of publishing another compiled book and since we already offer the “Best of MacTutor” Vol 1 and 2 for a reasonable price, “Basic School” will not become a volume of its own. Have heart though, I’ve used information from C, Pascal, and assembly language columns many times and so having the other language available in “Best of MacTutor” can actually help when writing BASIC programs too. I’ve also compiled a MacTutor Index HyperCard Stack which can help when you are trying to learn about a specific subject.

Dave Kelly

 

Community Search:
MacTech Search:

Software Updates via MacUpdate

Direct Mail 5.2.1 - Create and send grea...
Direct Mail is an easy-to-use, fully-featured email marketing app purpose-built for macOS. Create, send, and track great looking email campaigns that get results. Start your newsletter by selecting... Read more
Skim 1.4.32 - PDF reader and note-taker...
Skim is a PDF reader and note-taker for OS X. It is designed to help you read and annotate scientific papers in PDF, but is also great for viewing any PDF file. Skim includes many features and has a... Read more
ForkLift 3.1.1 - Powerful file manager:...
ForkLift is a powerful file manager and ferociously fast FTP client clothed in a clean and versatile UI that offers the combination of absolute simplicity and raw power expected from a well-executed... Read more
MarsEdit 4.0.5 - Quick and convenient bl...
MarsEdit is a blog editor for OS X that makes editing your blog like writing email, with spell-checking, drafts, multiple windows, and even AppleScript support. It works with with most blog services... Read more
MarsEdit 4.0.5 - Quick and convenient bl...
MarsEdit is a blog editor for OS X that makes editing your blog like writing email, with spell-checking, drafts, multiple windows, and even AppleScript support. It works with with most blog services... Read more
ffWorks 1.0.1 - Convert multimedia files...
ffWorks (was iFFmpeg), focused on simplicity, brings a fresh approach to the use of FFmpeg, allowing you to create ultra-high-quality movies without the need to write a single line of code on the... Read more
Civilization VI 1.0.7 - Next iteration o...
Sid Meier’s Civilization VI is the next entry in the popular Civilization franchise. Originally created by legendary game designer Sid Meier, Civilization is a strategy game in which you attempt to... Read more
Printopia 3.0.5 - Share Mac printers wit...
Run Printopia on your Mac to share its printers to any capable iPhone, iPad, or iPod Touch. Printopia will also add virtual printers, allowing you to save print-outs to your Mac and send to apps.... Read more
Hopper Disassembler 4.3.12- - Binary dis...
Hopper Disassembler is a binary disassembler, decompiler, and debugger for 32- and 64-bit executables. It will let you disassemble any binary you want, and provide you all the information about its... Read more
LaunchBar 6.9.4 - Powerful file/URL/emai...
LaunchBar is an award-winning productivity utility that offers an amazingly intuitive and efficient way to search and access any kind of information stored on your computer or on the Web. It provides... Read more

Latest Forum Discussions

See All

Why doesn't mobile gaming have its...
The Overwatch League is a pretty big deal. It's an attempt to really push eSports into the mainstream, by turning them into, well, regular sports. But slightly less sweaty. It's a lavish affair with teams from all around the world, and more... | Read more »
Give Webzen’s new billiard game PoolTime...
Best known for producing hugely popular MMO titles, South Korean publisher Webzen is now taking aim at a different genre altogether. PoolTime is a realistic eight ball pool simulator, allowing you to compete in real-time matches against players... | Read more »
Let Them Come Guide - How to survive aga...
Let Them Come is all about making it as far as possible against overwhelming odds. Check out some of these tips to help you last a little longer in your unwinnable fight: [Read more] | Read more »
All the best games on sale for iPhone an...
Happy last day of the week. I hope you've been having a good one. I have. I saw ten doggos today. So because I'm in a good mood, I thought I'd round up all of the best games that are currently on sale on the App Store. [Read more] | Read more »
The very best games that came out for iP...
We're getting to the end of the first real, full, proper week of 2018. And in that time we've seen some pretty awesome games landing on the App Store. Of course, we've seen some absolute duffers as well. The sort of games that you look at and... | Read more »
Rusty Lake Paradise (Games)
Rusty Lake Paradise 1.4 Device: iOS Universal Category: Games Price: $2.99, Version: 1.4 (iTunes) Description: Jakob, the oldest son of the Eilander family, is returning to Paradise island after his mother passed away. Since her... | Read more »
Antihero Guide - Sneaky tricks to get ah...
Games of Antihero start out small and streamlined, but they quickly turn into long strategic conquests as you fight for control of the Victorian-era streets. If you find yourself struggling in the skullduggery department, here are a few things you... | Read more »
Here's why Niantic pulling Pokemon...
If there's one thing that Pokemon GO did well, it was bringing people together. I still remember seeing groups of people around the marina near where I live in the weeks after the game came out, all of them trying to grab some water Pokemon. There... | Read more »
Let Them Come (Games)
Let Them Come 1.0 Device: iOS Universal Category: Games Price: $1.99, Version: 1.0 (iTunes) Description: | Read more »
Why the iPhone X isn't the best bet...
It's always tempting to have the best thing. As game players we're almost conditioned to chase the biggest processor, the most power, the next generation. Even mobile gaming isn't immune to this push to progress - there are games on the App Store... | Read more »

Price Scanner via MacPrices.net

Amazon offers Silver 13″ Apple MacBook Pros f...
Amazon has new Silver 2017 13″ #Apple #MacBook Pros on sale today for up to $150 off MSRP, each including free shipping: – 13″ 2.3GHz/128GB Silver MacBook Pro (MPXR2LL/A): $1199.99 $100 off MSRP – 13... Read more
Sale: 12″ 1.3GHz MacBooks on sale for $1499,...
B&H Photo has Space Gray and Rose Gold 12″ 1.3GHz #Apple MacBooks on sale for $100 off MSRP. Shipping is free, and B&H charges sales tax for NY & NJ residents only: – 12″ 1.3GHz Space... Read more
Apple offers Certified Refurbished 2017 iMacs...
Apple has a full line of Certified Refurbished iMacs available for up to $350 off original MSRP. Apple’s one-year warranty is standard, and shipping is free. The following models are available: – 27... Read more
13″ MacBook Airs on sale for $120-$100 off MS...
B&H Photo has 2017 13″ 128GB MacBook Airs on sale for $120 off MSRP. Shipping is free, and B&H charges sales tax for NY & NJ residents only: – 13″ 1.8GHz/128GB MacBook Air (MQD32LL/A): $... Read more
15″ Touch Bar MacBook Pros on sale for up to...
Adorama has Space Gray 15″ MacBook Pros on sale for $200 off MSRP. Shipping is free, and Adorama charges sales tax in NJ and NY only: – 15″ 2.8GHz MacBook Pro Space Gray (MPTR2LL/A): $2199, $200 off... Read more
21″ 3.4GHz 4K iMac on sale for $1399, $100 of...
Adorama has the 21″ 3.4GHz 4K #Apple #iMac on sale today for $1399. Their price is $100 off MSRP. Shipping is free, and Adorama charges sales tax in NJ and NY only: – 21″ 3.4GHz 4K iMac (MNE02LL/A... Read more
B&H offering 13″ Apple MacBook Pros for u...
B&H Photo has 13″ MacBook Pros on sale for up to $75-$120 off MSRP. Shipping is free, and B&H charges sales tax for NY & NJ residents only: – 13-inch 2.3GHz/128GB Space Gray MacBook Pro (... Read more
B&H continues to offer clearance 2016 15″...
B&H Photo has clearance 2016 15″ #MacBook Pros available for up to $800 off original MSRP. Shipping is free, and B&H charges NY & NJ sales tax only: – 15″ 2.7GHz Touch Bar MacBook Pro... Read more
The cheapest 15″ Apple MacBook Pro available...
B&H Photo has the 15″ 2.2GHz MacBook Pro available for $200 off MSRP including free shipping plus NY & NJ sales tax only: – 15″ 2.2GHz MacBook Pro (MJLQ2LL/A): $1799 $200 off MSRP Apple has... Read more
Beats sale continues: $30-$80 off headphones,...
Walmart has Beats by Dr. Dre on sale on their online store for $30-$80 off MSRP (#Beats #Walmart), depending on the item: – Powerbeats3 Wireless Earphones: $134, save $65 – BeatsX Earphones: $109,... Read more

Jobs Board

*Apple* Retail - Multiple Positions - Apple,...
Job Description:SalesSpecialist - Retail Customer Service and SalesTransform Apple Store visitors into loyal Apple customers. When customers enter the store, Read more
Site Reliability Engineer, *Apple* Pay - Ap...
# Site Reliability Engineer, Apple Pay Job Number: 113356036 Santa Clara Valley, California, United States Posted: 12-Jan-2018 Weekly Hours: 40.00 **Job Summary** Read more
UI Tools and Automation Engineer, *Apple* M...
# UI Tools and Automation Engineer, Apple Media Products Job Number: 86351939 Santa Clara Valley, California, United States Posted: 11-Jan-2018 Weekly Hours: 40.00 Read more
*Apple* Retail - Multiple Positions - Apple,...
Job Description: Sales Specialist - Retail Customer Service and Sales Transform Apple Store visitors into loyal Apple customers. When customers enter the store, Read more
UI Tools and Automation Engineer, *Apple* M...
# UI Tools and Automation Engineer, Apple Media Products Job Number: 113136387 Santa Clara Valley, California, United States Posted: 11-Jan-2018 Weekly Hours: 40.00 Read more
All contents are Copyright 1984-2011 by Xplain Corporation. All rights reserved. Theme designed by Icreon.