### Chapter 5 Code Examples

### MATLAB Code Examples

1 2 3 |
x=[1 2 3 5 7 7.5 8 10]; y=[2 6.5 7 7 5.5 4 6 8]; plot(x,y) |

1 2 3 |
year=[1988:1:1994]; sales=[127 130 136 145 158 178 211]; plot(year, sales, '--r*') |

1 2 3 |
x=[-2: 0.01: 4]; y=3.5 .^ (-0.5*x) .* cos(6*x); plot(x,y) |

1 2 3 4 5 |
x=[-2:0.01:4]; y=3*x.^3-26*x+6; yd=9*x.^2-26; ydd=18*x; plot(x,y,'-b',x,yd,'--r',x,ydd,':k') |

1 2 3 4 5 6 7 8 |
x=[10:0.1:22]; y=95000./x.^2; x_data=[10:2:22]; y_data=[950 640 460 340 250 180 140]; plot(x,y,'-') hold on plot(x_data,y_data,'ro--') hold off |

1 2 3 4 5 |
xlabel('DISTANCE (cm)') ylabel('INTENSITY (lux)') title('Light Intensity as a Function of Distance') axis([8 24 0 1200]) legend('Theory','Experiment') |

### Python Code Examples

1 2 3 |
x=np.array([1,2,3,5,7,7.5,8,10]) y=np.array([2,6.5,7,7,5.5,4,6,8]) plt.plot(x,y) |

1 2 3 |
year=np.arange(1988,1994.1,1) sales=np.array([127,130,136,145,158,178,211]) plot(year, sales, '--r*') |

1 2 3 |
x=np.arange(-2,4.001,.01) y=3.5**(-0.5*x)*np.cos(6*x) plt.plot(x,y) |

1 2 3 4 5 |
x=np.arange(-2,4.001,.01) y=3*x**3-26*x+6 yd=9*x**2-26 ydd=18*x plt.plot(x,y,'-b',x,yd,'--r',x,ydd,':k') |

1 2 3 4 5 6 |
x=np.arange(10,22.01,0.1) y=95000/x**2; x_data=np.arange(10,22.01,2) y_data=np.array([950,640,460,340,250,180,140]) plt.plot(x,y,'-',label='Theory') plt.plot(x_data,y_data,'ro--',label='Experiment') |

1 2 3 4 5 |
plt.xlabel('DISTANCE (cm)') plt.ylabel('INTENSITY (lux)') plt.title('Light Intensity as a Function of Distance') plt.axis([8,24,0,1200]) plt.legend() |

### Exercises

- Make two separate plots of the function

one for x and one for x - Plot the function

Notice that the function has a vertical asymptote at x=2. Plot the function by creating two vectors for the domain of x, the first vector with elements from -4 to 1.7, and the second vector with elements from 2.3 to 8. For each x vector, create a y vector with the corresponding values of y according to the function. To plot the function make two curves in the same plot. - You have an electrical circuit that includes a voltage source vs with an internal resistance rs and a load resistance . The power P dissipated in the load is given by:

Plot the power P as a function of for given that and - The Gateway Arch in St. Louis is shaped according to the equation:

Make a plot of the arch. Draw a horizontal line at ground level (y=0),

- In astronomy, the relationship between the relative temperature (temperature relative to the sun), relative luminosity and relative radius of a star is modeled by:

The HR (Hertzsprung-Russell) diagram is a plot of versus the temperature. The following data are given:Sun Spica Regulus Alioth Barndard’sStar Epsilon Indi Beta Crucis Temp (K) 5840 22400 13260 9400 3130 4280 28200 L/L 1 13400 150 108 0 0.15 43000 R/R 1 7.8 3.5 3.7 0.18 0.76 8 To compare the data with the model, plot a HR diagram. The diagram should have two sets of points. One uses the values of L/L from the table (use asterisk markers), and the other uses values of L/L that are calculated by the equation using R/R from the table (use circle markers). In the HR diagram both axes are logarithmic. In addition, the values of the temperature on the horizontal axis are decreasing from left to right. Label the axes and use a legend. (See https://en.wikipedia.org/wiki/Hertzsprung%E2%80%93Russell_diagram for a general explanation of the HR diagram.)

- The position x as a function of time that a particle moves along a straight line is given by:

The velocity v(t) is determined by the derivative of x(t) with respect to t, and the acceleration a(t) is determined by the derivative of v(t) with respect to t.

Derive the expressions for the velocity and acceleration of the particle, and make plots for the position, velocity, and acceleration as a function of time for t Time t is measured in seconds, and position x is measured in meters. Make three plots, one for position, one for velocity, and one for acceleration. Label the axes appropriately with correct units. - A resistor, R=2 and an inductor, L=1.7H are connected to a voltage source in series (RL circuit). When the voltage source applies a rectangular voltage pulse with an amplitude of V=24V and a duration of 0.5s, the current i(t) in the circuit as a function of time is given by:

for t

for t

Make a plot of the current as a function of time for t seconds.