% m-file to demonstrate how to use a user-written function.
% Brian Hornbuckle, January 23, 2007.

clear;

temperature = linspace(0,20,1000);      % temperature vector, deg C

sat_vapor_pressure = sat_vp(temperature);       % use function to compute
                                                % saturation vapor pressure
                                                
% note:  the m-file that defines the function MUST be in the SAME
% DIRECTORY as the m-file that uses the function.  
                                                
figure(1)

plot(temperature,sat_vapor_pressure,'b-');
legend('saturation vapor pressure', 2);
xlabel('temperature, \circC');  ylabel('water vapor pressure, kPa');
axis([0 20 0 3]);  grid on;