clc
clear

plot(sin(x))
t=linspace(0,1e-5,10000);

Harm = 50;

A=0.3;
T=1/10^6;
tr=10*10^-9;
tau=T-2*tr;
k=1/(3*10^7);
m=5*10^-7;
wo=2*pi/T;

a=(k*tr^2)/(2);
b=A*((tr+tau)-tr);
c=(((tr+tau)-(T/2))*(k*((tr-tau)+(T/2)-2*m)))/2;

a0=(1/T)*(a+b+c);

f = a0/2;

for n = 1:2:Harm
    a1 = (k*sin(n*wo*tr))/(n*wo);
    b1 = 0;
    c1 = 0;
    f = f + 2/T*(a1+b1+c1);
end

plot(t,f)