can anyone help me run this, or get this running
1 view (last 30 days)
Show older comments
Relly Syam
on 2 Aug 2021
Answered: Walter Roberson
on 2 Aug 2021
clear; clc;
format long e
tic %bagian
n=10 ;
syms c0 c1 c2 c3 c4 c5 c6 c7 c8 c9 c10 c11 t r0 r1 r2 r3 r4 r5 r6 r7 r8 r9 r10 r11 c0=0;
c1=1/11; c2=2/11; c3=3/11; c4=4/11; c5=5/11; c6=6/11; c7=7/11; c8=8/11; c9=9/11; c10=10/11; c11=1;
r0=0; r1=1; r2=2; r3=3; r4=4; r5=5; r6=6; r7=7; r8=8; r9=9; r10=10; r11=11;
EvalAt = [c0, c1, c2, c3, c4, c5, c6, c7, c8, c9, c10, c11];
ktemp = arrayfun(@(EA) euler([r0, r1, r2, r3, r4, r5, r6, r7, r8, r9, r10, r11], EA).', EvalAt, 'uniform', 0);
ptemp=arrayfun(@(EA) int((3^(EA-t))*euler([r0, r1, r2, r3, r4, r5, r6, r7, r8, r9, r10, r11],t),t,0,EA).', EvalAt, 'uniform', 0);
E = horzcat(ktemp{:}).' ;
K = horzcat(ptemp{:}).' ;
Ek=E-K ;
Inv_Ek= inv(Ek) ;
F=[c0*3^c0;c1*3^c1;c2*3^c2;c3*3^c3;c4*3^c4;c5*3^c5;c6*3^c6;c7*3^c7;c8*3^c8;c9*3^c9;c10*3^c10;c11*3^c11];
C=Ek\F ;
Ua=@(x)(C(1)*euler(0,x)+C(2)*euler(1,x)+C(3)*euler(2,x)+C(4)*euler(3,x)+C(5)*euler(4,x)+C(6)*euler(5,x)+C(7)*euler(6,x)+C(8)*euler(7,x)+C(9)*euler(8,x)+C(10)*euler(9,x)+C(11)*euler(10,x)+C(12)*euler(11,x)) ;
Ue=@(x)((3^EA)*(1-e^(-EA)));
uaa=zeros(12,1) ;
uee=zeros(12,1) ;
xx=zeros(12,1) ;
k=0;
for i=1:12
uaa(i)=Ua(k);
uee(i)=Ue(k);
xx(i)=k;
k=k+.1;
end
Uap= uaa;
Uex= uee;
y=(abs(uaa-uee));
[xx uee uaa y];
uee uaa y
%gambar plot(xx,uee,'k',xx,uaa,':o') grid on legend({'eksak','aproksimasi N=10'},'Location','Northwest')
%plot(xx,y,'k') %legend({'error N=10'},'Location','Northwest') %grid on xlabel('x'),ylabel('Eksak, Aproksimasi'); toc
0 Comments
Accepted Answer
Walter Roberson
on 2 Aug 2021
format long e
n=10 ;
syms c0 c1 c2 c3 c4 c5 c6 c7 c8 c9 c10 c11 t r0 r1 r2 r3 r4 r5 r6 r7 r8 r9 r10 r11
c0=0;
c1=1/11; c2=2/11; c3=3/11; c4=4/11; c5=5/11; c6=6/11; c7=7/11; c8=8/11; c9=9/11; c10=10/11; c11=1;
r0=0; r1=1; r2=2; r3=3; r4=4; r5=5; r6=6; r7=7; r8=8; r9=9; r10=10; r11=11;
EvalAt = [c0, c1, c2, c3, c4, c5, c6, c7, c8, c9, c10, c11];
ktemp = arrayfun(@(EA) euler([r0, r1, r2, r3, r4, r5, r6, r7, r8, r9, r10, r11], EA).', EvalAt, 'uniform', 0);
ptemp=arrayfun(@(EA) int((3^(EA-t))*euler([r0, r1, r2, r3, r4, r5, r6, r7, r8, r9, r10, r11],t),t,0,EA).', EvalAt, 'uniform', 0);
E = horzcat(ktemp{:}).' ;
K = horzcat(ptemp{:}).' ;
Ek=E-K ;
%do not calculate inv(Ek) !!
%{
Inv_Ek= inv(Ek) ;
%}
F=[c0*3^c0;c1*3^c1;c2*3^c2;c3*3^c3;c4*3^c4;c5*3^c5;c6*3^c6;c7*3^c7;c8*3^c8;c9*3^c9;c10*3^c10;c11*3^c11];
C = double(Ek)\F ; %<=== important change!!
Ua=@(x)(C(1)*euler(0,x)+C(2)*euler(1,x)+C(3)*euler(2,x)+C(4)*euler(3,x)+C(5)*euler(4,x)+C(6)*euler(5,x)+C(7)*euler(6,x)+C(8)*euler(7,x)+C(9)*euler(8,x)+C(10)*euler(9,x)+C(11)*euler(10,x)+C(12)*euler(11,x)) ;
Ue = @(EA)((3^EA)*(1-exp(-EA))); %<=== recheck this !!!
uaa=zeros(12,1) ;
uee=zeros(12,1) ;
xx=zeros(12,1) ;
k=0;
for i=1:12
uaa(i)=Ua(k);
uee(i)=Ue(k);
xx(i)=k;
k=k+.1;
end
Uap= uaa;
Uex= uee;
y=(abs(uaa-uee));
[xx, uee, uaa, y];
uee, uaa, y
%gambar plot(xx,uee,'k',xx,uaa,':o') grid on legend({'eksak','aproksimasi N=10'},'Location','Northwest')
%plot(xx,y,'k') %legend({'error N=10'},'Location','Northwest') %grid on xlabel('x'),ylabel('Eksak, Aproksimasi'); toc
0 Comments
More Answers (0)
See Also
Categories
Find more on Time Series in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!