function [y]=init(x) y=cos(4*%pi*x)-cos(6*%pi*x);