[MatLab] "error: 'b' undefined..." ?

[kaspar]

Ciao

Ho un esercizio da svolgere con Matlab (io uso Octave che dovrebbe essere lo stesso per quello che mi serve.)

Ho definito una funzione (che non è importante esattamente sapere cosa fa...)

Codice:

function I = trapezi(f, a, b, n)
  H = (b-a)/n;
  S = 0; x = a;
  while x <= b
    S = S + (f(x) + f(x+H));
    x = x+H;
  end
  I = S*H/2;
end


Nella "Command Window" faccio questo:

Codice:

>> f = @(x) x*exp(-x)*cos(2*x);
>> met_quadr = @(met, n) met(f, 0, 2*pi, n);
>> met_quadr(trapezi, 10);
error: 'b' undefined near line 3 column 8
error: called from
    trapezi at line 3 column 5


Perché quell'errore? b è il nome di un parametro della funzione trapezi...

[Indrjo Dedej]

Ecco la traduzione in Python:

Codice:

# hello.py

from math import exp, pi, cos

def trapezi(f, a, b, n):
  H, S, x = (b-a)/n, 0, a
  while x <= b:
    S, x = S+f(x)+f(x+H), x+H
  return S*H/2

f = lambda x : x*exp(-x)*cos(2*x)

met_quadr = lambda met, n : met(f, 0, 2*pi, n)


Codice:

$ python3
>>> from hello import *
>>> met_quadr(trapezi, 10)
-0.15287622420461677


Quindi penso proprio che

Codice:

...
>> met_quadr = @(met, n) met(f, 0, 2*pi, n);
...


non fa quel che tu vuoi.

Ho un esercizio da svolgere con Matlab (io uso Octave che dovrebbe essere lo stesso per quello che mi serve.)

Anche a me viene lo stesso errore; può darsi che sia colpa di Octave, non so, mentre in Matlab il codice gira.

[feddy]

@OP l'esponenziale in MatLab, come in Python, si scrive

Codice:

exp(x)

, non

Codice:

e^

[kaspar]

@OP l'esponenziale in MatLab, come in Python, si scrive

Codice:

exp(x)

, non

Codice:

e^

Su GNU Octave quella notazione per l'esponenziale va bene. Evidentemente in MatLab no. Correggo.

[...] penso proprio che

Codice:

...
>> met_quadr = @(met, n) met(f, 0, 2*pi, n);
...


non fa quel che tu vuoi.

Capisco. Dovrò scrivere qualche riga in più, pazienza.


Grazie.

[feddy]

https://octave.org/doc/v4.2.1/Exponents ... ithms.html

Qui sembra che l'esponenziale sia uguale a matlab,python,ecc. Non ho Octave sul PC al momento, quindi non ho verificato.