I wanted to rewrite my octave code to make it easier to read and make it a bit faster, just a simple vectorisation. I then run the two implementations and gets two very different results.
Did I fail my implementation?

I am trying to implement envelope tracking and pulse width modulation.
In this example, I have a 1 MHz sine wave. I sample this sine wave at 3MHz and each sample will determine a duty cycle for a pulse width signal with period of 0.33us.

I want to implement Romberg’s method in Octave.
Sadly calculating first line of matrix (T_0k method) is working extremely slow comparing to my own Python implementation.
I am using array istead of dictionary but in method mentioned before im not using any data structure.