Problema Matlab con fattoriale

[ireon]

Se io scrivo in Matlab le seguenti istruzioni:

syms x;
f = (x+1)^x-factorial(x);

Mi appare il messaggio d'errore "N must be a matrix of non negative integers"

Quindi come faccio a scrivere correttamente la funzione?

[g.schgor]

Prova a definire x come intero e limitarlo ai primi numeri (4 o 5)

[DirtyDeeds]

Prova con queste dichiarazioni:

Codice: Seleziona tutto

syms x positive;
assumeAlso(x,'integer');


[ireon]

Mi appare scritto "Undefined function 'assumeAlso' for input arguments of type 'sym'."

[dimaios]

Sei sicuro di cosa hai scritto ?
Prova a copiare sul workspace queste righe di codice ( cut & paste e premi enter ).

Codice: Seleziona tutto

clear all;
syms x positive;
assumeAlso(x,'integer');
f = (x+1)^x-factorial(x);


A me funziona correttamente.

[ireon]

Mi appare scritto il messaggio d'errore, posto l'immagine:

[dimaios]

Verifica se i comandi esistono nella tua versione di Matlab.
Secondo le release notes di Matlab risulterebbe che ....
R2012a
New Features, Bug Fixes, Compatibility Considerations
expand all
New Special Functions
New Vector Analysis Functions
Computations with Symbolic Functions
Assumptions on Variables

You can set assumptions on symbolic variables by using these functions:

assume sets assumptions on symbolic variables.

assumeAlso adds assumptions on symbolic variables without erasing the previous assumptions.

assumptions shows assumptions set on symbolic variables.

New Relational Operators Create Equations, Inequalities, and Relations
New Logical Operators Create Logical Expressions

Prova a digitare :

Codice: Seleziona tutto

help assumeAlso

e vedi se la funzione esiste nella tua versione di Matlab che dall'immagine risulta essere la 2011b.

[ireon]

Non esiste nella mia versione, mi appare scritto "assumeAlso not found." Non c'è un altro modo per sottrarre una quantità fattoriale incognita?

[dimaios]

Fai un tentativo con questo :

Codice: Seleziona tutto

syms x positive;
f = (x+1)^x-factorial(round(x));

[ireon]

Ho risolto con il seguente codice:

syms x;
f = (x+1)^x-sym('n!');