procedure Input_Data; begin tinc:= N; tinc:= 1.0/tinc; for i:= 0 to N-1 do begin D1[i]:= exp(-tinc*i)*sin(2*tinc*i+i/7.5)*sqr(cos(4*tinc*i)); D2[i]:= 0.0; end; end;