Скажем 0.1 сек на генерацию.

Снова меня терзают смутные сомнения, что посчитать "честно" 10к точек за такое время на дохлом компе нереально. Тогда останется только интерполяция степеней косинусов. Впрочем, дерзайте, мне лень оценивать выч.затраты.

Меня интересует как мне менять форму этой кривой использую приведенную формулу (или что получше).

Я привел ссылку на статью. Там есть формулы с понятными параметрами. Например, (30). Ваша же модель слишком жесткая. Надо отдельно вычислять X и Y для заданного фазового угла - это и есть параметрический способ задания функции. Можно взять X=H_x * [cos(alpha)]^3+H_c*sin(alpha), Y=M_s*cos(alpha), M_s - намагниченность насыщения, H_c - коэрцитивная сила, H_x - напряженность поля в точке насыщения (максимальное значение намагничивающего поля), 0<alpha<2*Pi - параметр. Пробуйте.

Ппоробую завтра видео сделать из того что есть.

Ненене, не надо видео. "Ша, уже никто никуда не идет!" (с)