Re: Two Superposed Plots

[Martin Senator]

Thanks to everyone who helped.

This file is to document a solution
that worked.


   1 + echo on;
   2 +
   3 + # Just issued the commands
   4 + # > page_screen_output=0; > diary on; and > echo on;
   5 + # file superpose.m ;
   6 + # purpose: to show how to superpose two plots in Octave.
   7 +
   8 +
   9 + omega = linspace (0, 2.5, 501) ;
  10 + a     = 1./sqrt( (-omega.^2 + 1 ).^2 + (0.64*omega).^2 ) ;
  11 + mat1 = [omega' a' ]  ;
  12 +
  13 + th_circ = linspace ( 0, 2*pi, 361 ) ;
  14 + x_circ  = cos ( th_circ ) + 0.4 ;
  15 + y_circ  = sin ( th_circ ) ;
  16 + mat_circ = [x_circ', y_circ'] ;
  17 +
  18 + titlestr = "Illustrating Superposition of Two Plots" ;
  19 + xlabelstr = "Forcing Frequency, {/Symbol w}" ;
  20 + ylabelstr = "Amplitude, a"
  21 ylabelstr = Amplitude, a
  22 +
  23 +    __gnuplot_set__ terminal postscript portrait enhanced color
  24 +    __gnuplot_set__ output 'plotsup.ps'
  25 +    __gnuplot_set__ nokey ;
  26 +    __gnuplot_set__ timestamp bottom  ;
  27 +
  28 +    __gnuplot_set__ xtics nomirror ;
  29 +    __gnuplot_set__ ytics nomirror ;
  30 +    __gnuplot_set__ xtics  ;
  31 +    __gnuplot_set__ ytics  ;
  32 +    __gnuplot_set__ x2tics  ;
  33 +    __gnuplot_set__ y2tics  ;
  34 +    __gnuplot_set__ xrange [0:2.5]  ;
  35 +    __gnuplot_set__ yrange [0:1.8]  ;
  36 +    __gnuplot_set__ x2range [-3:1.5]  ;
  37 +    __gnuplot_set__ y2range [-3.0:1.5]  ;
  38 +    __gnuplot_set__ size square  ;
  39 + #
  40 +    __gnuplot_plot__ mat1(:,1:2) axes x1y1\
  41 +             ,mat_circ(:,1:2) axes x2y2  ;
  42 + #
  43 + xlabel ( xlabelstr ) ;
  44 + ylabel ( ylabelstr ) ;
  45 + title  ( titlestr )  ;
  46 +    __gnuplot_set__ x2label "Displacement"  ;
  47 +    __gnuplot_set__ y2label "Velocity"  ;
  48 + #
  49 +      system('rm  plotsup.ps');
  50 +      __gnuplot_set__ output 'plotsup.ps'
  51 +      replot
  52 +      closeplot
  53 +
  54
  55 + diary off;

ms