使用S_TIDE工具包里的s_nodal_cal.m函數來計算分潮的交點因子和訂正角，調用程序如下：

[ff,uu]=s_nodal_cal({'K1';'O1';'M2'},[1920,01,01],[2020,01,01],365,10);

subplot(2,1,1);

plot(ff(1,:),'k');hold on;plot(ff(2,:),'r');plot(ff(3,:),'b')

legend('K1','O1','M2');

set(gca,'XTick',0:25:2020)

set(gca,'XTickLabel',{'1920','1945','1970','1995','2020'})

title('(a) nodal factor')

subplot(2,1,2);

plot(uu(1,:),'k');hold on;plot(uu(2,:),'r');plot(uu(3,:),'b')

set(gca,'XTick',0:25:2020)

set(gca,'XTickLabel',{'1920','1945','1970','1995','2020'})

title('(b) nodal angle')

結果如下：

需要注意的是輸入的分潮名必須用cell數組，365代表計算的時間間隔為1年（365天）。10代表緯度為北緯10度，這個參數不重要，只要不為0即可。