for i = 1:5e6:15e6
    fg.write_unsafe(['apply:sin ', num2str(i)])
    fg.opc
    fg.write_unsafe('func:shape ramp')
    fg.opc
end
fg.error

%%
% plot y=11*exp(-|x/6|)*sin(x) from x= -12.5 to x=9.4
x = linspace(-10*pi,10*pi,8000);
y = 11.*exp(-abs(x./6)).*sin(x);
plot(x,y,x+20*pi,y)

%%
tic
fg.frequency = 2000;
fg.waveform = 'sinus';
fg.unit = 'Vrms';
fg.frequency = 1e6;
toc

%%



a=-10*log(1+(2*pi*f*r*c)^2);
s=step(2*pi*f-1/(r*c))*(-20*log(2*pi*f*r*c));
p=-atan(2*pi*f*r*c);

%%
f1=100;
phase_diff1 = 0.5;
tp=((1:1000)./1000);
tp2=tp-((phase_diff1*1000/(f1*360))./1000); 
% normalizing - 360 is 1 complete phase shift

x=0.8*exp(2*pi*i*signal_freq*tp)+0.4*exp(2*pi*i*20*tp);
y=0.8*exp(2*pi*i*signal_freq*tp2)+0.4*exp(2*pi*i*20*tp);
%%
Fs = (wave1.stop-wave1.start)/wave1.length;
p1=(fft(real(wave1.waveform),1024));
p2=(fft(real(wave2.waveform),1024));

phase_diff=(angle(p1(f1)))-(angle(p2(f1)));


%%
tic

[data,raw] = transferFunction(sc,fg,500,5000,50,10);
figure;
subplot(2,1,1)
semilogx(data.frequency,data.magnitude)
subplot(2,1,2)
semilogx(data.frequency,data.phase);
toc