...

Plot σ_x vs. r_i

Calculate σ_x at O

Essentially, we need to repeat our previous calculation of σ_x at O for a range of r_ivalues. One way to do this is to loop over the previous calculation. There are two types of loops: for and while. Let's check the documentation to see what MATLAB provides for loops.

Help > Index (tab)

Search for loops.

There are two types of loops: for and while. Bring up the documentation for the for loop and look Enter for loop in the search bar and then select for which is the first entry up top. Look at the first example. This shows the format of the for loop. It is noted on this page that " You can often speed up the execution of MATLAB code by replacing for and while loops with vectorized code". That is what we'll do in Step 4. Let's stick it out with slowpoke for loop for now.

We'll create a new M-script file from beam.m and implement the for loop in that file.

Editor (tab) > File Save > Save As

Enter beam2.m for filename. Confirm that you are editing beam2.m (beam.m remains intact).

...

Check that your Workspace is cleared. Type help clear at the command line for help on this command. The command line help is a quick-and-easy way to check the documentation. But the information is not as extensive as the help documentation browser.

Clear the command window by entering clc

...

Our program reports sigma_x at the end of each iteration (since we didn't add a semi-colon after the sigma_x statement). From the output, we can deduce that sigma_x is a row vector, with the columns getting built up at each iteration. Note that a different way to look at the sigma_x array is to double-click on sigma_x in the Workspace. 

How do we know if the results are right? Let's spot-check the result for ri=1e-2 m against the previous step; this value of ri corresponds to k=6. What is your value of sigma_x(k) for k=6? You can check this from the output or by typing sigma_x(6) at the command line. My value is exactly the same as in Step 2 ... nailed that one! How about you?

Now that we have spot checked our result, add a semi-colon after the sigma_x statement to avoid cluttering up the output.

Plot σ_x vs. r_i

Let's pick MATLAB's brain on how to make a plot.

Help > ( contents tab) > matlab > getting started > graphics > using basic plotting functions

Home (tab) > Documentation (icon) > Search: Plotting Functions > Plotting Functions

Note that you can access this exact page without using the search field. Doing this will give you a sense of how to navigate the documentation. 

Home (tab) > Documentation (icon) > MATLAB > Graphics > 2-D and 3-D Plots > Plotting Basics > Plotting Functions

This page contains These help pages contain a lucid description of the extensive plotting capabilities available in MATLAB. When looking for plotting help, this would be the first place to turn to.

...

Click the Run icon in the Editor. The underscore symbol (_) gives subscripts and \sigma gives the corresponding Greek symbol. See

Help > Index (tab) > Greek letters and mathematical symbols

See my entire program here (right click and select save target as, or just left-click and copy-paste in the editor).

See how we are building up the program slowly and checking the results obsessively at each stage? This is a good programming approach and ends up being the most time-efficient especially for larger programs. So I'd highly recommend this approach.

Recall that "you can often speed up the execution of MATLAB code by replacing for and while loops with vectorized code". We are into speed, so let's vectorize our program.

Go to Step 4: Plot σ_x vs. r_i:  (Take 2See and rate the complete Learning Module)

Go to all MATLAB Learning Modules