<< Chapter < Page Chapter >> Page >
Problem Set for Publishing with MATLAB

Write a script to plot function y x x for 100 x 10 using increments of 100 . Publish your m-file to html.

The m-file content: % This script plots a graph of Graph of y=sin(x)/x clc % Clear screenx = pi/100:pi/100:10*pi; % Create a row vector y = sin(x)./x; % Calculate y as function of xplot(x,y),title('Graph of y=sin(x)/x'),xlabel('x'),ylabel('y'),grid The html output:

publish
The published html file.

A gas is expanded in an engine cylinder, following the law PV1.3=c. The initial pressure is 2550 kPa and the final pressure is 210 kPa. If the volume at the end of expansion is 0.75 m3, write a script to compute the work done by the gas and publish your solution to an html file. This is the same problem as this Problem you have solved before.

The m-file content: clc disp('A gas is expanded in an engine cylinder, following the law PV^1.3=c')disp('The initial pressure is 2550 kPa and the final pressure is 210 kPa.') disp('If the volume at the end of expansion is 0.75 m3,')disp('Compute the work done by the gas.') disp(' ') % Display blank linen=1.3; P_i=2550; % Initial pressureP_f=210; % Final pressure V_f=.75; % Final volumeV_i=(P_f*(V_f^n)/P_i)^(1/n); % Initial volume c=P_f*V_f^n;v=V_i:.001:V_f; % Creating a row vector for volume, v p=c./(v.^n); % Computing pressure for volumeWorkDone=trapz(v,p) % Integrating p*dv The html output:

publish
The published html file.

A force F acting on a body at a distance s from a fixed point is given by F 3 s 1 s 2 . Write a script to compute the work done when the body moves from the position where s=1 to that where s=10 and and publish your solution to an html file. Include a table of contents in the output file. This is the same problem as this Problem you have solved before.

The m-file content: %% Work done % This script computes the work done on an objectclc disp('A force F acting on a body at a distance s from a fixed point is given by')disp('F=3*s+(1/(s^2)) where s is the distance in meters') disp('Compute the total work done in moving')disp('From the position where s=1 to that where s=10.') disp(' ') % Display blank line%% Create a row vector for distance, s s=1:.001:10;%% Compute Force for s F=3.*s+(1./(s.^2)); % Computing Force for s%% Integrating F*ds over 1 to 10 meters. WorkDone=trapz(s,F) The html output:

publish
The published html file.

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, A brief introduction to engineering computation with matlab. OpenStax CNX. Nov 17, 2015 Download for free at http://legacy.cnx.org/content/col11371/1.11
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'A brief introduction to engineering computation with matlab' conversation and receive update notifications?

Ask