Plot a function represented as x, y numerical arrays.
Post the resulting image for the following input arrays (taken fromPython's Example section onTime a function):
x = {0, 1, 2, 3, 4, 5, 6, 7, 8, 9}; y = {2.7, 2.8, 31.4, 38.1, 58.0, 76.2, 100.5, 130.0, 149.3, 180.0};This task is intended as a subtask forMeasure relative performance of sorting algorithms implementations.
/* ARM assembly AARCH64 Raspberry PI 3B *//* program areaPlot64.s */ /*******************************************//* Constantes file *//*******************************************//* for this file see task include a file in language AArch64 assembly*/.include "../includeConstantesARM64.inc".equ HAUTEUR, 22.equ LARGEUR, 50.equ MARGEGAUCHE, 10/*******************************************//* Structures *//********************************************//* structure for points */ .struct 0point_posX: .struct point_posX + 8 point_posY: .struct point_posY + 8 point_end:/*******************************************//* Initialized data *//*******************************************/.dataszMessError: .asciz "Number of points too large !! \n"szCarriageReturn: .asciz "\n"szMessMovePos: .ascii "\033[" // cursor positionposY: .byte '0' .byte '6' .ascii ";"posX: .byte '0' .byte '3' .asciz "H*"szMessEchelleX: .asciz "Y^ X="szClear1: .byte 0x1B .byte 'c' // other console clear .byte 0szMessPosEch: .ascii "\033[" // scale cursor positionposY1: .byte '0' .byte '0' .ascii ";"posX1: .byte '0' .byte '0' .asciz "H"//x = {0, 1, 2, 3, 4, 5, 6, 7, 8, 9};//y = {2.7, 2.8, 31.4, 38.1, 58.0, 76.2, 100.5, 130.0, 149.3, 180.0}; /* areas points */tbPoints: .quad 0 // 1 .quad 27 // Data * 10 for integer operation .quad 1 // 2 .quad 28 .quad 2 // 3 .quad 314 .quad 3 // 4 .quad 381 .quad 4 // 5 .quad 580 .quad 5 // 6 .quad 762 .quad 6 // 7 .quad 1005 .quad 7 // 8 .quad 1300 .quad 8 // 9 .quad 1493 .quad 9 // 10 .quad 1800/*******************************************//* UnInitialized data *//*******************************************/.bsssZoneConv: .skip 30/*******************************************//* code section *//*******************************************/.text.global main main: // entry of program ldr x0,qAdrtbPoints // area address mov x1,10 // size mov x2,LARGEUR mov x3,HAUTEUR bl plotArea b 100f100: // standard end of the program mov x0, 0 // return code mov x8,EXIT // request to exit program svc 0 // perform the system callqAdrsZoneConv: .quad sZoneConvqAdrszCarriageReturn: .quad szCarriageReturnqAdrtbPoints: .quad tbPoints/************************************/ /* create graph *//************************************/ /* x0 contains area points address *//* x1 contains number points *//* x2 contains graphic weight *//* x3 contains graphic height *//* REMARK : no save x9-x20 registers */plotArea: stp x2,lr,[sp,-16]! // save registers stp x3,x4,[sp,-16]! // save registers cmp x1,x2 bge 99f mov x9,x0 mov x4,x1 ldr x10,qAdrposX ldr x11,qAdrposY mov x12,#0 // indice mov x13,point_end // element area size mov x17,0 // Y maxi mov x19,-1 // Y Mini1: //search mini maxi madd x14,x12,x13,x0 // load coord Y ldr x15,[x14,point_posY] cmp x15,x17 csel x17,x15,x17,hi // maxi ? cmp x15,x19 csel x19,x15,x19,lo // mini ? add x12,x12,#1 cmp x12,x1 // end ? blt 1b // no -> loop // compute ratio udiv x15,x17,x3 // ratio = maxi / height add x15,x15,1 // for adjust ldr x0,qAdrszClear1 // clear screen bl affichageMess udiv x20,x2,x4 // compute interval X = weight / number points mov x12,0 // indice2: // loop begin for display point madd x14,x12,x13,x9 // charge X coord point ldr x16,[x14,point_posX] mul x16,x20,x12 // interval * indice add x0,x16,MARGEGAUCHE // + left margin mov x1,x10 // conversion ascii and store bl convPos ldr x18,[x14,point_posY] // charge Y coord point udiv x18,x18,x15 // divide by ratio sub x0,x3,x18 // inversion position ligne mov x1,x11 // conversion ascii and store bl convPos ldr x0,qAdrszMessMovePos // display * at position X,Y bl affichageMess add x12,x12,1 // next point cmp x12,x4 // end ? blt 2b // no -> loop // display left scale // display Y Mini mov x0,0 ldr x1,qAdrposX1 bl convPos mov x0,HAUTEUR ldr x1,qAdrposY1 bl convPos ldr x0,qAdrszMessPosEch bl affichageMess mov x0,x19 ldr x1,qAdrsZoneConv bl conversion10 ldr x0,qAdrsZoneConv bl affichageMess // display Y Maxi mov x0,0 ldr x1,qAdrposX1 bl convPos mov x0,0 ldr x1,qAdrposY1 bl convPos ldr x0,qAdrszMessPosEch bl affichageMess mov x0,x17 ldr x1,qAdrsZoneConv bl conversion10 ldr x0,qAdrsZoneConv bl affichageMess // display average value mov x0,0 ldr x1,qAdrposX1 bl convPos mov x0,HAUTEUR/2 add x0,x0,#1 ldr x1,qAdrposY1 bl convPos ldr x0,qAdrszMessPosEch bl affichageMess lsr x0,x17,#1 ldr x1,qAdrsZoneConv bl conversion10 ldr x0,qAdrsZoneConv bl affichageMess // display X scale mov x0,0 ldr x1,qAdrposX1 bl convPos mov x0,HAUTEUR+1 ldr x1,qAdrposY1 bl convPos ldr x0,qAdrszMessPosEch bl affichageMess ldr x0,qAdrszMessEchelleX bl affichageMess mov x12,0 // indice mov x19,MARGEGAUCHE10: udiv x20,x2,x4 madd x0,x20,x12,x19 ldr x1,qAdrposX1 bl convPos mov x0,HAUTEUR+1 ldr x1,qAdrposY1 bl convPos ldr x0,qAdrszMessPosEch bl affichageMess madd x14,x12,x13,x9 // load X coord point ldr x0,[x14,point_posX] ldr x1,qAdrsZoneConv bl conversion10 ldr x0,qAdrsZoneConv bl affichageMess add x12,x12,1 cmp x12,x4 blt 10b ldr x0,qAdrszCarriageReturn bl affichageMess mov x0,0 // return code b 100f99: // error ldr x0,qAdrszMessError bl affichageMess mov x0,-1 // return code100: ldp x3,x4,[sp],16 // restaur 2 registers ldp x2,lr,[sp],16 // restaur 2 registers ret // return to address lr x30qAdrszMessMovePos: .quad szMessMovePosqAdrszClear1: .quad szClear1qAdrposX: .quad posXqAdrposY: .quad posYqAdrposX1: .quad posX1qAdrposY1: .quad posY1qAdrszMessEchelleX: .quad szMessEchelleXqAdrszMessPosEch: .quad szMessPosEchqAdrszMessError: .quad szMessError/************************************/ /* conv position in ascii and store at address *//************************************/ /* x0 contains position *//* x1 contains string address */convPos: stp x2,lr,[sp,-16]! // save registers stp x3,x4,[sp,-16]! // save registers mov x2,10 udiv x3,x0,x2 add x4,x3,48 // convert in ascii strb w4,[x1] // store posX msub x4,x3,x2,x0 add x4,x4,48 strb w4,[x1,1]100: ldp x3,x4,[sp],16 // restaur 2 registers ldp x2,lr,[sp],16 // restaur 2 registers ret // return to address lr x30/********************************************************//* File Include fonctions *//********************************************************//* for this file see task include a file in language AArch64 assembly */.include "../includeARM64.inc"1800 * * * *900 * * * *27 * *Y^ X= 0 1 2 3 4 5 6 7 8 9
INCLUDE "D2:REAL.ACT" ;from the Action! Tool KitDEFINE PTR="CARD"DEFINE BUF_SIZE="100"DEFINE REAL_SIZE="3"TYPE Settings=[ INT xMin,xMax,xStep,yMin,yMax,yStep INT xLeft,xRight,yTop,yBottom INT tickLength]BYTE ARRAY xs(BUF_SIZE),ys(BUF_SIZE)BYTE count=[0]PTR FUNC GetXPtr(BYTE i)RETURN (xs+3*i)PTR FUNC GetYPtr(BYTE i)RETURN (ys+3*i)PROC AddPoint(CHAR ARRAY xstr,ystr) REAL POINTER p p=GetXPtr(count) ValR(xstr,p) p=GetYPtr(count) ValR(ystr,p) count==+1RETURNPROC InitData() AddPoint("0.0","2.7") AddPoint("1.0","2.8") AddPoint("2.0","31.4") AddPoint("3.0","38.1") AddPoint("4.0","58.0") AddPoint("5.0","76.2") AddPoint("6.0","100.5") AddPoint("7.0","130.0") AddPoint("8.0","149.3") AddPoint("9.0","180.0")RETURNINT FUNC GetXPos(Settings POINTER s INT x) INT res res=x*(s.xRight-s.xLeft)/(s.xMax-s.xMin)+s.xLeftRETURN (res)INT FUNC GetYPos(Settings POINTER s INT y) INT res res=y*(s.yTop-s.yBottom)/(s.yMax-s.yMin)+s.yBottomRETURN (res)INT FUNC GetXPosR(Settings POINTER s REAL POINTER x) REAL nom,denom,div,tmp INT res IntToReal(s.xRight-s.xLeft,tmp) RealMult(tmp,x,nom) IntToReal(s.xMax-s.xMin,denom) RealDiv(nom,denom,div) res=RealToInt(div)+s.xLeftRETURN (res)INT FUNC GetYPosR(Settings POINTER s REAL POINTER y) REAL nom,denom,div,tmp INT res IntToReal(s.yBottom-s.yTop,tmp) RealMult(tmp,y,nom) IntToReal(s.yMax-s.yMin,denom) RealDiv(nom,denom,div) res=-RealToInt(div)+s.yBottomRETURN (res)BYTE FUNC AtasciiToInternal(CHAR c) BYTE c2 c2=c&$7F IF c2<32 THEN RETURN (c+64) ELSEIF c2<96 THEN RETURN (c-32) FIRETURN (c)PROC CharOut(INT x,y CHAR c) BYTE i,j,v PTR addr addr=$E000+AtasciiToInternal(c)*8; FOR j=0 TO 7 DO v=Peek(addr) i=8 WHILE i>0 DO IF v&1 THEN Plot(x+i,y+j) FI v=v RSH 1 i==-1 OD addr==+1 ODRETURNPROC TextOut(INT x,y CHAR ARRAY text) BYTE i FOR i=1 TO text(0) DO CharOut(x,y,text(i)) x==+8 ODRETURNPROC DrawAxes(Settings POINTER s) INT i,x,y CHAR ARRAY t(10) Plot(s.xLeft,s.yTop) DrawTo(s.xLeft,s.yBottom) DrawTo(s.xRight,s.yBottom) FOR i=s.xMin TO s.xMax STEP s.xStep DO x=GetXPos(s,i) Plot(x,s.yBottom) DrawTo(x,s.yBottom+s.tickLength) StrI(i,t) TextOut(x-t(0)*4,s.yBottom+s.tickLength+1,t) OD FOR i=s.yMin TO s.yMax STEP s.yStep DO y=GetYPos(s,i) Plot(s.xLeft-s.tickLength,y) DrawTo(s.xLeft,y) StrI(i,t) TextOut(s.xLeft-s.tickLength-1-t(0)*8,y-4,t) ODRETURNPROC DrawPoint(INT x,y) Plot(x-1,y-1) DrawTo(x+1,y-1) DrawTo(x+1,y+1) DrawTo(x-1,y+1) DrawTo(x-1,y-1)RETURNPROC DrawSeries(Settings POINTER s) INT i,x,y,prevX,prevY REAL POINTER p FOR i=0 TO count-1 DO p=GetXPtr(i) x=GetXPosR(s,p) p=GetYPtr(i) y=GetYPosR(s,p) DrawPoint(x,y) IF i>0 THEN Plot(prevX,prevY) DrawTo(x,y) FI prevX=x prevY=y ODRETURNPROC DrawPlot(Settings POINTER s) DrawAxes(s) DrawSeries(s)RETURNPROC Main() BYTE CH=$02FC,COLOR1=$02C5,COLOR2=$02C6 Settings s Graphics(8+16) Color=1 COLOR1=$0C COLOR2=$02 InitData() s.xMin=0 s.xMax=9 s.xStep=1 s.yMin=0 s.yMax=180 s.yStep=20 s.xLeft=30 s.xRight=311 s.yTop=8 s.yBottom=177 s.tickLength=3 DrawPlot(s) DO UNTIL CH#$FF OD CH=$FFRETURNScreenshot from Atari 8-bit computer
Like C, this is often outsourced to another program like gnuplot, but is also possible with GtkAda.
withGtk.Main;withGtk.Window;useGtk.Window;withGtk.Widget;useGtk.Widget;withGtk.Handlers;useGtk.Handlers;withGlib;useGlib;withGtk.Extra.Plot;useGtk.Extra.Plot;withGtk.Extra.Plot_Data;useGtk.Extra.Plot_Data;withGtk.Extra.Plot_Canvas;useGtk.Extra.Plot_Canvas;withGtk.Extra.Plot_Canvas.Plot;useGtk.Extra.Plot_Canvas.Plot;procedurePlotCoordsispackageHandleris newCallback(Gtk_Widget_Record);Window:Gtk_Window;Plot:Gtk_Plot;PCP:Gtk_Plot_Canvas_Plot;Canvas:Gtk_Plot_Canvas;PlotData:Gtk_Plot_Data;x,y,dx,dy:Gdouble_Array_Access;procedureExitMain(Object:accessGtk_Widget_Record'Class)isbeginDestroy(Object);Gtk.Main.Main_Quit;endExitMain;beginx:=newGdouble_Array'(0.0,1.0,2.0,3.0,4.0,5.0,6.0,7.0,8.0,9.0);y:=newGdouble_Array'(2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0);Gtk.Main.Init;Gtk_New(Window);Set_Title(Window,"Plot coordinate pairs with GtkAda");Gtk_New(PlotData);Set_Points(PlotData,x,y,dx,dy);Gtk_New(Plot);Add_Data(Plot,PlotData);Autoscale(Plot);Show(PlotData);Hide_Legends(Plot);Gtk_New(PCP,Plot);Show(Plot);Gtk_New(Canvas,500,500);Show(Canvas);Put_Child(Canvas,PCP,0.15,0.15,0.85,0.85);Add(Window,Canvas);Show_All(Window);Handler.Connect(Window,"destroy",Handler.To_Marshaller(ExitMain'Access));Gtk.Main.Main;endPlotCoords;
File: Plot_coordinate_pairs.a68
#!/usr/bin/algol68g-full --script ## -*- coding: utf-8 -*- #PR READ "prelude/errata.a68" PR;PR READ "prelude/exception.a68" PR;PR READ "prelude/math_lib.a68" PR;CO REQUIRED BY "prelude/graph_2d.a68" CO MODE GREAL= REAL; # single precision # FORMAT greal repr = $g(-3,0)$;PR READ "prelude/graph_2d.a68" PR;[]REAL x = (0, 1, 2, 3, 4, 5, 6, 7, 8, 9);[]REAL y = (2.7, 2.8, 31.4, 38.1, 58.0, 76.2, 100.5, 130.0, 149.3, 180.0);test:( REF GRAPHDD test graph = INIT LOC GRAPHDD; type OF window OF test graph := "gif"; # or gif, ps, X, pnm etc # title OF test graph := "Plot coordinate pairs"; sub title OF test graph := "Algol68"; interval OF (axis OF test graph)[x axis] := (0, 8); label OF (axis OF test graph)[x axis] := "X axis"; interval OF (axis OF test graph)[y axis] := (0, 200); label OF (axis OF test graph)[y axis] := "Y axis"; PROC curve = (POINTYIELD yield)VOID: FOR i TO UPB x DO yield((x[i],y[i])) OD; (begin curve OF (METHODOF test graph))(~); (add curve OF (METHODOF test graph))(curve, (red,solid)); (end curve OF (METHODOF test graph))(~));PR READ "postlude/exception.a68" PR
Image - Link, since uploads seem to be disabled currently.
(AutoHotkey1.1+)
#SingleInstance,Force#NoEnvSetBatchLines,-1OnExit,ExitFileOut:=A_Desktop"\MyNewFile.png"Font:="Arial"x:=[0,1,2,3,4,5,6,7,8,9]y:=[2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0]; Uncomment if Gdip.ahk is not in your standard library; #Include, Gdip.ahkif(!pToken:=Gdip_Startup()){MsgBox,48,Gdipluserror!,Gdiplusfailedtostart.PleaseensureyouhaveGdiplusonyoursystem.ExitApp}If(!Gdip_FontFamilyCreate(Font)){MsgBox,48,Fonterror!,Thefontyouhavespecifieddoesnotexistonyoursystem.ExitApp}pBitmap:=Gdip_CreateBitmap(900,900),G:=Gdip_GraphicsFromImage(pBitmap),Gdip_SetSmoothingMode(G,4),pBrush:=Gdip_BrushCreateSolid(0xff000000),Gdip_FillRectangle(G,pBrush,-3,-3,906,906),Gdip_DeleteBrush(pBrush),pPen1:=Gdip_CreatePen(0xffffcc00,2),pPen2:=Gdip_CreatePen(0xffffffff,2),pPen3:=Gdip_CreatePen(0xff447821,1),pPen4:=Gdip_CreatePen(0xff0066ff,2),Gdip_DrawLine(G,pPen2,50,50,50,850),Gdip_DrawLine(G,pPen2,50,850,850,850),FontOptions1:="x0 y870 Right cbbffffff r4 s16 Bold",Gdip_TextToGraphics(G,0,FontOptions1,Font,40,20)Loop,%x.MaxIndex()-1{Offset1:=50+(x[A_Index]*80),Offset2:=Offset1+80,Gdip_DrawLine(G,pPen1,Offset1,850-(y[A_Index]*4),Offset1+80,850-(y[A_Index+1]*4))}Loop,%x.MaxIndex(){Offset1:=50+((A_Index-1)*80),Offset2:=Offset1+80,Offset3:=45+(x[A_Index]*80),Offset4:=845-(y[A_Index]*4),Gdip_DrawLine(G,pPen2,45,Offset1,55,Offset1),Gdip_DrawLine(G,pPen2,Offset2,845,Offset2,855),Gdip_DrawLine(G,pPen3,50,Offset1,850,Offset1),Gdip_DrawLine(G,pPen3,Offset2,50,Offset2,850),Gdip_DrawLine(G,pPen4,Offset3,Offset4,Offset3+10,Offset4+10),Gdip_DrawLine(G,pPen4,Offset3,Offset4+10,Offset3+10,Offset4),FontOptions1:="x0 y"(Offset1-7)" Right cbbffffff r4 s16 Bold",FontOptions2:="x"(Offset2-7)" y870 Left cbbffffff r4 s16 Bold",Gdip_TextToGraphics(G,220-(A_Index*20),FontOptions1,Font,40,20),Gdip_TextToGraphics(G,A_Index,FontOptions2,Font,40,20)}Gdip_DeletePen(pPen1),Gdip_DeletePen(pPen2),Gdip_DeletePen(pPen3),Gdip_DeletePen(pPen4),Gdip_SaveBitmapToFile(pBitmap,FileOut),Gdip_DisposeImage(pBitmap),Gdip_DeleteGraphics(G)Run,%FileOutExit:Gdip_Shutdown(pToken)ExitApp
DIMx(9),y(9)x()=0,1,2,3,4,5,6,7,8,9y()=2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0ORIGIN100,100VDU23,23,2;0;0;0;VDU5FORx=1TO9GCOL7:LINE100*x,720,100*x,0GCOL0:PLOT0,-10,-4:PRINT;x;NEXTFORy=20TO180STEP20GCOL7:LINE900,4*y,0,4*yGCOL0:PLOT0,-212,20:PRINTy;NEXTLINE0,0,0,720LINE0,0,900,0GCOL4FORi%=0TO9IFi%=0THENMOVE100*x(i%),4*y(i%)ELSEDRAW100*x(i%),4*y(i%)ENDIFNEXT
We could use thesuite provided byRaster graphics operations, but those functions lack a facility to draw text.
#include<stdio.h>#include<stdlib.h>#include<math.h>#include<plot.h>#define NP 10doublex[NP]={0,1,2,3,4,5,6,7,8,9};doubley[NP]={2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0};voidminmax(double*x,double*y,double*minx,double*maxx,double*miny,double*maxy,intn){inti;*minx=*maxx=x[0];*miny=*maxy=y[0];for(i=1;i<n;i++){if(x[i]<*minx)*minx=x[i];if(x[i]>*maxx)*maxx=x[i];if(y[i]<*miny)*miny=y[i];if(y[i]>*maxy)*maxy=y[i];}}/* likely we must play with this parameter to make the plot looks better when using different set of data */#define YLAB_HEIGHT_F 0.1#define XLAB_WIDTH_F 0.2#define XDIV (NP*1.0)#define YDIV (NP*1.0)#define EXTRA_W 0.01#define EXTRA_H 0.01#define DOTSCALE (1.0/150.0)#define MAXLABLEN 32#define PUSHSCALE(X,Y) pl_fscale((X),(Y))#define POPSCALE(X,Y) pl_fscale(1.0/(X), 1.0/(Y))#define FMOVESCALE(X,Y) pl_fmove((X)/sx, (Y)/sy)intmain(){intplotter,i;doubleminx,miny,maxx,maxy;doublelx,ly;doublexticstep,yticstep,nx,ny;doublesx,sy;charlabs[MAXLABLEN+1];plotter=pl_newpl("png",NULL,stdout,NULL);if(plotter<0)exit(1);pl_selectpl(plotter);if(pl_openpl()<0)exit(1);/* determines minx, miny, maxx, maxy */minmax(x,y,&minx,&maxx,&miny,&maxy,NP);lx=maxx-minx;ly=maxy-miny;pl_fspace(floor(minx)-XLAB_WIDTH_F*lx,floor(miny)-YLAB_HEIGHT_F*ly,ceil(maxx)+EXTRA_W*lx,ceil(maxy)+EXTRA_H*ly);/* compute x,y-ticstep */xticstep=(ceil(maxx)-floor(minx))/XDIV;yticstep=(ceil(maxy)-floor(miny))/YDIV;pl_flinewidth(0.25);/* compute scale factors to adjust aspect */if(lx<ly){sx=lx/ly;sy=1.0;}else{sx=1.0;sy=ly/lx;}pl_erase();/* a frame... */pl_fbox(floor(minx),floor(miny),ceil(maxx),ceil(maxy));/* labels and "tics" */pl_fontname("HersheySerif");for(ny=floor(miny);ny<ceil(maxy);ny+=yticstep){pl_fline(floor(minx),ny,ceil(maxx),ny);snprintf(labs,MAXLABLEN,"%6.2lf",ny);FMOVESCALE(floor(minx)-XLAB_WIDTH_F*lx,ny);PUSHSCALE(sx,sy);pl_label(labs);POPSCALE(sx,sy);}for(nx=floor(minx);nx<ceil(maxx);nx+=xticstep){pl_fline(nx,floor(miny),nx,ceil(maxy));snprintf(labs,MAXLABLEN,"%6.2lf",nx);FMOVESCALE(nx,floor(miny));PUSHSCALE(sx,sy);pl_ftextangle(-90);pl_alabel('l','b',labs);POPSCALE(sx,sy);}/* plot data "point" */pl_fillcolorname("red");pl_filltype(1);for(i=0;i<NP;i++){pl_fbox(x[i]-lx*DOTSCALE,y[i]-ly*DOTSCALE,x[i]+lx*DOTSCALE,y[i]+ly*DOTSCALE);}pl_flushpl();pl_closepl();}
No one would use the previous code to produce a plot (that looksthis way; instead, normally we produce data through a program, then we plot the data using e.g.gnuplot or other powerful tools; the result (with gnuplot and without enhancement) could looklike this instead.
Following code creates a plot in EPS format, with auto scaling and line/symbol/color controls. Plotting function loosely follows Matlab command style. Not thorough by any means, just to give an idea on how this kind of things can be coded.
#include<stdio.h>#include<math.h>#include<string.h>#define N 40doublex[N],y[N];voidminmax(doublex[],intlen,double*base,double*step,int*nstep){inti;doublediff,minv,maxv;*step=1;minv=maxv=x[0];for(i=1;i<len;i++){if(minv>x[i])minv=x[i];if(maxv<x[i])maxv=x[i];}if(minv==maxv){minv=floor(minv);maxv=ceil(maxv);if(minv==maxv){minv--;maxv++;}}else{diff=maxv-minv;while(*step<diff)*step*=10;while(*step>diff)*step/=10;if(*step>diff/2)*step/=5;elseif(*step>diff/5)*step/=2;}*base=floor(minv/*step)**step;*nstep=ceil(maxv/*step)-floor(minv/*step);}/* writes an eps with 400 x 300 dimention, using 12 pt font */#define CHARH 12#define CHARW 6#define DIMX 398#define DIMY (300 - CHARH)#define BOTY 20.intplot(doublex[],doubley[],intlen,char*spec){intnx,ny,i;doublesx,sy,x0,y0;charbuf[100];intdx,dy,lx,ly;doubleofs_x,ofs_y,grid_x;minmax(x,len,&x0,&sx,&nx);minmax(y,len,&y0,&sy,&ny);dx=-log10(sx);dy=-log10(sy);ly=0;for(i=0;i<=ny;i++){sprintf(buf,"%g\n",y0+i*sy);if(strlen(buf)>ly)ly=strlen(buf);}ofs_x=ly*CHARW;printf("%%!PS-Adobe-3.0\n%%%%BoundingBox: 0 0 400 300\n""/TimesRoman findfont %d scalefont setfont\n""/rl{rlineto}def /l{lineto}def /s{setrgbcolor}def ""/rm{rmoveto}def /m{moveto}def /st{stroke}def\n",CHARH);for(i=0;i<=ny;i++){ofs_y=BOTY+(DIMY-BOTY)/ny*i;printf("0 %g m (%*.*f) show\n",ofs_y-4,ly,dy,y0+i*sy);if(i)printf("%g %g m 7 0 rl st\n",ofs_x,ofs_y);}printf("%g %g m %g %g l st\n",ofs_x,BOTY,ofs_x,ofs_y);for(i=0;i<=nx;i++){sprintf(buf,"%g",x0+i*sx);lx=strlen(buf);grid_x=ofs_x+(DIMX-ofs_x)/nx*i;printf("%g %g m (%s) show\n",grid_x-CHARW*lx/2,BOTY-12,buf);if(i)printf("%g %g m 0 7 rl st\n",grid_x,BOTY);}printf("%g %g m %g %g l st\n",ofs_x,BOTY,grid_x,BOTY);if(strchr(spec,'r'))printf("1 0 0 s\n");elseif(strchr(spec,'b'))printf("0 0 1 s\n");elseif(strchr(spec,'g'))printf("0 1 0 s\n");elseif(strchr(spec,'m'))printf("1 0 1 s\n");if(strchr(spec,'o'))printf("/o { m 0 3 rm 3 -3 rl -3 -3 rl -3 3 rl closepath st} def "".5 setlinewidth\n");if(strchr(spec,'-')){for(i=0;i<len;i++){printf("%g %g %s ",(x[i]-x0)/(sx*nx)*(DIMX-ofs_x)+ofs_x,(y[i]-y0)/(sy*ny)*(DIMY-BOTY)+BOTY,i?"l":"m");}printf("st\n");}if(strchr(spec,'o'))for(i=0;i<len;i++){printf("%g %g o ",(x[i]-x0)/(sx*nx)*(DIMX-ofs_x)+ofs_x,(y[i]-y0)/(sy*ny)*(DIMY-BOTY)+BOTY);}printf("showpage\n%%EOF");return0;}intmain(){inti;for(i=0;i<N;i++){x[i]=(double)i/N*3.14159*6;y[i]=-1337+(exp(x[i]/10)+cos(x[i]))/100;}/* string parts: any of "rgbm": color; "-": draw line; "o": draw symbol */plot(x,y,N,"r-o");return0;}
#include<windows.h>#include<string>#include<vector>//--------------------------------------------------------------------------------------------------usingnamespacestd;//--------------------------------------------------------------------------------------------------constintHSTEP=46,MWID=40,MHEI=471;constfloatVSTEP=2.3f;//--------------------------------------------------------------------------------------------------classvector2{public:vector2(){x=y=0;}vector2(floata,floatb){x=a;y=b;}voidset(floata,floatb){x=a;y=b;}floatx,y;};//--------------------------------------------------------------------------------------------------classmyBitmap{public:myBitmap():pen(NULL),brush(NULL),clr(0),wid(1){}~myBitmap(){DeleteObject(pen);DeleteObject(brush);DeleteDC(hdc);DeleteObject(bmp);}boolcreate(intw,inth){BITMAPINFObi;ZeroMemory(&bi,sizeof(bi));bi.bmiHeader.biSize=sizeof(bi.bmiHeader);bi.bmiHeader.biBitCount=sizeof(DWORD)*8;bi.bmiHeader.biCompression=BI_RGB;bi.bmiHeader.biPlanes=1;bi.bmiHeader.biWidth=w;bi.bmiHeader.biHeight=-h;HDCdc=GetDC(GetConsoleWindow());bmp=CreateDIBSection(dc,&bi,DIB_RGB_COLORS,&pBits,NULL,0);if(!bmp)returnfalse;hdc=CreateCompatibleDC(dc);SelectObject(hdc,bmp);ReleaseDC(GetConsoleWindow(),dc);width=w;height=h;returntrue;}voidclear(BYTEclr=0){memset(pBits,clr,width*height*sizeof(DWORD));}voidsetBrushColor(DWORDbClr){if(brush)DeleteObject(brush);brush=CreateSolidBrush(bClr);SelectObject(hdc,brush);}voidsetPenColor(DWORDc){clr=c;createPen();}voidsetPenWidth(intw){wid=w;createPen();}voidsaveBitmap(stringpath){BITMAPFILEHEADERfileheader;BITMAPINFOinfoheader;BITMAPbitmap;DWORDwb;GetObject(bmp,sizeof(bitmap),&bitmap);DWORD*dwpBits=newDWORD[bitmap.bmWidth*bitmap.bmHeight];ZeroMemory(dwpBits,bitmap.bmWidth*bitmap.bmHeight*sizeof(DWORD));ZeroMemory(&infoheader,sizeof(BITMAPINFO));ZeroMemory(&fileheader,sizeof(BITMAPFILEHEADER));infoheader.bmiHeader.biBitCount=sizeof(DWORD)*8;infoheader.bmiHeader.biCompression=BI_RGB;infoheader.bmiHeader.biPlanes=1;infoheader.bmiHeader.biSize=sizeof(infoheader.bmiHeader);infoheader.bmiHeader.biHeight=bitmap.bmHeight;infoheader.bmiHeader.biWidth=bitmap.bmWidth;infoheader.bmiHeader.biSizeImage=bitmap.bmWidth*bitmap.bmHeight*sizeof(DWORD);fileheader.bfType=0x4D42;fileheader.bfOffBits=sizeof(infoheader.bmiHeader)+sizeof(BITMAPFILEHEADER);fileheader.bfSize=fileheader.bfOffBits+infoheader.bmiHeader.biSizeImage;GetDIBits(hdc,bmp,0,height,(LPVOID)dwpBits,&infoheader,DIB_RGB_COLORS);HANDLEfile=CreateFile(path.c_str(),GENERIC_WRITE,0,NULL,CREATE_ALWAYS,FILE_ATTRIBUTE_NORMAL,NULL);WriteFile(file,&fileheader,sizeof(BITMAPFILEHEADER),&wb,NULL);WriteFile(file,&infoheader.bmiHeader,sizeof(infoheader.bmiHeader),&wb,NULL);WriteFile(file,dwpBits,bitmap.bmWidth*bitmap.bmHeight*4,&wb,NULL);CloseHandle(file);delete[]dwpBits;}HDCgetDC()const{returnhdc;}intgetWidth()const{returnwidth;}intgetHeight()const{returnheight;}private:voidcreatePen(){if(pen)DeleteObject(pen);pen=CreatePen(PS_SOLID,wid,clr);SelectObject(hdc,pen);}HBITMAPbmp;HDChdc;HPENpen;HBRUSHbrush;void*pBits;intwidth,height,wid;DWORDclr;};//--------------------------------------------------------------------------------------------------classplot{public:plot(){bmp.create(512,512);}voiddraw(vector<vector2>*pairs){bmp.clear(0xff);drawGraph(pairs);plotIt(pairs);HDCdc=GetDC(GetConsoleWindow());BitBlt(dc,0,30,512,512,bmp.getDC(),0,0,SRCCOPY);ReleaseDC(GetConsoleWindow(),dc);//bmp.saveBitmap( "f:\\rc\\plot.bmp" );}private:voiddrawGraph(vector<vector2>*pairs){HDCdc=bmp.getDC();bmp.setPenColor(RGB(240,240,240));DWORDb=11,c=40,x;RECTrc;chartxt[8];for(x=0;x<pairs->size();x++){MoveToEx(dc,40,b,NULL);LineTo(dc,500,b);MoveToEx(dc,c,11,NULL);LineTo(dc,c,471);wsprintf(txt,"%d",(pairs->size()-x)*20);SetRect(&rc,0,b-9,36,b+11);DrawText(dc,txt,lstrlen(txt),&rc,DT_RIGHT|DT_VCENTER|DT_SINGLELINE);wsprintf(txt,"%d",x);SetRect(&rc,c-8,472,c+8,492);DrawText(dc,txt,lstrlen(txt),&rc,DT_CENTER|DT_VCENTER|DT_SINGLELINE);c+=46;b+=46;}SetRect(&rc,0,b-9,36,b+11);DrawText(dc,"0",1,&rc,DT_RIGHT|DT_VCENTER|DT_SINGLELINE);bmp.setPenColor(0);bmp.setPenWidth(3);MoveToEx(dc,40,11,NULL);LineTo(dc,40,471);MoveToEx(dc,40,471,NULL);LineTo(dc,500,471);}voidplotIt(vector<vector2>*pairs){HDCdc=bmp.getDC();HBRUSHbr=CreateSolidBrush(255);RECTrc;bmp.setPenColor(255);bmp.setPenWidth(2);vector<vector2>::iteratorit=pairs->begin();inta=MWID+HSTEP*static_cast<int>((*it).x),b=MHEI-static_cast<int>(VSTEP*(*it).y);MoveToEx(dc,a,b,NULL);SetRect(&rc,a-3,b-3,a+3,b+3);FillRect(dc,&rc,br);it++;for(;it<pairs->end();it++){a=MWID+HSTEP*static_cast<int>((*it).x);b=MHEI-static_cast<int>(VSTEP*(*it).y);SetRect(&rc,a-3,b-3,a+3,b+3);FillRect(dc,&rc,br);LineTo(dc,a,b);}DeleteObject(br);}myBitmapbmp;};//--------------------------------------------------------------------------------------------------intmain(intargc,char*argv[]){ShowWindow(GetConsoleWindow(),SW_MAXIMIZE);plotpt;vector<vector2>pairs;pairs.push_back(vector2(0,2.7f));pairs.push_back(vector2(1,2.8f));pairs.push_back(vector2(2.0f,31.4f));pairs.push_back(vector2(3.0f,38.1f));pairs.push_back(vector2(4.0f,58.0f));pairs.push_back(vector2(5.0f,76.2f));pairs.push_back(vector2(6.0f,100.5f));pairs.push_back(vector2(7.0f,130.0f));pairs.push_back(vector2(8.0f,149.3f));pairs.push_back(vector2(9.0f,180.0f));pt.draw(&pairs);system("pause");return0;}//--------------------------------------------------------------------------------------------------
(use'(incantercorestatscharts))(defx(range010))(defy'(2.72.831.438.158.076.2100.5130.0149.3180.0))(view(xy-plotxy))
Boost.Process is part ofDelphiBoostLib.
programPlot_coordinate_pairs;{$APPTYPE CONSOLE}usesSystem.SysUtils,Boost.Process;varx:TArray<Integer>;y:TArray<Double>;beginx:=[0,1,2,3,4,5,6,7,8,9];y:=[2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0];varplot:=TPipe.Create('gnuplot -p',True);plot.WriteA('unset key; plot''-'''#10);forvari:=0toHigh(x)doplot.WriteA(format('%d %f'#10,[x[i],y[i]]));plot.writeA('e'#10);writeln('Press enter to close');Readln;plot.Kill;plot.Free;readln;end.
DuckDB interfaces nicely with programs likeYouPlot that can accept tabular data on their standard input.The YouPlot options used here are:
The execution pipeline could be setup as follows:
![[1]](/mt/?noimg=&dark=on&url=http%3a%2f%2fi.imgur.com%2f0RQaxNl.png){kind=link}
duckdb -c " set variable x = [0, 1, 2, 3, 4, 5, 6, 7, 8, 9]; set variable y = [2.7, 2.8, 31.4, 38.1, 58.0, 76.2, 100.5, 130.0, 149.3, 180.0]; COPY (select unnest(getvariable('x')) as x, unnest(getvariable('y')) as y) TO '/dev/stdout' WITH (FORMAT 'csv', DELIM '\t', HEADER);" | uplot scatter -H -t "Co-ordinate Pairs"x[] = [ 0 1 2 3 4 5 6 7 8 9 ]y[] = [ 2.7 2.8 31.4 38.1 58.0 76.2 100.5 130.0 149.3 180.0 ]# glinewidth 0.5gline 10 5 10 97gline 10 5 95 5gtextsize 3n = len x[]for i to n : max = higher y[i] maxglinewidth 0.1sty = max div 9stx = x[n] div 9for i range0 10 gline 10 5 + i * 10 95 5 + i * 10 gtext 2 4 + i * 10 i * sty gline i * 9 + 10 5 i * 9 + 10 97 gtext i * 9 + 10 1 i * stx.gcolor 900glinewidth 0.5gpenupfor i = 1 to n x = x[i] * 9 / stx + 10 y = y[i] / sty * 10 + 5 glineto x y.
Resulting imagehere.
(lib'plot)(defineys#(2.72.831.438.158.076.2100.5130.0149.3180.0))(define(fn)[ysn])(plot-sequencef9)→(("x:auto"09)("y:auto"2198))(plot-grid120)(plot-text" Rosetta plot coordinate pairs"010"white")
UsingEplot to produce PNG.
-module(plot_coordinate_pairs).-export([task/0,to_png_file/3]).task()->Xs=[0,1,2,3,4,5,6,7,8,9],Ys=[2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0],File="plot_coordinate_pairs",to_png_file(File,Xs,Ys).to_png_file(File,Xs,Ys)->PNG=egd_chart:graph([{File,lists:zip(Xs,Ys)}]),file:write_file(File++".png",PNG).
The result looks likethis.
Using theF# for Visualization library:
#r@"C:\Program Files\FlyingFrog\FSharpForVisualization.dll"letx=Seq.mapfloat[|0;1;2;3;4;5;6;7;8;9|]lety=[|2.7;2.8;31.4;38.1;58.0;76.2;100.5;130.0;149.3;180.0|]openFlyingFrog.GraphicsPlot([Data(Seq.zipxy)],(0.0,9.0))
USING:accessorsassocscolors.constantskernelsequencesuiui.gadgetsui.gadgets.chartsui.gadgets.charts.lines;chartnew{{0 9}{0 180}}>>axeslinenewCOLOR:blue>>color9<iota>{2.7 2.8 31.4 38.1 58 76.2 100.5 130 149.3 180}zip>>dataadd-gadget"Coordinate pairs"open-window
Fōrmulæ programs are not textual, visualization/edition of programs is done showing/manipulating structures but not text. Moreover, there can be multiple visual representations of the same program. Even though it is possible to have textual representation —i.e. XML, JSON— they are intended for storage and transfer purposes more than visualization and edition.
Programs in Fōrmulæ are created/edited online in itswebsite.
Inthis page you can see and run the program(s) related to this task and their results. You can also change either the programs or the parameters they are called with, for experimentation, but remember that these programs were created with the main purpose of showing a clear solution of the task, and they generally lack any kind of validation.
Solution
DimAsIntegeri,x(9),y(9)Fori=0To9x(i)=iNextiy(0)=2.7y(1)=2.8y(2)=31.4y(3)=38.1y(4)=58.0y(5)=76.2y(6)=100.5y(7)=130.0y(8)=149.3y(9)=180.0Locate22,4Fori=0To9Locate22,((i*4)+2):PrintiNextiFori=0To20Step2Locate(21-i),0:Print(i*10)NextiColor14Fori=0To9Locate(21-(y(i)/10)),(x(i)*4)+2:Print"."NextiSleep
unsetkey# Only one data set, so the key is uninformativeplot'-'# '-' can be replaced with a filename, to read data from that file.02.712.8231.4338.1468.0576.26100.57130.08149.39180.0e
Output is the same as for the gnuplot example on this page.
A program can of course supply commands and data to gnuplot as prepared files. For the spirit of controlling plotting with the native language however, this example shows how commands and data can be prepared programatically and supplied to gnuplot through stdin.
packagemainimport("fmt""log""os/exec")var(x=[]int{0,1,2,3,4,5,6,7,8,9}y=[]float64{2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0})funcmain(){g:=exec.Command("gnuplot","-persist")w,err:=g.StdinPipe()iferr!=nil{log.Fatal(err)}iferr=g.Start();err!=nil{log.Fatal(err)}fmt.Fprintln(w,"unset key; plot '-'")fori,xi:=rangex{fmt.Fprintf(w,"%d %f\n",xi,y[i])}fmt.Fprintln(w,"e")w.Close()g.Wait()}
packagemainimport("log""github.com/gonum/plot""github.com/gonum/plot/plotter""github.com/gonum/plot/plotutil""github.com/gonum/plot/vg")var(x=[]int{0,1,2,3,4,5,6,7,8,9}y=[]float64{2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0})funcmain(){pts:=make(plotter.XYs,len(x))fori,xi:=rangex{pts[i]=struct{X,Yfloat64}{float64(xi),y[i]}}p,err:=plot.New()iferr!=nil{log.Fatal(err)}iferr=plotutil.AddScatters(p,pts);err!=nil{log.Fatal(err)}iferr:=p.Save(3*vg.Inch,3*vg.Inch,"points.svg");err!=nil{log.Fatal(err)}}
Using JFreeChart and Groovy Swing Builder
importgroovy.swing.SwingBuilderimportjavax.swing.JFrameimportorg.jfree.chart.ChartFactoryimportorg.jfree.chart.ChartPanelimportorg.jfree.data.xy.XYSeriesimportorg.jfree.data.xy.XYSeriesCollectionimportorg.jfree.chart.plot.PlotOrientationdefchart={x=[0,1,2,3,4,5,6,7,8,9]y=[2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0]defseries=newXYSeries('plots')[x,y].transpose().each{x,y->series.addx,y}deflabels=["Plot Demo","X","Y"]defdata=newXYSeriesCollection(series)defoptions=[false,true,false]defchart=ChartFactory.createXYLineChart(*labels,data,PlotOrientation.VERTICAL,*options)newChartPanel(chart)}newSwingBuilder().edt{frame(title:'Plot coordinate pairs',defaultCloseOperation:JFrame.EXIT_ON_CLOSE,pack:true,show:true){widget(chart())}}
gnuplot is a package fromHackageDB.
importGraphics.Gnuplot.Simplepnts=[2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0]doPlot=plotPathStyle[(Title"plotting dots")](PlotStylePoints(CustomStyle[]))(zip[0..]pnts)
REAL :: n=10, x(n), y(n)x = (0, 1, 2, 3, 4, 5, 6, 7, 8, 9)y = (2.7, 2.8, 31.4, 38.1, 58.0, 76.2, 100.5, 130.0, 149.3, 180.0)WINDOW(WINdowhandle=wh, Width=-300, Height=-300, X=1, TItle='Rosetta')AXIS(WINdowhandle=wh, Title='x values', Yaxis, Title='y values')LINE(X=x, Y=y, SymbolDiameter=2)
linkprintf,numbersproceduremain()x:=[0.,1.,2.,3.,4.,5.,6.,7.,8.,9.]y:=[2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0]Plot(x,y,600,400)end$definePOINTR2# Point Radius$definePOINTC"red"# Point Colour$defineGRIDC"grey"# grid colour$defineAXISC"black"# axis/label colour$defineBORDER60# per side border$defineTICKS5.# grid ticks per axis$defineAXISFH20# font height for axis labelsprocedurePlot(x,y,cw,ch)/cw:=700# default dimensions/ch:=400uw:=cw-BORDER*2# usable dimensionsuh:=ch-BORDER*2wparms:=["Plot","g",sprintf("size=%d,%d",cw,ch),"bg=white"]# base window parmsdx:=sprintf("dx=%d",BORDER)# grid origindy:=sprintf("dy=%d",BORDER)&window:=open!wparms|stop("Unable to open window")X:=scale(x,uw)# scale data to usable spaceY:=scale(y,uh,"invert")WAttrib(dx,dy)# set origin=grid & draw grideveryx:=(X.tickfromtoX.ticktobyX.tick)*X.tickscaledo{ifx=0thenFg(AXISC)elseFg(GRIDC)DrawLine(x,Y.tickfrom*Y.tickscale,x,Y.tickto*Y.tickscale)}everyy:=(Y.tickfromtoY.ticktobyY.tick)*Y.tickscaledo{ify=uhthenFg(AXISC)elseFg(GRIDC)DrawLine(X.tickfrom*X.tickscale,y,X.tickto*X.tickscale,y)}Fg(POINTC)# draw data points ....everyi:=1to*X.scaleddoFillCircle(X.scaled[i],Y.scaled[i],POINTR)Fg(AXISC)# label gridWAttrib(dx,"dy=0")# label X axisFont(sprintf("Helvetica,%d",AXISFH))ytxt:=ch-BORDER+1+(WAttrib("ascent")-WAttrib("descent"))/2everyx:=X.tickscale*(xv:=X.tickfromtoX.ticktobyX.tick)doDrawString(x-TextWidth(xv)/2,ytxt+integer(AXISFH*1.5),xv)WAttrib("dx=0",dy)# label Y axiseveryy:=Y.tickscale*(yv:=Y.tickfromtoY.ticktobyY.tick)doDrawString(BORDER/2-TextWidth(yv)/2,ytxt-BORDER-y,yv)WriteImage(sprintf("PlotPoints-%d.gif",&now))# save imageWAttrib("dx=0","dy=0")# close off nicelyFont("Helvetica,10")DrawString(10,ch-5,"Right click to exit")untilEvent()==&rpress# wait for left mouse buttonclose(&window)endrecordscaledata(low,high,range,pix,raw,scaled,tick,tickfrom,tickto,tickscale)procedurescale(data,pix,opts[])P:=scaledata(pmin:=min!data,pmax:=max!data,prange:=real(pmax-pmin),pix,data,q:=[])/ticks:=TICKSP.tick:=ceil(prange/(10^(k:=floor(log(prange,10))))*(10^k)/ticks)P.tickfrom:=P.tick*floor(pmin/P.tick)P.tickto:=P.tick*ceil(pmax/P.tick)P.tickscale:=real(pix)/(P.tickto-P.tickfrom)everyput(q,integer((!data-P.tickfrom)*P.tickscale))if!opts=="invert"then# invert is for yeveryq[i:=1to*q]:=pix-q[i]returnPend
printf.icn provides formattingnumbers.icn provides floor,ceil
require'plot'X=:i.10Y=:2.72.831.438.158.076.2100.5130.0149.3180.0'dot; pensize 2.4'plotX;Y
If you eliminate the left argument to plot, the dots will be connected instead of being isolated.
importjava.awt.*;importjava.awt.event.*;importjava.awt.geom.*;importjavax.swing.JApplet;importjavax.swing.JFrame;publicclassPlot2dextendsJApplet{double[]xi;double[]yi;publicPlot2d(double[]x,double[]y){this.xi=x;this.yi=y;}publicstaticdoublemax(double[]t){doublemaximum=t[0];for(inti=1;i<t.length;i++){if(t[i]>maximum){maximum=t[i];}}returnmaximum;}publicstaticdoublemin(double[]t){doubleminimum=t[0];for(inti=1;i<t.length;i++){if(t[i]<minimum){minimum=t[i];}}returnminimum;}publicvoidinit(){setBackground(Color.white);setForeground(Color.white);}publicvoidpaint(Graphicsg){Graphics2Dg2=(Graphics2D)g;g2.setRenderingHint(RenderingHints.KEY_ANTIALIASING,RenderingHints.VALUE_ANTIALIAS_ON);g2.setPaint(Color.black);intx0=70;inty0=10;intxm=670;intym=410;intxspan=xm-x0;intyspan=ym-y0;doublexmax=max(xi);doublexmin=min(xi);doubleymax=max(yi);doubleymin=min(yi);g2.draw(newLine2D.Double(x0,ym,xm,ym));g2.draw(newLine2D.Double(x0,ym,x0,y0));for(intj=0;j<5;j++){intinterv=4;g2.drawString(""+(j*(xmax-xmin)/interv+xmin),j*xspan/interv+x0-10,ym+20);g2.drawString(""+(j*(ymax-ymin)/interv+ymin),x0-20-(int)(9*Math.log10(ymax)),ym-j*yspan/interv+y0-5);g2.draw(newLine2D.Double(j*xspan/interv+x0,ym,j*xspan/interv+x0,ym+5));g2.draw(newLine2D.Double(x0-5,j*yspan/interv+y0,x0,j*yspan/interv+y0));}for(inti=0;i<xi.length;i++){intf=(int)((xi[i]-xmin)*xspan/(xmax-xmin));inth=(int)(((ymax-ymin)-(yi[i]-ymin))*yspan/(ymax-ymin));g2.drawString("o",x0+f-3,h+14);}for(inti=0;i<xi.length-1;i++){intf=(int)((xi[i]-xmin)*xspan/(xmax-xmin));intf2=(int)((xi[i+1]-xmin)*xspan/(xmax-xmin));inth=(int)(((ymax-ymin)-(yi[i]-ymin))*yspan/(ymax-ymin));inth2=(int)(((ymax-ymin)-(yi[i+1]-ymin))*yspan/(ymax-ymin));g2.draw(newLine2D.Double(f+x0,h+y0,f2+x0,h2+y0));}}publicstaticvoidmain(Stringargs[]){JFramef=newJFrame("ShapesDemo2D");f.addWindowListener(newWindowAdapter(){publicvoidwindowClosing(WindowEvente){System.exit(0);}});double[]r={0,1,2,3,4,5,6,7,8,9};double[]t={2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.09};JAppletapplet=newPlot2d(r,t);f.getContentPane().add("Center",applet);applet.init();f.pack();f.setSize(newDimension(720,480));f.show();}}
Using R (non-interactive mode)
jq is designed to interoperate with other tools, and in this section we illustrate how jq can be used with R in a simple pipeline: jq will produce a stream of CSV data that will be piped into R operating in non-interactive mode. Assuming the jq and R programs are respectively in plot.jq and plot.R, the pipeline would look like this:
jq-n-M-r-fplot.jq|RCMDBATCHplot.RThe above would produce the plot as a .pdf file.
plot.jq
# NOTE: This definition of transpose can be omitted# if your version of jq includes transpose as a builtin.## transpose a possibly jagged matrix, quickly; # rows are padded with nulls so the result is always rectangular.def transpose: if . == [] then [] else . as $in | (map(length) | max) as $max | length as $length | reduce range(0; $max) as $j ([]; . + [reduce range(0;$length) as $i ([]; . + [ $in[$i][$j] ] )] ) end;def x: [0, 1, 2, 3, 4, 5, 6, 7, 8, 9];def y: [2.7, 2.8, 31.4, 38.1, 58.0, 76.2, 100.5, 130.0, 149.3, 180.0];def plot(x;y): "A,B", ( [x,y] | transpose | map( @csv ) | .[]);plot(x;y)
plot.R
mydata<-read.table(file("stdin"),header=TRUE,sep=",")x=mydata$A# x-axisy=mydata$B# y-axisplot(x,y,# plot the variablesmain="Scatterplot Example",xlab="x-axis label",# x-axis labelylab="y-axis label")# y-axis label
Using Plots library with PlotlyJS as backend:
usingPlotsplotlyjs()x=[0,1,2,3,4,5,6,7,8,9]y=[2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0]p=scatter(x,y)savefig(p,"/tmp/testplot.png")
// Version 1.2.31importorg.jfree.chart.ChartFactoryimportorg.jfree.chart.ChartPanelimportorg.jfree.data.xy.XYSeriesimportorg.jfree.data.xy.XYSeriesCollectionimportorg.jfree.chart.plot.PlotOrientationimportjavax.swing.JFrameimportjavax.swing.SwingUtilitiesimportjava.awt.BorderLayoutfunmain(args:Array<String>){valx=intArrayOf(0,1,2,3,4,5,6,7,8,9)valy=doubleArrayOf(2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0)valseries=XYSeries("plots")(0untilx.size).forEach{series.add(x[it],y[it])}vallabels=arrayOf("Plot Demo","X","Y")valdata=XYSeriesCollection(series)valoptions=booleanArrayOf(false,true,false)valorient=PlotOrientation.VERTICALvalchart=ChartFactory.createXYLineChart(labels[0],labels[1],labels[2],data,orient,options[0],options[1],options[2])valchartPanel=ChartPanel(chart)SwingUtilities.invokeLater{valf=JFrame()with(f){defaultCloseOperation=JFrame.EXIT_ON_CLOSEadd(chartPanel,BorderLayout.CENTER)title="Plot coordinate pairs"isResizable=falsepack()setLocationRelativeTo(null)isVisible=true}}}
Similar to Groovy entry.
1)defineX&Y:{defX012345789}->X{defY2.72.831.438.158.076.2100.5130.0149.3180.0}->Y2)defineafunctionreturningasequenceofSVGpoints{defcurve{lambda{:curve:kx:ky}{S.map{{lambda{:curve:kx:ky:i}{*:kx{S.get:i{{car:curve}}}}{*:ky{S.get:i{{cdr:curve}}}}}:curve:kx:ky}{S.serie0{-{S.length{X}}1}}}}}3)drawapolylineinaSVGcontext{svg{@width="580"height="300"style="background:#eee"}{g{AXES580300}{polyline{@points="{curve {cons X Y} 30 0.9}"stroke="#000"fill="transparent"stroke-width="1"}}}}where{defAXES{lambda{:w:h}{@transform="translate({/ :w 2},{/ :h 2}) scale(1,-1)"}{line{@x1="-{/ :w 2}:w"y1="0"x2="{/ :w 2}"y2="0"stroke="red"fill="transparent"}}{line{@x1="0"y1="-{/ :h 2}"x2="0"y2="{/ :h 2}"stroke="green"fill="transparent"}}}}4)theresultcanbeseeninhttp://lambdaway.free.fr/lambdawalks/?view=plot4
First version writes directly to LB's console window.
'Plotting coordinate pairs MainWin - StyleFor i = 0 To 9 x(i) = iNext iy(0) = 2.7y(1) = 2.8y(2) = 31.4y(3) = 38.1y(4) = 58.0y(5) = 76.2y(6) = 100.5y(7) = 130.0y(8) = 149.3y(9) = 180.0Locate 4, 22For i = 0 To 9 Locate ((i * 4) + 2), 22 Print iNext iFor i = 0 To 20 Step 2 Locate 0, (21 - i) Print (i * 10)Next iFor i = 0 To 9 Locate (x(i) * 4) + 2, (21 - (y(i)/ 10)) Print "."Next iEnd
The second version uses the more typical graphic window approach, and is written to enable easy adaptation to other data sets.
nomainwinDATA 0, 1, 2, 3, 4, 5, 6, 7, 8, 9DATA 2.7, 2.8, 31.4, 38.1, 58.0, 76.2, 100.5, 130.0, 149.3, 180.0For i = 0 To 9 READ tmp: x( i) = tmpNext iFor i = 0 To 9 READ tmp: y( i) = tmpNext i'Plotting coordinate pairsWindowHeight = 500WindowWidth = 430Open "Plot coordinate pairs" For Graphics_nsb_nf As #gwin#gwin "trapclose [quit.gwin]"#gwin "Color Black; Down"'25, 418 is 0,0global offsetX, offsetY, scaleX, scaleYoffsetX = 25: offsetY = 418scaleX = 40: scaleY = 2maxX = 9: maxY = 200#gwin "line "; sx( maxX);" "; sy( 0);" "; sx( 0); " "; sy( 0)#gwin "goto "; sx( 0); " "; sy( maxY)For x = 0 To 9 #gwin "place ";sx(x);" ";sy(0) #gwin "Go -18" #gwin "|"; xNext#gwin "turn 90"For y = 0 To 200 Step 20 #gwin "place ";sx(0);" ";sy(y) #gwin "Go -5" #gwin "place ";0;" ";sy(y) #gwin "|"; yNext#gwin "size 3"For i = 0 To 9 #gwin "set ";sx(x(i));" ";sy(y(i))Next i#gwin "Flush"Wait[quit.gwin] Close #gwin End'x, y to screen x, yfunction sx(x) sx = offsetX +x*scaleXend functionfunction sy(y) sy = offsetY-y*scaleY 'y is invertedend function
Displaying the plot with vector graphics
on plotGraphic local tCoordinates local x = "0, 1, 2, 3, 4, 5, 6, 7, 8, 9" local y = "2.7, 2.8, 31.4, 38.1, 58.0, 76.2, 100.5, 130.0, 149.3, 180.0" if there is a graphic "graph" then delete graphic "graph" repeat with i = 1 to the number of items of x put item i of x into item 1 of line i of tCoordinates put item i of y into item 2 of line i of tCoordinates end repeat create graphic "graph" set the style of graphic "graph" to "polygon" set the points of graphic "graph" to tCoordinates repeat with i = 1 to the number of lines of tCoordinates put the top of grc "graph" + the height of grc "graph" - item 2 of line i of tCoordinates into item 2 of line i of tCoordinates end repeat set the points of graphic "graph" to tCoordinates set the height of graphic "graph" to 200 set the width of graphic "graph" to 300 set the loc of grc "graph" to the loc of this cardend plotGraphic
Displaying the plot with the Line Graph widget
on plotLineGraph local tCoordinates local x = "0, 1, 2, 3, 4, 5, 6, 7, 8, 9" local y = "2.7, 2.8, 31.4, 38.1, 58.0, 76.2, 100.5, 130.0, 149.3, 180.0" if there is a widget "graph" then delete widget "graph" repeat with i = 1 to the number of items of x put item i of x into item 1 of line i of tCoordinates put item i of y into item 2 of line i of tCoordinates end repeat create widget "graph" as "com.livecode.widget.linegraph" set the graphData of widget "graph" to tCoordinates set the height of widget "graph" to 250 set the width of widget "graph" to 350 set the loc of widget "graph" to the loc of this cardend plotLineGraph
w_width=love.graphics.getWidth()w_height=love.graphics.getHeight()x={0,1,2,3,4,5,6,7,8,9}y={2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0}origin={24,24}points={}x_unit=w_width/x[10]/2y_unit=w_height/10--add points to an array properly formatted for the line functionfori=1,10,1dotable.insert(points,(x[i]*x_unit)+origin[1])table.insert(points,(w_height-(y[i]*2))-origin[2])endfunctionlove.draw()--draw axes and gridlove.graphics.setColor(0,0.8,0)--draw x axislove.graphics.line(origin[1],w_height-origin[2],w_width,w_height-origin[2])--draw y axislove.graphics.line(origin[1],w_height-origin[2],origin[1],origin[2])--draw gridfori=1,20,1dolove.graphics.line(origin[1],(w_height-origin[2])-(i*y_unit),w_width,(w_height-origin[2])-(i*y_unit))love.graphics.line(origin[1]+(i*x_unit),origin[2],origin[1]+(i*x_unit),w_height-origin[2])end--draw line plotlove.graphics.setColor(0.8,0,0)love.graphics.line(points)--draw labelslove.graphics.setColor(0.8,0.8,0.8)fori=0,9,1do--draw x axis labelslove.graphics.print(i,(x_unit*i)+origin[1],love.graphics.getHeight()-origin[2])--draw y axis labelslove.graphics.print(i*y_unit/2,origin[1],((love.graphics.getHeight()-i*y_unit)-origin[2]))endend
Last statements used for copy the console screen to clipboard
Result imagehere
Module Pairs { \\ written in version 9.5 rev. 13 \\ use Gdi+ antialiasing (not work with Wine in Linux, but we get no error) smooth on Const center=2, right=3, left=1, blue=1, angle=0, dotline=3 Const size9pt=9, size11pt=11 Cls ,0 ' use current background color, set split screen from line 0 Cursor 0,3 Report center, "Coordinate pairs" x = (0, 1, 2, 3, 4, 5, 6, 7, 8, 9) y = (2.7, 2.8, 31.4, 38.1, 58.0, 76.2, 100.5, 130.0, 149.3, 180.0) dx=scale.x/2/len(x) dy=dx 'ratio 1:1 graH=dy*len(x) Basex=scale.x/4 Basey=(scale.y+graH)/2 Move Basex, Basey \\ draw use relative coordinates Draw 0,-graH \\ Step just move graphic cursor Step 0, graH Draw scale.x/2 Step -scale.x/2 \\ scX is 1, not used max=Y#max() \\ Auto scale for Y, using 0 for start of axis Y scY=-graH/((max+5^log(max) ) div 100)/100 \\ make vertical axis using dots with numbers center per dx j=1 For i=basex+dx to basex+dx*x#max() Step dx Move i, basey Step 0, twipsy*10 Legend format$("{0}",array(x,j)), "courier", size9pt, angle, center Width 1, dotline { draw 0, -graH-twipsy*10,7} j++ Next i \\ the same for horizontal axis HalfTextHeight=Size.y("1","courier", size9pt)/2 For i=basey-dy to basey-dy*x#max() Step dy Move basex, i Step -twipsx*10 Width 1, dotline { draw scale.x/2+twipsx*10,,7} Move basex-100, i+HalfTextHeight Legend format$("{0}",(i-basey)/scY), "courier", size9pt, angle, left Next i ex=each(x) : ey=each(y) \\ start from first point. We use Draw to for absolute coordinates Move array(x,0)*dx+Basex, array(y,0)*scy+Basey While ex, ey { Width 2 { Draw to array(ex)*dx+Basex, array(ey)*scy+Basey, blue } } \\ second pass for marks and labels ex=each(x) : ey=each(y) While ex, ey { Move array(ex)*dx+Basex, array(ey)*scy+Basey Step -75, -75 Pen 12 {draw 150: draw 0,150 : draw -150 : draw 0,-150} Pen 13 { Step 200, -200 Legend format$("({0}-{1})",array(ex),array(ey) ), "courier bold", size11pt, angle, right } } \\ screenshot to clipboard Screenshot$="" Move 0,0 Copy scale.x, scale.y to Screenshot$ Clipboard Screenshot$ a$=key$}Pairsx:=Vector([0,1,2,3,4,5,6,7,8,9]):y:=Vector([2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0]):plot(x,y,style="point");
x={0,1,2,3,4,5,6,7,8,9};y={2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0};ListPlot[{x,y}//Transpose]
>>x=[0,1,2,3,4,5,6,7,8,9];>>y=[2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0];>>plot(x,y,'.-')
".."(m,n):=makelist(i,i,m,n);infix("..")$x:0..9$y:[2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0]$wxplot2d(['discrete,x,y],[style,[points,5,1,1]],[gnuplot_term,png],[gnuplot_out_file,"qsort-range-10-9.png"])$
There exists two libraries providing a Nim interface to “gnuplot” (the GNU plotting program), which are named “gnuplot” and “gnuplotlib”. We have chosen the second one as it is be more complete and more convenient to use.
The library launches “gnuplot” which does the plotting. From “gnuplot”, it is possible to save the drawing into a PDF, a SVG or an image (BMP, PNG) file.
importgnuplotletx=[0,1,2,3,4,5,6,7,8,9]y=[2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0]withGnuPlot:plot(x,y,"Coordinate pairs")
This library doesn’t use an external process to does the plotting. It uses a syntax mostly compliant with “ggplot2” syntax.
importggplotnimletx=[0,1,2,3,4,5,6,7,8,9]y=[2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0]letdf=seqsToDf(x,y)# Build a dataframe.df.ggplot(aes("x","y"))+ggtitle("Coordinate pairs")+geomLine()+themeOpaque()+ggsave("coordinate_pairs.png")
#load"graphics.cma"openGraphicsletroundx=int_of_float(floor(x+.0.5))letx=[0;1;2;3;4;5;6;7;8;9]andy=[2.7;2.8;31.4;38.1;58.0;76.2;100.5;130.0;149.3;180.0]let()=open_graph"";List.iter2(funxy->(* scale to fit in the window *)let_x=x*30and_y=round(y*.2.0)inplot_x_y)xy;ignore(wait_next_event[Key_pressed]);close_graph();;;
Using theArchimedes library, one can write:
moduleA=Archimedeslety=[|2.7;2.8;31.4;38.1;58.0;76.2;100.5;130.0;149.3;180.0|]let()=letvp=A.init[]inA.Axes.boxvp;A.set_colorvpA.Color.red;A.Array.yvpy;A.closevp
x=[0,1,2,3,4,5,6,7,8,9];y=[2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0];plot(x,y,"o");pause;
; define input arrays(definex'(0123456789))(definey'(2.72.831.438.158.076.2100.5130.0149.3180.0)); render(import(libgl2))(glOrtho010020001)(gl:set-renderer(lambda(mouse)(glClearGL_COLOR_BUFFER_BIT)(glColor3f010)(glBeginGL_LINE_STRIP)(mapglVertex2fxy)(glEnd)))
plothraw(vx, vy)
usesPlotWPF;beginvarx:=|0.0,1,2,3,4,5,6,7,8,9|;vary:=|2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0|;varplot:=newMarkerGraphWPF(x,y);plot.addlinegraph(x,y);end.
useGD::Graph::points;@data=([0,1,2,3,4,5,6,7,8,9],[2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0],);$graph=GD::Graph::points->new(400,300);openmy$fh,'>',"qsort-range-10-9.png";binmode$fh;print$fh$graph->plot(\@data)->png;close$fh;
useImager;useImager::Plot;@x=(0,1,2,3,4,5,6,7,8,9);@y=(2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0);$plot=Imager::Plot->new(Width=>400,Height=>300,GlobalFont=>'PATH_TO_TTF_FONT',);$plot->AddDataSet(X=>\@x,Y=>\@y,style=>{marker=>{size=>2,symbol=>'circle',color=>Imager::Color->new('red'),},},);$img=Imager->new(xsize=>500,ysize=>400,);$img->box(filled=>1,color=>'white');$plot->Render(Image=>$img,Xoff=>50,Yoff=>350);$img->write(file=>'qsort-range-10-9.png');
//Aamrun, 26th June 2022intx[]={0,1,2,3,4,5,6,7,8,9};floaty[]={2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0};size(300,300);surface.setTitle("Rosetta Plot");stroke(#ff0000);for(inti=0;i<x.length;i++){ellipse(x[i],y[i],3,3);}
Output same as BBC BASIC, you can run this onlinehere.
---- demo\rosetta\Plot_coordinate_pairs.exw-- ======================================--withjavascript_semanticsincludepGUI.eincludeIupGraph.econstantx={0,1,2,3,4,5,6,7,8,9},y={2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0}functionget_data(Ihandlegraph)return{{x,y,CD_BLUE}}endfunctionIupOpen()Ihandlegraph=IupGraph(get_data),dlg=IupDialog(graph,`TITLE="Plot coordinate pairs"`)IupSetAttributes(dlg,"RASTERSIZE=320x240,MINSIZE=320x200")IupSetAttributes(graph,"XTICK=1,XMIN=0,XMAX=9")IupSetAttributes(graph,"YTICK=20,YMIN=0,YMAX=180")IupShow(dlg)ifplatform()!=JSthenIupMainLoop()IupClose()endif
(load "@lib/ps.l")(scl 1)(de plot (PsFile DX DY Lst) (let (SX (length Lst) SY (apply max Lst) N 0 Val) (out PsFile (psHead (+ DX 20) (+ DY 40)) (font (9 . "Helvetica")) (if (or (=0 SX) (=0 SY)) (window 60 12 DX DY (font 24 ,"Not enough Data") ) (setq Lst # Build coordinates (let X -1 (mapcar '((Y) (cons (*/ (inc 'X) DX SX) (- DY (*/ Y DY SY)) ) ) Lst ) ) ) (color 55 95 55 # Background color (let (X (+ DX 40) Y (+ DY 40)) (poly T 0 0 X 0 X Y 0 Y 0 0) ) ) (window 20 20 DX DY # Plot coordinates (poly NIL 0 0 0 DY (- DX 20) DY) (color 76 24 24 (poly NIL (caar Lst) (cdar Lst) (cdr Lst)) ) ) (window 4 4 60 12 (ps (format SY *Scl))) (for X SX (window (+ 6 (*/ (dec X) DX SX)) (+ 24 DY) 30 12 (ps (format (dec X)) 0) ) ) ) (page) ) ) )(plot "plot.ps" 300 200 (2.7 2.8 31.4 38.1 58.0 76.2 100.5 130.0 149.3 180.0))(call 'display "plot.ps")
/x[0123456789]def/y[2.72.831.438.158.076.2100.5130.0149.3180.0]def/i1defnewpathx0gety0getmovetoxlength1sub{xigetyigetlineto/ii1adddef}repeatstroke
StructurePlotDatax.iy.fEndStructureGlobali,x,y.f,max_x,max_y,min_x=#MAXLONG,min_y=Infinity()Definecount=(?serie_y-?serie_x)/SizeOf(Integer)-1GlobalDimMyData.PlotData(count)Restoreserie_xFori=0TocountRead.ixMyData(i)\x=xIfx>max_x:max_x=x:EndIfIfx<min_x:min_x=x:EndIfNextRestoreserie_yFori=0TocountRead.fyMyData(i)\y=yIfy>max_y:max_y=y:EndIfIfy<min_y:min_y=y:EndIfNextProcedureUpdatePlot(Win,w,h)Staticgblm=20,gtrm=5;graph's bottom-left and top-right marginProtectedcount=ArraySize(MyData())Ifw>gblmAndh>gblmAndcount>0SetWindowTitle(Win,"PureBasic Plot "+Str(w)+"x"+Str(h))Protectedgw=w-gblm,gh=h-gblm;graph's width and heightProtectedi,yf.f,xf.fyf=(gh-gtrm)/max_yxf=(gw-gtrm)/max_xCreateImage(0,w,h)ProtectedOutputID=ImageOutput(0)StartDrawing(OutputID)DrawingMode(#PB_2DDrawing_Transparent);-DrawgridFori=0Tocounty=gh-max_y*i/count*yfLineXY(gblm,y,w-gtrm,y,$467E3E);Y-scaleDrawText(1,y-5,RSet(StrD(i/count*max_y,1),5))x=gblm+max_x*i/count*xfy=gh;X-ScaleLineXY(x,y,x,gtrm,$467E3E)Ifi:DrawText(x-5,y+2,Str(i)):EndIfNext;-DrawcurveProtectedox=gblm,oy=gh,x,yFori=0Tocountx=gblm+MyData(i)\x*xfy=gh-MyData(i)\y*yfLineXY(ox,oy,x,y,$0133EE)ox=x:oy=yNextStopDrawing()ImageGadget(0,0,0,w,h,ImageID(0))EndIfEndProcedureDefineWin=OpenWindow(#PB_Any,0,0,600,400,"",#PB_Window_SystemMenu|#PB_Window_SizeGadget)IfWinSmartWindowRefresh(Win,1)UpdatePlot(Win,WindowWidth(Win),WindowHeight(Win))RepeatDefineevent=WaitWindowEvent()SelecteventCase#PB_Event_SizeWindowUpdatePlot(Win,WindowWidth(Win),WindowHeight(Win))EndSelectUntilevent=#PB_Event_CloseWindow;SavetheplotiftheuserwantstoIfMessageRequester("Question","Save it?",#PB_MessageRequester_YesNo)=#PB_MessageRequester_YesDefineFile$=SaveFileRequester("Save as","PB.png","PNG (*.png)|*.png",0)UsePNGImageEncoder()SaveImage(0,File$,#PB_ImagePlugin_PNG)EndIfEndIfDataSectionserie_x:Data.i0,1,2,3,4,5,6,7,8,9serie_y:Data.f2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0EndDataSection
Interactive session:
>>>x=[0,1,2,3,4,5,6,7,8,9]>>>y=[2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0]>>>importpylab>>>pylab.plot(x,y,'bo')>>>pylab.savefig('qsort-range-10-9.png')
See some other examples:
fromvisualimport*fromvisual.graphimport*plot1=gdisplay(title='VPython Plot-Demo',xtitle='x',ytitle='y (click and drag mouse to see coordinates)',foreground=color.black,background=color.white,x=0,y=0,width=400,height=400,xmin=0,xmax=10,ymin=0,ymax=200)f1=gdots(color=color.red)# create plot-objectf1.plot(pos=(0,2.7),color=color.blue)# add a single pointf1.plot(pos=[(1,2.8),# add a list of points(2,31.4),(3,38.1),(4,58.0),(5,76.2),(6,100.5),(7,130.0),(8,149.3),(9,180.0)])label(display=plot1.display,text="Look here",pos=(6,100.5),xoffset=30,yoffset=-20)
R has several different plotting paradigms. First we define the data.
x<-c(0,1,2,3,4,5,6,7,8,9)y<-c(2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0)
plot(x,y)
library(lattice)xyplot(y~x)
library(ggplot2)qplot(x,y)
Racket has a built-in plotting library
#langracket(requireplot)(definex(build-list10values))(definey(list2.72.831.438.158.076.2100.5130.0149.3180.0))(plot-new-window?#t)(plot(points(mapvectorxy)))
This opens a new window with this image (with interactive zooming)
And this
#langracket(requireplot)(definex(build-list10values))(definey(list2.72.831.438.158.076.2100.5130.0149.3180.0))(plot-new-window?#t)(plot(lines(mapvectorxy)))
opens a new window with this image
(formerly Perl 6)
Generate an SVG image file.
useSVG;useSVG::Plot;my@x =0..9;my@y = (2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0);saySVG.serialize:SVG::Plot.new(width =>512,height =>512,x =>@x,x-tick-step => {1 },min-y-axis =>0,values => [@y,],title =>'Coordinate Pairs',).plot(:lines);
See Plot coordinate pairs/REXX for the$PLOT program.
Example usage:
/*REXX program plots X,Y coördinate pairs of numbers with plain (ASCII) characters.*/x=0123456789y=2.72.831.438.158.076.2100.5130.0149.3180.0$=doj=1forwords(x)/*build a list suitable for $PLOT subr.*/$=$word(x,j)','word(y,j)/*add this X,Y coördinate to the $ list*/end/*j*//*$≡ 0,2.7 1,2.8 2,31.4 3,38.1 ··· */call'$PLOT'$/*invoke the REXX program: $PLOT */exitrc/*stick a fork in it, we're all done. */
│180 ∙││││││││ ∙│││││ ∙│││││││ ∙│││││││ ∙││││ ∙│││││ ∙││ ∙││││││∙0 ∙ 9└──────────────────────────────────────────────────────────────────────────────
/*REXX program plots X,Y coördinate pairs of numbers with plain (ASCII) characters.*/x=0123456789y=2.72.831.438.158.076.2100.5130.0149.3180.0$=doj=1forwords(x)/*build a list suitable for $PLOT subr.*/$=$word(x,j)','word(y,j)/*add this X,Y coördinate to the $ list*/end/*j*//*$≡ 0,2.7 1,2.8 2,31.4 3,38.1 ··· */call'$PLOT'$'(LABELDatapoints'/*invoke the REXX program: $PLOT */exitrc/*stick a fork in it, we're all done. */
│180 (9,180)∙││││││││││ ∙(8,149.3)││││││ ∙(7,130)││││││││││ ∙(6,100.5)│││││││ ∙(5,76.2)││││││ ∙(4,58)│││││││ ∙(3,38.1)││ ∙(2,31.4)││││││││∙(0,2.7) ∙(1,2.8)└──────────────────────────────────────────────────────────────────────────────────────────────
![[2]](/mt/?noimg=&dark=on&url=http%3a%2f%2fi43.tinypic.com%2f2a689yw.png){kind=link}
# Project : Plot coordinate pairsload "guilib.ring"paint = nullnew qapp { win1 = new qwidget() { setwindowtitle("Plot coordinate pairs") setgeometry(100,100,1024,900) label1 = new qlabel(win1) { setgeometry(10,10,1024,900) settext("") } new qpushbutton(win1) { setgeometry(50,50,100,30) settext("draw") setclickevent("draw()") } show() } exec() }func draw p1 = new qpicture() color = new qcolor() { setrgb(0,0,255,255) } pen = new qpen() { setcolor(color) setwidth(1) } paint = new qpainter() { begin(p1) setpen(pen) old = 0 yold = 0 xnew = 0 ynew = 0 x2 = [0, 1, 2, 3, 4, 5, 6, 7, 8, 9] y2 = [2.7, 2.8, 31.4, 38.1, 58.0, 76.2, 100.5, 130.0, 149.3, 180.0] for x = 1 to 9 drawline(100*x,720,100*x,0) drawtext(100*x,750,string(x)) next for y = 20 to 180 step 20 drawline(900,4*y,0,4*y) drawtext(0,720-4*y,string(y)) next drawline(0,0,0,720) drawline(0,0,900,0) for i = 1 to 10 if i=1 xold = 100*x2[i] yold = 720-4*y2[i] else xnew = 100*x2[i] ynew = 720-4*y2[i] drawline(xold,yold,xnew,ynew) xold = xnew yold = ynew ok next endpaint() } label1 { setpicture(p1) show() } returnOutput:
https://www.dropbox.com/s/q6tra0cqoty4pya/Plot.jpg?dl=0
≪ → x y ≪ ERASE@ clear graphics display x 0 + ≪ MIN ≫ STREAM x ≪ MAX ≫ STREAM XRNG@ set x range y 0 + ≪ MIN ≫ STREAM y ≪ MAX ≫ STREAM YRNG@ set y range 1 x SIZEFOR j x j GET y j GET R→C@ generate coordinates pairIF j 1 >THEN SWAP OVER LINEEND@ draw a line from previous pairNEXT DROP { (0,0) 10 "x" "y" } AXES DRAX LABEL@ draw axes and labels { } PVIEW@ display graphics≫ ≫ 'PLOTXY' STO
{0 1 2 3 4 5 6 7 8 9} {2.7 2.8 31.4 38.1 58.0 76.2 100.5 130.0 149.3 180.0}PLOTXY
require'gnuplot'x=[0,1,2,3,4,5,6,7,8,9]y=[2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0]Gnuplot.opendo|gp|Gnuplot::Plot.new(gp)do|plot|plot.data<<Gnuplot::DataSet.new([x,y])do|ds|ds.with="linespoints"ds.notitleendendend
importscala.swing.Swing.pair2Dimensionimportscala.swing.{MainFrame,Panel,Rectangle}importjava.awt.{Color,Graphics2D,geom}objectPlotCoordPairsextendsscala.swing.SimpleSwingApplication{//min/max of display-x resp. yval(dx0,dy0)=(70,30)val(dxm,dym)=(670,430)val(prefSizeX,prefSizeY)=(720,480)privatedefui=newPanel{importmath._valxmax={valf1=pow(10,log10(xs.max).toInt)valf2=if(f1<10)10elseround(xs.max/f1)*f1if(f2>=xs.max)f2else(round(xs.max/f1)+1)*f1}valymax={valf1=pow(10,log10(ys.max).toInt)valf2=if(f1<10)10elseround(ys.max/f1)*f1if(f2>=ys.max)f2else(round(ys.max/f1)+1)*f1}val(xinterv,yinterv)=(xmax/xs.size,ymax/xs.size)caseclassCoord(x:Double,y:Double){val(dx,dy)=((x/xmax*(dxm-dx0)+dx0).toInt,(dym-y/ymax*(dym-dy0)).toInt)}valpcentre=Coord(0,0)valpxmax=Coord(xmax,0)valpymax=Coord(0,ymax)background=Color.whitepreferredSize=(prefSizeX,prefSizeY)//axes:vala_path=newgeom.GeneralPatha_path.moveTo(pxmax.dx,pxmax.dy)a_path.lineTo(pcentre.dx,pcentre.dy)//x-axisa_path.lineTo(pymax.dx,pymax.dy)//y-axis// interval ticks:xs.map(i=>Coord(i*xinterv,0)).map(p=>{a_path.moveTo(p.dx,p.dy)a_path.lineTo(p.dx,p.dy+5)})xs.map(i=>Coord(0,i*yinterv)).map(p=>{a_path.moveTo(p.dx,p.dy)a_path.lineTo(p.dx-5,p.dy)})//grid:valg_path=newgeom.GeneralPath(1toxs.size).map(i=>Coord(i*xinterv,0)).map(p=>{g_path.moveTo(p.dx,p.dy);g_path.lineTo(Coord(p.x,ymax).dx,Coord(p.x,ymax).dy)})(1toxs.size).map(i=>Coord(0,i*yinterv)).map(p=>{g_path.moveTo(p.dx,p.dy);g_path.lineTo(Coord(xmax,p.y).dx,Coord(xmax,p.y).dy)})//labeling:valxlabels=(0toxs.size).map(i=>{valp=Coord(i*xinterv,0)Triple(p.x.toInt.toString,p.dx-3,p.dy+20)})valylabels=(0toxs.size).map(i=>{valp=Coord(0,i*yinterv)Triple(p.y.toInt.toString,p.dx-30,p.dy+5)})//curve:valpath=newgeom.GeneralPathvalcurve=xs.map(i=>Coord(xs(i),ys(i)))path.moveTo(curve.head.dx,curve.head.dy)curve.map(p=>path.lineTo(p.dx,p.dy))//...flag all function values:valrects=curve.map(p=>newRectangle(p.dx-3,p.dy-3,6,6))overridedefpaintComponent(g:Graphics2D)={super.paintComponent(g)g.setColor(Color.lightGray)g.draw(g_path)g.setColor(Color.black)g.draw(a_path)xlabels.map(t=>g.drawString(t._1,t._2,t._3))ylabels.map(t=>g.drawString(t._1,t._2,t._3))g.draw(path)rects.map(g.draw(_))}}valxs=0to9valys:List[Double]=List(2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0)deftop=newMainFrame{title="Rosetta Code >>> Task: Plot coordinate pairs | Language: Scala"contents=ui}}
-->x=[0,1,2,3,4,5,6,7,8,9];-->y=[2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0];-->plot2d(x,y)
require('GD::Graph::points')vardata=[[0,1,2,3,4,5,6,7,8,9],[2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0],]vargraph=%s'GD::Graph::points'.new(400,300)vargd=graph.plot(data)varformat='png'File("qsort-range.#{format}").write(gd.(format),:raw)
clearinput x y02.712.8231.4338.1458.0576.26100.57130.08149.39180.0endlines y xgraph export image.png
This solution does not use existing plotting packages, but constructs the graphics from bare-metal Tk code.
packagerequireTk# The actual plotting engineprocplotxy{canvasxsys}{globalxfacyfacsetmaxx[tcl::mathfunc::max{*}$xs]setmaxy[tcl::mathfunc::max{*}$ys]setxfac[expr{[winfowidth$canvas]*0.8/$maxx}]setyfac[expr{[winfoheight$canvas]*0.8/$maxy}]scale$canvasx0$maxx$xfacscale$canvasy0$maxy$yfacforeachx$xsy$ys{dot$canvas[expr{$x*$xfac}][expr{$y*$yfac}]-fillred}}# Rescales the contents of the given canvasprocscale{canvasdirectionfromtofac}{setf[expr{$from*$fac}]sett[expr{$to*$fac}]switch--$direction{x{setf[expr{$from*$fac}]sett[expr{$to*$fac}]$canvascreateline$f0$t0$canvascreatetext$f0-anchornw-text$from$canvascreatetext$t0-anchorn-text$to}y{setf[expr{$from*-$fac}]sett[expr{$to*-$fac}]$canvascreateline0$f0$t$canvascreatetext0$f-anchorse-text$from$canvascreatetext0$t-anchore-text$to}}}# Helper to make points, which are otherwise not a native item typeprocdot{canvasxyargs}{setid[$canvascreateoval[expr{$x-3}][expr{-$y-3}]\[expr{$x+3}][expr{-$y+3}]]$canvasitemconfigure$id{*}$args}pack[canvas.c-backgroundwhite]updatesetxs{0123456789}setys{2.72.831.438.158.076.2100.5130.0149.3180.0}plotxy.c$xs$ys.cconfig-scrollregion[.cbboxall].cmoveall2020# Save image (this is the only part that requires an external library)packagerequireImgsetim[imagecreatephoto-data.c]$imwriteplotxy.png-formatPNG
Of course, if we were generating an encapsulated postscript version, we would be able to do that directly.
Note also that in Tk 8.6, there is no need for an external library to write PNG images; the capability is directly supported.
FnOffPlotsOffNewPlot 1, 1, x, yZoomData
Ursala doesn't plot anything directly, but has libraries to generate LaTeX code for 2D and 3D graphics. The output file has to be run through LaTeX or included into a LaTeX document. Here's the way to do it just as a quick check (all default settings and dots connected with straight lines).
#import std#import flo#import fit#import plox = <0., 1., 2., 3., 4., 5., 6., 7., 8., 9.>y = <2.7, 2.8, 31.4, 38.1, 58.0, 76.2, 100.5, 130.0, 149.3, 180.0>#output dot'tex' latex_document+ plotmain = visualization[curves: <curve[points: ~&p/x y]>]
(output)
Here's one way you might do it if you were interested in publicationquality graphics. The dots are connected with a cubic spline interpolatingfunction sampled at 200 points, and the axes are nicelylabeled.
main =visualization[ abscissa: axis[ variable: 'problem size', hats: printf/*'%0.0f' ari10/0. 9.], ordinates: ~&iNC axis[ variable: 'execution time ($\mu$s)', hats: printf/*'%0.1f' ari6/0. 180.], curves: < curve[ points: ^(~&,chord_fit0@p/x y)* ari200/0. 9., attributes: {'linecolor': 'lightgray'}], curve[ scattered: true, points: ~&p/x y, attributes: {'linecolor': 'black'}]>](output)
Using Excel
PrivateSubplot_coordinate_pairs(xAsVariant,yAsVariant)DimchrtAsChartSetchrt=ActiveSheet.Shapes.AddChart.ChartWithchrt.ChartType=xlLine.HasLegend=False.HasTitle=True.ChartTitle.Text="Time".SeriesCollection.NewSeries.SeriesCollection.Item(1).XValues=x.SeriesCollection.Item(1).Values=y.Axes(xlValue,xlPrimary).HasTitle=True.Axes(xlValue,xlPrimary).AxisTitle.Characters.Text="microseconds"EndWithEndSubPublicSubmain()x=[{0,1,2,3,4,5,6,7,8,9}]y=[{2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0}]plot_coordinate_pairsx,yEndSub
import"graphics"forCanvas,ImageData,Colorimport"dome"forWindowimport"math"forPointclassPlotCoordinates{constructnew(width,height){Window.title="Plot coordinates"Window.resize(width,height)Canvas.resize(width,height)_width=width_height=height}init(){varx=[0,1,2,3,4,5,6,7,8,9]vary=[2.7,2.8,31.4,38.1,58.0,76.2,100.5,130.0,149.3,180.0]plotCoordinates(x,y)}plotCoordinates(x,y){varn=x.count// draw axesCanvas.line(40,_height-40,_width-40,_height-40,Color.blue,2)Canvas.line(40,_height-40,40,40,Color.blue,2)varlength=40*nvardiv=length/10varj=0for(iin0..9){varp=Point.new(40+j,_height-40)Canvas.print(i.toString,p.x-4,p.y+4,Color.white)j=j+div}j=divfor(iin1..9){varp=Point.new(10,_height-40-j)vars=(i*20).toStringif(s.count==2)s=" "+sCanvas.print(s,p.x,p.y,Color.white)j=j+div}Canvas.print("X",_width-44,_height-36,Color.green)Canvas.print("Y",20,40,Color.green)// plot pointsvarxStart=40varxScale=40varyStart=40varyScale=2varpoints=List.filled(n,null)for(iin0...n){points[i]=Point.new(xStart+x[i]*xScale,_height-yStart-y[i]*yScale)}Canvas.circlefill(points[0].x,points[0].y,3,Color.red)for(iin1...n){Canvas.line(points[i-1].x,points[i-1].y,points[i].x,points[i].y,Color.red)Canvas.circlefill(points[i].x,points[i].y,3,Color.red)}}update(){}draw(alpha){}}varGame=PlotCoordinates.new(500,500)
XPL0 does not provide a library routine for plotting graphs. An issuewith this particular task is how general to make the code. This is prettyspecific.
include c:\cxpl\codes; \intrinsic 'code' declarationsdef ScrW=640, ScrH=480, VidMode=$101;def Sx = ScrW/10, \pixels per horz grid line Sy = ScrH/10, \pixels per vert grid line Ox = (3+1+1)*8+2, \offset for horz grid: allow room for "180.0" Oy = ScrH-20; \offset for vert grid: allow room for labelsint X, DataX;real Y, DataY, Gain;def Brown=6, LCyan=11;[DataX:= [0, 1, 2, 3, 4, 5, 6, 7, 8, 9]; DataY:= [2.7, 2.8, 31.4, 38.1, 58.0, 76.2, 100.5, 130.0, 149.3, 180.0];SetVid(VidMode);for X:= 0 to 9 do \draw grid [Move(X*Sx+Ox, Oy); Line(X*Sx+Ox, Oy-9*Sy, Brown); \vert lines Move(Ox, Oy-X*Sy); Line(9*Sx+Ox, Oy-X*Sy, Brown); \horz lines ];Format(3,1); Attrib(LCyan); \label gridY:= 0.0;for X:= 0 to 9 do [Move(X*Sx+Ox-3, Oy+6); IntOut(6, X); \X axis Move(0, Oy-X*Sy-7); RlOut(6, Y); \Y axis Y:= Y + 20.0; ];Gain:= float(Sy)/20.0;Move(DataX(0)*Sx+Ox, Oy-Fix(DataY(0)*Gain)); \plot pointsfor X:= 1 to 9 do Line(DataX(X)*Sx+Ox, Oy-Fix(DataY(X)*Gain), LCyan);if ChIn(1) then []; \wait for keySetVid(3); \restore text]
x = [0, 1, 2, 3, 4, 5, 6, 7, 8, 9];y = [2.7, 2.8, 31.4, 38.1, 58.0, 76.2, 100.5, 130.0, 149.3, 180.0];window, 0;plmk, y, x;window, 1;plg, y, x, marks=0;
Solution using gnuplot. Output is the same as for the gnuplot example on this page.
A program can of course supply commands and data to gnuplot as prepared files. For the spirit of controlling plotting with the native language however, this example shows how commands and data can be prepared programmatically and supplied to gnuplot through stdin.
#<<<cmd:=0'|#set term wxt # X11unset key # Only one data set, so the key is uninformative plot '-' # '-' can be replaced with a filename, to read data from that file. 0 2.7 1 2.8 2 31.4 3 38.1 4 68.0 5 76.2 6 100.5 7 130.0 8 149.3 9 180.0e|;#<<<gnuplot:=System.popen("gnuplot","w");gnuplot.write(cmd); gnuplot.flush();ask("Hit return to finish"); gnuplot.close();