Integrales dobles para octave, usando comandos de matlab
Publicado por Juan123 (3 intervenciones) el 26/11/2015 03:23:45
He estado tratando de hallar la manera de resolver una integral doble en octave, ya que según noté el comando "int" de matlab no funciona en octave (pero si la función "quad").
El algoritmo que tengo para matlab es:
...La cual no tiene ningún problema en matlab, pero en octave, buscando en internet "int" equivale a "quad", sin embargo quad no resuelve cuando los límites de integración son variables.
Por ejemplo si declaro en octave:
>> fun = @(x,y) 1./(1+x.^2+y^2);
>> Integral = dblquad(fun,−1,1,0,1)
Halla el valor de la integral de la función fun, sin problemas.
Pero si en vez de poner 0 y 1 en los límites de la 1era integral pusiera x y x+2, (osea para a=x, b=x+2, en el código de matlab) por ejemplo. ¿Cómo declaro esos límites de integración variables en octave??..Lo he intentado declarando global x y, etc. pero no llego a nada, también modificando el código, pero siempre me sale: error: quadcc: lower limit of integration (A) must be a single real scalar.
He intentado y buscado varias veces resolver este problema, estaré muy agradecido por su ayuda. Saludos cordiales :(
El algoritmo que tengo para matlab es:
1
2
3
4
5
6
7
8
9
10
11
12
syms x y
i=input('Ingrese el subíndice i de la integral I a calcular: i=');
f=input('Ingrese la función f a integrar: f(x,y)=','s');
F=inline(f,'x','y');
disp('Ingrese los extremos del intervalo' );
a=input('Desde (y) :');
b=input('Hasta (y) :');
a1=input('Desde (x) :');
b1=input('Hasta (x) :');
fprintf('La integral pedida es I_%d es',i);
I=int(int(f,y,a,b),x,a1,b1)
...La cual no tiene ningún problema en matlab, pero en octave, buscando en internet "int" equivale a "quad", sin embargo quad no resuelve cuando los límites de integración son variables.
Por ejemplo si declaro en octave:
>> fun = @(x,y) 1./(1+x.^2+y^2);
>> Integral = dblquad(fun,−1,1,0,1)
Halla el valor de la integral de la función fun, sin problemas.
Pero si en vez de poner 0 y 1 en los límites de la 1era integral pusiera x y x+2, (osea para a=x, b=x+2, en el código de matlab) por ejemplo. ¿Cómo declaro esos límites de integración variables en octave??..Lo he intentado declarando global x y, etc. pero no llego a nada, también modificando el código, pero siempre me sale: error: quadcc: lower limit of integration (A) must be a single real scalar.
He intentado y buscado varias veces resolver este problema, estaré muy agradecido por su ayuda. Saludos cordiales :(
Valora esta pregunta
0