static public static void figure( int  figno)

static public static void subplot( int  p)

static public void title(String titleStr)

static public void xlabel(String xLabelStr)

static public void ylabel(String yLabelStr)

static public void zlabel(String zLabelStr) 

static public PlotPanel plot(double [] x, double [] y, double [] z, int low, int high, Color color, String

 static public PlotPanel splinePlot( double [] xvals, double [] yvals, int NP, Color color, String name)
 
      

val N=3000
val axis = linspace(0.0, 10.0, N)
val x1 = 2.3*sin(4.5*axis)
val x2 = -4.5*cos(axis*7.8)+67.8*sin(16.7*axis)
val x3 = x1-7.8*x2+0.7*sin(3.4*axis)

figure(1)
subplot(3,1,1); plot(axis, x1); title("First signal x1")
subplot(3,1,2); plot(axis, x2); title("Second signal x2")
subplot(3,1,3); plot(axis, x3); title("Third signal x3")

var N=4000

var A = rand(N,N)    // create a large random matrix in 2-D storage representation, i.e. as Java's 2-D arrays

tic()   // start timer
var AA = A*A  // multiply using native OpenBLAS
var tm=toc()
println("performing native multiplication in $tm secs")

var N=4000

var A = rand1d(N,N)    // create a large random matrix in 1-D row-major storage representation
tic()   // start timer
var AA = A*A  // multiply using native OpenBLAS
var tm=toc()
println("performing native multiplication in $tm secs")

 fun size(M: Mat1D): IntArray  
 fun fill1d(nrows: Int, ncols: Int, value: Double): Mat1D 
 fun rand1d(nrows: Int, ncols: Int): Mat1D 
 fun ones1d(nrows: Int, ncols: Int): Mat1D 
 fun zeros1d(nrows: Int, ncols: Int): Mat1D 

public static Matrix read(File f)
public static Matrix read(String filename) 
public static void write(File f, Matrix M)
public static void write(String filename,Matrix M) 
public  static void mapf(Matrix M, computeFunction c)
public static void resize(Matrix M, int m, int n)
public static Matrix find(Matrix M)

 public  static Matrix sum(Matrix M)
 
 public static Matrix sin(Matrix M)  // similarly all trigonometric i.e. cos, tan, atan etc
 public  static Matrix rand(int m,int  n)
 
public static Matrix inc(double begin, double pitch, double end)
public static Matrix  linspace(double startv,  double endv, int nP)
public static Matrix  logspace(double  startOrig,  double endOrig,  int nP)

 

jfigure(1)

var t = inc(0.0, 0.01, 10.0)


var x = sin(t)

var lineSpecs="."

jplot(t,x, lineSpecs)

println(t*8.8)

jfigure(2)

jsubplot(222)

var x11 = sin(t*4.56)

jplot(t,x11)

jhold(true)

lineSpecs = ":g"

jplot(t,sin(5.0 * x11 ), lineSpecs)

jsubplot(223)

lineSpecs = ":r"

jplot(t,x11, lineSpecs)



// create a new figure and preform a plot at subplot 3,2,1

  jfigure(3)
  
 jsubplot(3,2,1)
 
 var t2 = inc(0.0, 0.01, 10.0);  
 
 var x2 = sin(4.56 *t2 )+cos(t2 * 8.7)
 
 jplot(t2,x2, ".-")
 
jsubplot(3,2,3)

var x3 = cos(t2 * 7.8)+sin(8.9 * t2) * 6.7

jplot(t2, x3)

jlineColor(1, Color.RED)

jsubplot(3,2,5)

var x4 = cos(t2 * 7.66)+sin(3.23 *t2 )*8.665

jplot(t2, x4+x3)

jlineColor(1, Color.GREEN)

jsubplot(3,2,6)

jplot(t2, x4 * 8.77+x3)

jtitle("x4*8.77+x3")

jlineColor(1, Color.BLUE)

// the Ikeda map

var c1 = 0.4

var c2 = 0.9

var c3 = 9.0

var rho = 0.85

var x = 0.5

var y = 0.5

var niters = 200000

var xall = DoubleArray(niters)

var yall = DoubleArray(niters)

tic()

for ( k in 0..niters-1 )  {

  var xp = x; var yp=y
  
  var taut = c1-c3/(1.0 + xp * xp+yp * yp)
  
  x = rho + c2 * xp * cos(taut)-yp * sin(taut)
  
  y = c2*(xp * sin(taut)+yp * cos(taut))
  
  xall[k] = x
  
  yall[k] = y
  
}

var tm = toc()

scatterPlotsOn()

figure(1)

plot(xall, yall, "time = "+tm)