Feeds
已回答
RK4 Function Code
clc a=0;b=2;N=10;alph=0.5; h=(b-a)/N;t(1)=a;w(1)=alph; f=@(t,y) y-t^2+1; for i=2:N+1 k1= h*f(t(i-1),w(i-1)); k2= h...
RK4 Function Code
clc a=0;b=2;N=10;alph=0.5; h=(b-a)/N;t(1)=a;w(1)=alph; f=@(t,y) y-t^2+1; for i=2:N+1 k1= h*f(t(i-1),w(i-1)); k2= h...
2 years 前 | 0
已回答
Euler's method
clc a=0;b=2;N=10; alph=0.5; h=(b-a)/N; t(1)=a; w(1)=alph; f = @(t,y) y-t^2+1; for i=2:N+1 t(i)=a+(i+1)*h; w(i)=w(i...
Euler's method
clc a=0;b=2;N=10; alph=0.5; h=(b-a)/N; t(1)=a; w(1)=alph; f = @(t,y) y-t^2+1; for i=2:N+1 t(i)=a+(i+1)*h; w(i)=w(i...
2 years 前 | 0
已回答
Composite Simpson's 1/3 rule code error
clc a=0;b=2;n=8; h=(b-a)/n; f=@(x) exp(2*x)*sin(3*x); XI0= f(a) +f(b); XI1 = 0; XI2 = 0; for i = 1:n-1 x = a+i*h; ...
Composite Simpson's 1/3 rule code error
clc a=0;b=2;n=8; h=(b-a)/n; f=@(x) exp(2*x)*sin(3*x); XI0= f(a) +f(b); XI1 = 0; XI2 = 0; for i = 1:n-1 x = a+i*h; ...
2 years 前 | 0
已回答
Secant
clc p0 = 0; p1 = 1; tol = 10^-6; N=12; i = 2; f = @(x) %write function; q0 = f(p0); q1 = f(p1); while i<=N p = p1- (q1...
Secant
clc p0 = 0; p1 = 1; tol = 10^-6; N=12; i = 2; f = @(x) %write function; q0 = f(p0); q1 = f(p1); while i<=N p = p1- (q1...
2 years 前 | 0
