Ayuda plotting!!
Publicado por Mary (27 intervenciones) el 17/07/2015 14:34:59
Hola!! Estoy resolviendo una EDO como veis ahí, que depende un parametro k, que va desde 0 hasta un valor que yo introduzco. QUiero dibujar la solucion en funcion de k, pero me da error por length!!!! Ayuda!!
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
%Me pide frecuencias w1 y w2
%Además va a pedirme hasta que k queremos evaluar (Klim)
w1=input('Introduce la frecuencia 1: ');
w2=input('Introduce la frecuencia 2: ');
klim=input('Introduce el valor de k que quieras limite: ');
a=(w2-w1);
x0(1)=0;
x0(2)=0;
tspan=[0 20];
for k =0:1:klim
x0=zeros(1,2);
%funcion que queremos integrar
f=@(t,x)[a-k*sin(x)];
%ode45 integra funcion
[t,x]=ode45(f,tspan,x0);
plot(k,x)
end
Valora esta pregunta
0
