Conversion from MuPAD to MATLAB

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.

製品

リリース

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