Hi all! In my work I have to define the 5th-derivative of the Gaussian pulse and I do that as follow:
----------------------------------
5TH-DERIVATIVE OF GAUSSIAN PULSE
(define 2PIsqr (sqrt (* 2 pi)))
(define ((GM_5 sigma on-time) t1)
(begin
(let ((t (- t1 (/ on_time 2))))
(if (and (<= t1 on_time) (>= t1 0))
(* 10 (expt t 3)) (* 2PIsqr (expt sigma 9)))
(/ (expt t 5) (* 2PIsqr (expt sigma 11)))
(/ (* 15 t) (* 2PIsqr (expt sigma 7))))
(exp(/ (- 0 (expt t 2)) (* 2 (expt sigma 2)))))
(+ 0 0))
)))
------------------------------------
I know sigma=51ps and I want this pulse be centered in
0.5 ns, so being my time unit 3.33ps :
-----------------------------------
(define w 15.3)
(set! sources (list (make source
(src (make custom-src (src-func (GM_5 w 300))
(end-time 300)))
(component Ey)
(center 0 Cs_y 0))))
-----------------------------------
I try using the function GM_5 and it run correctly,
but now that I complete my example it doesn't run
correctly no more....where am I going wrong?
I enclose my .ctl file.
Thank you in advance!!
Federica
___________________________________
L'email della prossima generazione? Puoi averla con la nuova Yahoo! Mail:
http://it.docs.yahoo.com/nowyoucan.html
planar_elliptical.ctl
Description: 2551654273-planar_elliptical.ctl
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
