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.
darova
darova 2020 年 3 月 27 日
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.

カテゴリ

ヘルプ センター および File Exchange で Symbolic Math Toolbox についてさらに検索

製品

リリース

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