%Test of solve_inhomogeneous using step function
dt = 5e-8
N = 100

E = f_step(N); % Applied EMF
[V,t] = solve_inhomogeneous(@D_lowpass, E, dt); % Solution V(t)

Va = ifftn(anti_alias(fftn(V)));

figure(1)
hold off
plot(t,E,'k-','linewidth',2)
hold on
plot(t,V,'b-','linewidth',2)
plot(t,Va,'g-','linewidth',2)
legend('E(t)','V(t)','V(t) anti-aliased')
print('test3.pdf','-dpdf')