; Parameters:
l_s := 1060 nm  { center wavelength }
T_range := 25 ps  { width of temporal range }
N_t := 2^10  { number of grid points }
dt := T_range / N_t  { temporal resolution }
E0 := 10 nJ  { initial pulse energy }
tau0 := 1 ps  { initial pulse duration }
chirp0 := 0 GHz / ps  { initial chirp }
 
calc  { evaluate the following composite expression }
  begin
  set_pulse_grid(T_range, N_t, l_s);  { define the pulse grid }
  startpulse_G(E0, tau0, chirp0);  { generate the start pulse }
  end

calc  { define the fibers }
  begin
    { active fiber: }
    set_device(1);
    set_fiber(L_active, 20, 'Yb');
    add_ring(r_co, N_Yb);
    pump := addinputchannel(P_pump_in, l_p, 'I_p', 0, forward);
    signal_active_fw := addinputchannel(0, l_s, 'I_s', 0, forward);
    signal_active_bw := addinputchannel(0, l_s, 'I_s', 0, backward);
    set_GVD(GVD_af);
    set_n2(n2_f);  { nonlinear index }
    finish_fiber();
    { passive fiber: }
    set_device(2);
    set_fiber(L_passive, 5, '-');
    signal_passive_fw := addinputchannel(0, l_s, 'I_s', 0, forward);
    signal_passive_bw := addinputchannel(0, l_s, 'I_s', 0, backward);
    set_GVD(GVD_pf);
    set_n2(n2_f);  { nonlinear index }
    finish_fiber();
  end;

DoResonatorRoundTrip() := 
  { Simulate one resonator round-trip of the pulse. }
  begin
    global allow all;
    pp_multiply_expr_f('r_FBG%(f)');
      { fiber Bragg grating as output coupler }
    pp_fiber(1, signal_active_fw);  { active fiber }
    pp_fiber(2, signal_passive_fw);  { passive fiber }
    pp_sat_abs(dR_S, tau_S, E_sat_S);  { SESAM }
    pp_loss(loss_S);  { loss in SESAM }
    pp_fiber(2, signal_passive_bw);  { passive fiber }
    pp_fiber(1, signal_active_bw);  { active fiber }
    pp_center(1);  { center the pulse }
    calc_dyn(0, T_rt, T_rt);
      { dynamical simulation of gain recovery }
  end

請注意，我們的示例基于線性諧振器，在到達另一端（SESAM，可飽和吸收器）之后，我們一直回到初始位置。當然，我們函數中使用的各種參數必須在其他地方定義。

最后一個函數調用在一個往返時間內模擬光纖的泵送，以補充存儲的能量。

此后，很容易模擬一些往返行程并在這些之后顯示腔內脈沖的一些參數：

在這里，我們還存儲了所有脈沖，以便以后可以調用它們進行其他計算和生成圖表，還可以使用交互式脈沖顯示窗口檢查它們。

這里顯示的簡單方法實際上存在一個問題：通常需要進行大量往返操作，直到光纖的增益穩定在最終水平。本質上，存在弛豫振蕩只有經過大量往返之后，這些信息才會消失。如果您想學習這些，那就很好了。但是，如果要快速找到穩定狀態，則可以應用一些其他技巧：使用幾行代碼（此處未顯示），您可以在每次往返過程中根據脈沖能量重新計算光纖的增益-只需好像增益可以立即做出響應。這樣，您可以有效地抑制張弛振蕩，并更快地獲得最終結果。另一種可能性是人為地將增益飽和度大幅度提高（例如100），同時將有效脈沖重復頻率降低相同因子。這將非常有效地加速弛豫振蕩，從而極大地節省了計算時間。當然，

另一個問題是您可能不確定要達到穩定狀態需要往返多少次。因此，您可能需要定義一個模擬往返的函數，直到基本的脈沖屬性（能量，持續時間，帶寬）不再實質改變為止：

"Variation of Pump Power"
 
x: 400, 600
"pump power (mW)", @x
y: 0, 20
frame
 
! for x := CS_x1 to CS_x2 step 20 do
  begin
    set_device(1);  { active fiber }
    set_P_in(pump, x * mW);
    if FindSteadyState(0.01, 1000) then
    begin
      point(x + i * E_p() / nJ, "R");
        { filled rectangle for pulse energy in nJ }
      point(x + i * tau_p() / ps, "t");
        { open triangle for pulse duration in ps }
    end;
  end

• 考慮一個由種子激光器的連續脈沖序列饋送的放大器系統。同樣，您可以模擬脈沖能量如何向穩態發展。您還可以研究，例如，如果您在脈沖之間有隨機的時間間隔，或者只是缺少一些脈沖，或者是一系列脈沖串，或者泵浦功率或種子脈沖的特性發生波動等，那么它是如何波動的。

• 甲再生放大器工作使得一個有時注入種子脈沖，然后讓它在演變一定數量往返，然后噴射它。在抽運放大器晶體一段時間后，可以執行下一個放大周期。在這里，可以輕松實現一種功能，該功能可以在給定數量的諧振器往返行程下模擬完整的放大周期。同樣，可以使用一種函數來模擬這種放大循環，直到達到穩態為止。有些系統會出現分叉甚至混亂的情況，請參閱一些在線示例案例。

• pp_dispersion（）用于應用色散       高達4階
• pp_multiply_expr_f（）用于將脈沖的頻域幅度乘以與頻率相關的因子，
  例如，用于施加任意色散（相移）或帶通濾波器；可以在時域中完成類似
  的操作，例如對于光調制器
• pp_prism_pair（）和pp_grating_pair（）用于分別使用給定參數
  （例如間距，每毫米的行數等）應用棱鏡對或光柵對的色散。
• pp_compress（）用于自動優化的色散壓縮器（最高4階）
• pp_noise（）用于添加一些隨機噪聲
• pp_SPM（）用于具有自相位調制的光學元件
• pp_add_pulse（）用于相干地添加脈沖，例如在干涉儀的模擬中或附加
  脈沖模式鎖定中