can anyone help me graph numerical solutions and exact solutions

2 vues (au cours des 30 derniers jours)
Relly Syam
Relly Syam le 19 Mai 2021
Commenté : Relly Syam le 19 Mai 2021
%
clear;
clc;
format long g
syms c0 c1 c2 c3 c4 c5 c6 c7 c8 c9 c10 t r0 r1 r2 r3 r4 r5 r6 r7 r8 r9 r10
c0=0; c1=1/10; c2=2/10; c3=3/10; c4=4/10; c5=5/10; c6=6/10; c7=7/10; c8=8/10; c9=9/10; c10=10/10; r0=0; r1=1; r2=2; r3=3; r4=4; r5=5; r6=6; r7=7; r8=8; r9=9; r10=10;
EvalAt = [c0, c1, c2, c3, c4, c5, c6, c7, c8, c9, c10];
ktemp = arrayfun(@(EA) euler([r0, r1, r2, r3, r4, r5, r6, r7, r8, r9, r10], EA).', EvalAt, 'uniform', 0);
ptemp=arrayfun(@(EA) int((t-EA)*euler([r0, r1, r2, r3, r4, r5, r6, r7, r8, r9, r10],t),t,0,EA).', EvalAt, 'uniform', 0);
E = horzcat(ktemp{:}).'
K = horzcat(ptemp{:}).'
Ek=E-K
Inv_Ek= inv(Ek)
F=[1; 1; 1; 1; 1; 1; 1; 1; 1; 1; 1];
C=inv(Ek)*F
%solusi aproximasinya
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))
Ue=@(x)(cos(x))
uaa=[];
xx=[];
k=0;
for i=1:11
uaa(i)=Ua(k);
uee(i)=Ue(k);
xx(i)=k;
k=k+.1;
end
Uap= uaa'
Uex= uee'
Err=abs(uaa'-uee')
y=(abs(uaa-uee));
[xx uee uaa y];
uee;
uaa;
y;
plot(xx,uaa,'o',xx,uee,'r')
grid on
plot(xx,y)
grid on
toc
  3 commentaires
Relly Syam
Relly Syam le 19 Mai 2021
i just found a graph like this, and i think it's just a graph error solution. Could you please help me graph the exact and numerical solutions
Relly Syam
Relly Syam le 19 Mai 2021
I got it now, thank you very much

Connectez-vous pour commenter.

Réponses (0)

Catégories

En savoir plus sur Networks dans Help Center et File Exchange

Community Treasure Hunt

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

Start Hunting!

Translated by