function [x,y] = oscdata(n)
% Initialization of oscillatory data
x = linspace(0,n*pi/2,n+1);
y = sin(x);  