Я работаю над уравнениями в частных производных для теплообмена. Есть 3 переменные, у, время и п. Я должен интегрировать функцию, как
f = sin(n*pi*y)dy on y between [0 L]
и вычислим сумму ряда
f_Integration*e^(n*pi*t)*sin(n*pi*y) with n [1..Inf].
Для каждого шага я пытался использовать функцию handle:
f = @(y,n)...
или symfun:
f(y,n) = ...
Я настроил:
syms n y t;
Я интегрирую с int:
int(ftoIntegrate(y,n), y , 0, L)
Я вычисляю сумму с помощью symsum:
symsum(ftoSum,n,1,Inf)
Если я использую plot или fplot, я ничего не наблюдаю. Должен ли я использовать векторы? Я не уверен, какой подход является правильным для того, что я хотел бы сделать.