В системе компьютерной алгебры MAPLE надо реализовать процедуру вычисления коэффициентов и носителей следующего степенного ряда:
Вот, что получилось:
Код:
restart;
with(LinearAlgebra);
with(VectorCalculus);
with(plots);
alpha := Vector([m[1], m[2]]);
beta:= Vector([d-m[1], d-m[2]]);
z[1] := -nu_1;
z[2] := -nu_2;
z := Vector([z[1], z[2]]);
t1 := DotProduct(alpha, z);
t2 := DotProduct(beta, z);
d := 3; m[1] := 1; m[2] := 2; j := 1; n := 0; t := 5; for k from 0 to t do n := n+k end do;
p1 := 0; for k from 0 to t do for nu_1 from 0 to k do nu_2 := k-\nu_1;
p1 := p1+{(-1)^(nu_1 +nu_2)}/{(nu_1! * nu_2! * d^(nu_1+nu_2]-1))}* mul(-d* s+m_1*nu_1+m_2 * nu_2+1, s = 1 .. nu_1+nu_2-1)\cdot x_1^(nu_1) *x_2^(nu_2) end do end do;
print(p1/d);
Отсюда получился ряд:
Теперь мне надо как-то изобразить на графике носители ряда. То есть, на графике должны быть точки (0,0), (0,1), (1,0) и т.д. Я думала реализовать это через запись в матрицу и потом уже построить график этой матрицы. Но никак не получается сделать запись в матрицу в цикле. MAPLE либо не хочет создавать матрицу, либо для каждой пары
создает новую матрицу.
Не подскажите, как можно это решить?