% Set up initial data for fish4
  nfish = 10; d=20; n=10;
  disp('Initialize fish for simulation')
  nfish = input('Input the number of fish to be simulated(10 is reasonable): ');
  rand('seed',0);
  fisho = 10*(rand(nfish,1) + sqrt(-1)*rand(nfish,1));
  d = input('Input the size of the ocean (number of grid points on a side)(20 is reasonable): ');
  n = input('Input the number of time steps to simulate(10 is reasonable): ');