Problemas para ejecutar un runge kutta
Publicado por Ray (1 intervención) el 05/06/2017 15:29:31
Buenos días, estoy haciendo un código para programar un método de Runge Kutta, pero me marca un error: 'In an assignment A(I) = B, the number of elements in B and I must be the same.'
mi código es el siguiente:
alguién puede encontrar el problema
mi código es el siguiente:
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
clc;
clear all;
N=10; a=0; b=2.5; Yini=3;
F='-1.2*Y + 7*exp(-0.3*X)';
Yp='70/9*exp(-0.3*Xpart)- 43/9*exp(-1.2*Xpart)';
FXY=inline(F);
Ypart=inline(Yp);
X(1)=a; Y(1)=Yini;
h=(b-a)/N;
for i=1:N
k1(i+1) = FXY(X(i),Y(i));
k2(i+1) = FXY(X(i)+0.5*h,Y(i)+ 0.5*h*k1);
k3(i+1) = FXY((X(i)+0.5*h),(Y(i)+ 0.5*h*k2));
k4(i+1) = FXY((X(i)+h),(Y(i)+ k3*h));
Y(i+1) = Y(i) + (1/6)*(k1(i)+2*k2(i)+2*k3(i)+k4(i))*h;
end
Xpart= a:h:b;
Yparti= Ypart(Xpart);
plot(X,Y,'--r', Xpart, Yparti)
alguién puede encontrar el problema
Valora esta pregunta
0
