PROGRAM QA1
C CALIFORNIA COMPUTER PRODUCTS SAMPLE PROGRAM NUMBER 2      V097
CALCOMP  DRIVER PROGRAM INCH VERSION FOR FORTRAN 77.
CALCOMP HCBS FOR  906 CONTROLLERS                          SEPT, 1981 *
C             PROPRIETARY SOFTWARE NOTIFICATION
C  THIS SOFTWARE IS THE PROPRIETARY AND CONFIDENTIAL PROPERTY OF
C    CALIFORNIA COMPUTER PRODUCTS, INC. AND IS LICENSED FOR USE ON THE
C    DESIGNATED EQUIPMENT ON WHICH IT WAS ORIGINALLY INSTALLED AND
C    CANNOT BE MODIFIED, DUPLICATED OR COPIED IN ANY FORM WITHOUT PRIOR
C    WRITTEN CONSENT OF CALIFORNIA COMPUTER PRODUCTS, INC.
C  IF SUPPLIED UNDER A GOVERNMENT CONTRACT, THE FOLLOWING APPLIES: 
C                 RESTRICTED RIGHTS LEGEND
C  USE, DUPLICATION OR DISCLOSURE BY THE GOVERNMENT IS SUBJECT TO
C    RESTRICTIONS AS SET FORTH IN PARAGRAPH (B)(3)(B) OF THE RIGHTS
C    IN TECHNICAL DATA AND COMPUTER SOFTWARE CLAUSE IN DAR 7-104.9
C    (A).  CALIFORNIA COMPUTER PRODUCTS, INC., 2411 WEST LAPALMA
C    AVENUE, ANAHEIM, CALIFORNIA 92803.
      DATA LDEV /40 /
      CHARACTER * 10 ICHR1
      CHARACTER * 14 ICHR2,ICHR3
      CHARACTER * 5  ICHR4,ICHR6
      CHARACTER * 4  ICHR5,ICHR7
      CHARACTER * 29 ICHR8
      CHARACTER * 32 ICHR9,ICHR10,ICHR11,ICHR12
      CHARACTER * 1 IBCD
      ICHR1='WIDTH (FT)'
      ICHR2='THICKNESS (IN)'
      ICHR3='PRESSURE (PSI)'
      ICHR4='THK= '
      ICHR5=' IN.'
      ICHR6='WTH= '
      ICHR7=' FT.'
      ICHR8='CRITICAL BUCKLING PRESSURE OF'
      ICHR9='HYPERBOLIC PARABOLOID SHELLS FOR'
      ICHR10='FIXED WIDTH VS VARYING THICKNESS'
      ICHR11='FIXED THICKNESS VS VARYING WIDTH'
      ICHR12='PREPARED ON A CALCOMP PLOTTER'
      CALL C_PLOTS(0,0,LDEV)
C ESTABLISH AN ORIGIN SO NEGATIVE VALUES UP TO -0.5 MAY BE USED
      CALL C_PLOT(0.5,0.5,-3)
C PLOT X-AXIS FOR WIDTH
      X=0.0
      DO 210  I=1,10
        CALL C_PLOT(X,0.0,3)
        X=X+1.0
        CALL C_PLOT(X,0.0,2)
        CALL C_PLOT(X,-.1,2)
 210    CALL C_NUMBER(X,-0.25,0.1,5.0*X,0.0,-1)
      CALL C_SYMBOL(4.0,-0.40,0.12,IBCD,1,0.0,-1)
      CALL C_SYMBOL(4.2,-0.45,0.14,ICHR1,INTEQ,0.0,10)
      CALL C_PLOT (0.0,0.5,-3)
C PLOT X-AXIS FOR THICKNESS
      X=0.0
      DO 220  I=1,5
        CALL C_PLOT(X,0.0,3)
        X=X+1.0
        CALL C_PLOT(X,0.0,2)
        CALL C_PLOT(X,-.1,2)
        CALL C_PLOT(X,0.0,2)
        X=X+1.0
        CALL C_PLOT(X,0.0,2)
        CALL C_PLOT(X,-.1,2)
 220    CALL C_NUMBER(X,-0.25,0.1,X,0.0,-1)
      CALL C_SYMBOL(3.7,-0.40,0.12,IBCD,7,0.0,-1)
      CALL C_SYMBOL(4.0,-0.45,0.14,ICHR2,INTEQ,0.0,14)
C PLOT Y-AXIS
      Y=0.0
      DO 230  I=1,9
        CALL C_PLOT(0.0,Y,3)
        Y=Y+1.0
        CALL C_PLOT(0.0,Y,2)
        CALL C_PLOT(-.1,Y,2)
 230    CALL C_NUMBER(-.15,Y-.2,0.1,1000.*Y,90.0,0)
      CALL C_SYMBOL(-0.30,3.5,0.14,ICHR3,INTEQ,90.0,14)
      THICK=3.0
      WIDTH=25.0
      DO 260  I=1,3
        TSQR=THICK*THICK
        WSQR=WIDTH*WIDTH
        PSI=100.99*TSQR
        CALL C_SYMBOL(0.6,PSI/1000.0,0.1,ICHR4,INTEQ,0.0,5)
        CALL C_NUMBER(999.0,999.0,0.10,THICK,0.0,0)
        CALL C_SYMBOL(999.0,999.0,0.10,ICHR5,INTEQ,0.0,4)
        CALL C_SYMBOL( 2.0 ,999.0,0.12,IBCD,1,0.0,-1)
        DO 240 J=10,50
          WX=J
          PSI=10099.0*TSQR/(WX*WX)
 240      CALL C_PLOT(WX/5.0,PSI/1000.0,2)
        PSI=10099.0*81.0/WSQR
        CALL C_SYMBOL(9.2,PSI/1000.0,0.1,ICHR6,INTEQ,0.0,5)
        CALL C_NUMBER(999.0,999.0,0.10,WIDTH,0.0,0)
        CALL C_SYMBOL(999.0,999.0,0.10,ICHR7,INTEQ,0.0,4)
        CALL C_SYMBOL( 9.0 ,999.0,0.12,IBCD,7,0.0,-1)
        DO 250 J=5,50
          TX=J
          TX=(50.0-TX)/5.0
          PSI=10099.0*TX*TX/WSQR
 250      CALL C_PLOT(TX,PSI/1000.0,2)
        THICK=THICK+3.0
 260    WIDTH=WIDTH-5.0
      CALL C_SYMBOL(3.3,8.5,.14,ICHR8,INTEQ,0.,29)
      CALL C_SYMBOL(3.1,8.2,.14,ICHR9,INTEQ,0.,32)
      CALL C_SYMBOL(3.1,7.9,.14,ICHR10,INTEQ,0.,32)
      CALL C_SYMBOL(3.1,7.6,.14,ICHR11,INTEQ,0.,32)
      CALL C_SYMBOL(3.3,7.0,.14,ICHR12,INTEQ,0.,29)
      CALL C_PLOT(0.0,0.0,999)
      STOP
      END