Can anybody help me to convert the following MuPAD functions to MATLAB Symbolic Math?
//Peano Th and kernel
//pp - truncated power
pp:=(z,n)->piecewise([z>=0,z^n],[Otherwise,0]):
//PeanoKernel - compute Peano kernel given L, dex and var
PeanoKernel:=proc(L,d,t)
local K;
begin
K:=1/d!*eval(L(x->pp(x-t,d)));
return(K);
end_proc:
//PeanoCorr - Corollary to Peano Theorem
PeanoCorr:=proc(L,d)
local K;
begin
K:=L(x->x^(d+1));
K:=1/(d+1)!*factor(simplify(K))*(D@@(d+1))(f)(`ξ`);
return(K);
end_proc:
//PeanoKernel:=proc(L,d,t) //versiune veche
// local K,ff;
//begin
// ff:=x-->piecewise([x-t>=0,(x-t)^d],[Otherwise,0]);
// K:=1/d!*eval(L(x->ff(x)));
// return(K);
//end_proc:
//PeanoEstimation - Estimation based on Peano's theorem
PeanoEstimation:=proc(L,d,a,b)
local K;
begin
K:=PeanoKernel(L,d,t);
M[d+1]*int(abs(K),t=a..b);
end_proc:
The automatic conversion yields a lot of errors.

댓글 수: 4

darova
darova 2020년 3월 26일
I think faster be to write new code. DO you how original formulas?
Radu Trimbitas
Radu Trimbitas 2020년 3월 26일
I tried but I had a lot of errors.
If L is a linear and continuous functional, and L(f)=0 for each polynomial f of degree at most d
, where K is the Peano Kernel
, and is z for and 0 otherwise.
I intend to use these formulas to compute rests for various approximation formulas.
The code works fine in MuPAD, but Mathworks removed the MuPAD Notebook Interface in 2020 a prerelease.
What about this
syms t
K = 1/factorial(d)*L( (x-t)^d );
L = integral(K*f^(d+1),a,b)
Radu Trimbitas
Radu Trimbitas 2020년 3월 28일
f^(d+1) means the derivative of order (d+1) of f

댓글을 달려면 로그인하십시오.

답변 (1개)

Cris LaPierre
Cris LaPierre 2020년 3월 28일

0 개 추천

Have you seen this page? It recommends trying the convertMuPADNotebook function. It also points to the resource page Convert MuPAD Notebooks to MATLAB Live Scripts.

제품

릴리스

R2019b

질문:

2020년 3월 25일

답변:

2020년 3월 28일

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by