int main () { int n = 100, i, ic; double fpi = 3.1415926 / 180.0, step, x; float xray[100], y1ray[100], y2ray[100]; step = 360. / (n - 1); for (i = 0; i < n; i++) { xray[i] = (float) (i * step); x = xray[i] * fpi; y1ray[i] = (float) sin (x); y2ray[i] = (float) cos (x); } metafl ("cons"); scrmod ("revers"); disini (); pagera (); complx (); axspos (450, 1800); axslen (2200, 1200); name ("X-axis", "x"); name ("Y-axis", "y"); labdig (-1, "x"); ticks (9, "x"); ticks (10, "y"); titlin ("Demonstration of CURVE", 1); titlin ("SIN(X), COS(X)", 3); ic = intrgb (0.95,0.95,0.95); axsbgd (ic); graf (0.0, 360.0, 0.0, 90.0, -1.0, 1.0, -1.0, 0.5); setrgb (0.7, 0.7, 0.7); grid (1, 1); color ("fore"); height (50); title (); color ("red"); curve (xray, y1ray, n); color ("green"); curve (xray, y2ray, n); disfin (); return 0; }
int main ( int argc, char *argv[] ) /******************************************************************************/ /* Purpose: CIRCLE_INOUT uses DISLIN to draw points in and out of a circle. Licensing: This code is distributed under the GNU LGPL license. Modified: 01 May 2011 Author: John Burkardt Reference: Helmut Michels, The Data Plotting Software DISLIN - version 10.4, Shaker Media GmbH, January 2010, ISBN13: 978-3-86858-517-9. */ { int i; int j; int m; int n_circle = 100; int n_in; int n_out; int pat; float pi = 3.14159265; float r; float s; float t; float *x; float *xy_in; float *xy_out; float *y; printf ( "\n" ); printf ( "CIRCLE_INOUT:\n" ); printf ( " C version\n" ); printf ( " Use DISLIN routines to make a scatterplot.\n" ); /* Read the data. */ r4mat_header_read ( "circle_in.txt", &m, &n_in ); xy_in = r4mat_data_read ( "circle_in.txt", m, n_in ); r4mat_header_read ( "circle_out.txt", &m, &n_out ); xy_out = r4mat_data_read ( "circle_out.txt", m, n_out ); /* Specify the format of the output file. */ metafl ( "png" ); /* Indicate that new data overwrites old data. */ filmod ( "delete" ); /* Specify the name of the output graphics file. */ setfil ( "circle_inout.png" ); /* Choose the page size and orientation. 'USA' is 2160 plot units wide and 2790 plot units high. 'P' requests PROFILE rather than LANDSCAPE orientation. */ setpag ( "usap" ); /* For PNG output, use reverse the default black background to white. */ scrmod ( "reverse" ); /* Open DISLIN. */ disini ( ); /* Plot a border around the page. */ pagera ( ); /* Use the COMPLEX font. */ complx ( ); /* Set the axis origin 230 plot units to the right, and 2560 plot units DOWN. */ axspos ( 230, 2560 ); /* Define the X and Y sizes of the axis system in plot units. */ axslen ( 1700, 1700 ); /* Label the X and Y axes. */ name ( "<--- X --->", "X" ); name ( "<--- Y --->", "Y" ); /* Relate the physical coordinates to the axes. */ graf ( 0.0, 1.0, 0.0, 0.1, 0.0, 1.0, 0.0, 0.1 ); /* Define axis system titles. */ titlin ( "Random points inside/outside the unit circle", 1 ); /* Draw the title. */ title ( ); /* Add a grid, with one grid line for every tick mark in the X and Y axes. */ grid ( 1, 1 ); /* Select the shading pattern. */ pat = 16; shdpat ( pat ); /* Set the color to blue. */ color ( "blue" ); /* At every data point, draw a circle of radius 0.01. */ for ( i = 0; i < n_in; i++ ) { rlcirc ( xy_in[0+i*2], xy_in[1+i*2], 0.01 ); } /* Set the color to red. */ color ( "red" ); /* At every data point, draw a circle of radius 0.01. */ for ( i = 0; i < n_out; i++ ) { rlcirc ( xy_out[0+i*2], xy_out[1+i*2], 0.01 ); } /* Draw a red circle. */ x = ( float * ) malloc ( n_circle * sizeof ( float ) ); y = ( float * ) malloc ( n_circle * sizeof ( float ) ); for ( i = 0; i < n_circle; i++ ) { t = pi * ( float ) ( i ) / ( float ) ( n_circle - 1 ); x[i] = cos ( t ); y[i] = sin ( t ); } thkcrv ( 10 ); curve ( x, y, n_circle ); /* End this plot. */ endgrf ( ); /* Free memory. */ free ( x ); free ( xy_in ); free ( xy_out ); free ( y ); /* Close DISLIN. */ disfin ( ); /* Terminate. */ printf ( "\n" ); printf ( "CIRCLE_INOUT:\n" ); printf ( " Normal end of execution.\n" ); return 0; }
int main ( int argc, char *argv[] ) /******************************************************************************/ /* Purpose: MAIN demonstrates the use of bar graphs. Modified: 09 April 2011 Reference: Helmut Michels, The Data Plotting Software DISLIN - version 10.4, Shaker Media GmbH, January 2010, ISBN13: 978-3-86858-517-9. */ { static char cbuf[25]; static char *ctit = "Bar Graphs (BARS)"; int i; int nya = 2700; static float x[9] = { 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0 }; static float y[9] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 }; static float y1[9] = { 1.0, 1.5, 2.5, 1.3, 2.0, 1.2, 0.7, 1.4, 1.1 }; static float y2[9] = { 2.0, 2.7, 3.5, 2.1, 3.2, 1.9, 2.0, 2.3, 1.8 }; static float y3[9] = { 4.0, 3.5, 4.5, 3.7, 4.0, 2.9, 3.0, 3.2, 2.6 }; printf ( "\n" ); printf ( "DISLIN_EX05:\n" ); printf ( " C version:\n" ); printf ( " Demonstrate the display of data in bar graphs.\n" ); /* Specify the format of the output file. */ metafl ( "png" ); /* Specify that if a file already exists of the given name, the new data should overwrite the old. */ filmod ( "delete" ); /* Specify the name of the output graphics file. */ setfil ( "dislin_ex05.png" ); /* Choose the page size and orientation. */ setpag ( "usap" ); /* For PNG output, reverse the default black background to white. */ scrmod ( "reverse" ); /* Open DISLIN. */ disini ( ); /* Plot a border around the page. */ pagera ( ); /* Use the COMPLEX font. */ complx ( ); ticks ( 1, "x" ); intax ( ); axslen ( 1600, 700 ); titlin ( ctit, 3 ); legini ( cbuf, 3, 8 ); leglin ( cbuf, "FIRST", 1 ); leglin ( cbuf, "SECOND", 2 ); leglin ( cbuf, "THIRD", 3 ); legtit ( " " ); shdpat ( 5L ); for ( i = 1; i <= 3; i++ ) { if ( 1 < i ) { labels ( "none", "x" ); } axspos ( 300, nya-(i-1)*800 ); graf ( 0.0, 10.0, 0.0, 1.0, 0.0, 5.0, 0.0, 1.0 ); if ( i == 1 ) { bargrp ( 3, 0.15 ); color ( "red" ); bars ( x, y, y1, 9 ); color ( "green" ); bars ( x, y, y2, 9 ); color ( "blue" ); bars ( x, y, y3, 9 ); color ( "fore" ); reset ( "bargrp" ); } else if ( i == 2 ) { height ( 30 ); labels ( "delta", "bars" ); labpos ( "center", "bars" ); color ( "red" ); bars ( x, y, y1, 9 ); color ( "green" ); bars ( x, y1, y2, 9 ); color ( "blue" ); bars ( x, y2, y3, 9 ); color ( "fore" ); reset ( "height" ); } else if ( i == 3 ) { labels ( "second", "bars" ); labpos ( "outside", "bars" ); color ( "red" ); bars ( x, y, y1, 9 ); color ( "fore" ); } if ( i < 3 ) { legend ( cbuf, 7 ); } else { height ( 50 ); title ( ); } endgrf ( ); } /* Close DISLIN. */ disfin ( ); /* Terminate. */ printf ( "\n" ); printf ( "DISLIN_EX05:\n" ); printf ( " Normal end of execution.\n" ); return 0; }
int main ( int argc, char *argv[] ) /******************************************************************************/ /* Purpose: ORBITAL uses DISLIN to display a contour plot of Z(X,Y) data. Licensing: This code is distributed under the GNU LGPL license. Modified: 20 May 2011 Author: John Burkardt Reference: Helmut Michels, The Data Plotting Software DISLIN - version 10.4, Shaker Media GmbH, January 2010, ISBN13: 978-3-86858-517-9. */ { int i; int j; int k; int level; int level_num; float level_value; int m; int n; int nn; float *x; float xmax; float xmin; float *xyz; float *y; float ymax; float ymin; float *z; float zmax; float zmin; printf ( "\n" ); printf ( "ORBITAL:\n" ); printf ( " C version:\n" ); printf ( " Use DISLIN to make a contour plot of Z(X,Y) data.\n" ); /* Read the data. */ r4mat_header_read ( "orbital.txt", &m, &nn ); xyz = r4mat_data_read ( "orbital.txt", m, nn ); /* Split the data. The contouring routine expects that data is along fixed X and Y coordinates, and so the X and Y data is to be given as vectors, not arrays. */ n = 101; x = ( float * ) malloc ( n * sizeof ( float ) ); y = ( float * ) malloc ( n * sizeof ( float ) ); z = ( float * ) malloc ( n * n * sizeof ( float ) ); k = 0; for ( i = 0; i < n; i++ ) { x[i] = xyz[0+k*3]; k = k + 1; } xmin = r4vec_min ( n, x ); xmax = r4vec_max ( n, x ); k = 0; for ( i = 0; i < n; i++ ) { y[i] = xyz[1+k*3]; k = k + n; } ymin = r4vec_min ( n, y ); ymax = r4vec_max ( n, y ); /* Z is a table. The first dimension should contain values for constant Y. */ k = 0; for ( i = 0; i < n; i++ ) { for ( j = 0; j < n; j++ ) { z[i+j*n] = xyz[2+k*3]; k = k + 1; } } zmax = r4mat_max ( n, n, z ); zmin = r4mat_min ( n, n, z ); /* Specify the format of the output file. */ metafl ( "png" ); /* Indicate that new data overwrites old data. */ filmod ( "delete" ); /* Specify the name of the output graphics file. */ setfil ( "orbital.png" ); /* Choose the page size and orientation. 'USA' is 2160 plot units wide and 2790 plot units high. 'P' requests PORTRAIT orientation. */ setpag ( "usap" ); /* For PNG output, reverse the default black background to white. */ scrmod ( "reverse" ); /* Open DISLIN. */ disini ( ); /* Plot a border around the page. */ pagera ( ); /* Use the SIMPLX font. */ simplx ( ); /* Set the axis origin in plot units to the right, and plot units DOWN. */ axspos ( 230, 2500 ); /* Define the X and Y sizes of the axis system in plot units. */ axslen ( 1700, 1700 ); /* Label the X and Y axes. */ name ( "X axis", "X" ); name ( "Y axis", "Y" ); /* Relate the physical coordinates to the axes, and specify tick marks. */ graf ( xmin, xmax, xmin, 1.0, ymin, ymax, ymin, 1.0 ); /* BEGIN LEVEL 2 COMMANDS. */ /* Define the title. */ titlin ( "Orbital contour plot", 1 ); title ( ); /* Set color to "blue". */ color ( "blue" ); /* Draw the contour plot. */ level_num = 10; for ( level = 1; level <= level_num; level++ ) { level_value = ( ( float ) ( level_num + 1 - level ) * zmin + ( float ) ( level ) * zmax ) / ( float ) ( level_num + 1 ); contur ( x, n, y, n, z, level_value ); } /* End this graph. */ endgrf ( ); /* RETURN FROM LEVEL 2 TO LEVEL 1. */ /* Close DISLIN. */ disfin ( ); /* Free memory. */ free ( x ); free ( xyz ); free ( y ); free ( z ); /* Terminate. */ printf ( "\n" ); printf ( "ORBITAL:\n" ); printf ( " Normal end of execution.\n" ); return 0; }