0% encontró este documento útil (0 votos)
48 vistas15 páginas

Programas

El documento presenta comandos básicos de Mathematica como impresiones, operaciones aritméticas, asignación de variables, funciones, derivadas e integrales. También incluye estructuras de control como if, for y break para programación básica.
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
0% encontró este documento útil (0 votos)
48 vistas15 páginas

Programas

El documento presenta comandos básicos de Mathematica como impresiones, operaciones aritméticas, asignación de variables, funciones, derivadas e integrales. También incluye estructuras de control como if, for y break para programación básica.
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd

Hernández Calderón Paola Yvette

Comando Básicos de Mathematica


Cuestiones básicas y comentarios
(* Este es un comentario,
y como puedes ver se coloca entre paréntesis y asteriscos *)

Print["Así se puede imprimir en pantalla"] ;


escribe
Print[
escribe
"Si se finaliza una instrucción con punto y coma, no se despliega en pantalla,
a excepción de las impresiones en pantalla"] ;
Print["Para ejecutar un bloque de instrucciones es necesario
escribe
presionar SHIFT + ENTER"] ;

Así se puede imprimir en pantalla

Si se finaliza una instrucción con punto y coma, no


se despliega en pantalla, a excepción de las impresiones en pantalla
Para ejecutar un bloque de instrucciones es necesario presionar SHIFT + ENTER

Operadores aritméticos
Print["De esta forma se realizan las operaciones básicas"]
escribe
2+2
3-1
4  16
3*1
(*La siguiente operación no se muestra porque se terminó con punto y coma*)
2 + 2;

De esta forma se realizan las operaciones básicas


4

1
4
3

Printed by Wolfram Mathematica Student Edition


2 Métodos.nb

Asignación de variables
(* Para declarar una variable sólo es necesario asignarle un valor,
no es necesario definir su tipo, como en otros lenguajes de programación. *)
a = 1;
b = 2;
(* Sólo aparece el valor de c porque las
otras instrucciones se terminaron con punto y coma.*)
c = a+b

(* Para limpiar las variables se utiliza el siguiente comando*)


Clear[a, b]
borra

Funciones y sus gráficas


(* Definición de dos funciones que dependen de "x" *)
f[x_] := x ^ 2 - 2 x + 1;
g[x_] := Exp[x];
exponencial

(* Gráfica de f(x) *)
Plot[f[x], {x, - 2, 4}, PlotStyle → Magenta]
representación gráfica estilo de repr⋯ magenta
(* Gráfica de f(x) y g(x) *)
Plot[{f[x], g[x]}, {x, - 2, 4}, PlotStyle → {Blue, Magenta}]
representación gráfica estilo de repre⋯ azul magenta

-2 -1 1 2 3 4

Printed by Wolfram Mathematica Student Edition


Métodos.nb 3

20

15

10

-2 -1 1 2 3 4

(* Si colocas el cursor sobre un comando o palabra reservada,


aparece un pequeño cuadro emergente. *)

Printed by Wolfram Mathematica Student Edition


4 Métodos.nb

(* Y si das click en el signo que tiene una 'i',


se desplegará una ventana de ayuda, en la que se muestran diversos ejemplos. *)

Printed by Wolfram Mathematica Student Edition


Métodos.nb 5


Derivadas e integrales
Derivadas
(* Ya se definieron dos funciones, f(x) y g(x)*)
Print["Se calculan las primeras tres derivadas de f(x) con la notación f'(x)"]
escribe
f '[x]
f ''[x]
f '''[x]
Print["Se calculan las primeras tres derivadas de f(x) con el comando D[f(x)]"]
escribe deriva
D[f[x], x]
deriva
D[f[x], {x, 2}]
deriva
D[f[x], {x, 3}]
deriva

Se calculan las primeras tres derivadas de f(x) con la notación f'(x)


-2 + 2 x

Se calculan las primeras tres derivadas de f(x) con el comando D[f(x)]


-2 + 2 x

Printed by Wolfram Mathematica Student Edition


6 Métodos.nb

Integrales
Print"Se calcula la integral indefinida de ∫ f(x)ⅆx"
escribe
Integrate[f[x], x]
integra
2
Print"Se calcula la integral definida de  f(x)ⅆx"
escribe 0

Integrate[f[x], {x, 0, 2}]


integra

Se calcula la integral indefinida de ∫f(x)ⅆx


x3
x - x2 +
3
2
Se calcula la integral definida de  f(x)ⅆx
0

2
3

Comandos básicos de programación


if

Print["Sintaxis if → If[condición,verdadero,falso]"];
escribe si
valor = 10;
If[valor ≥ 10, Print["Es mayor o igual a 10"], Print["Es menor a 10"]];
si escribe escribe

Sintaxis if → If[condición,verdadero,falso]

Es mayor o igual a 10

if-else
Print["Sintaxis if-else → If[condición,verdadero,falso,ninguna de las anteriores]"];
escribe si
valor = "cosa";
If[valor ≥ 10, Print["Es mayor o igual a 10"],
si escribe
Print["Es menor a 10"], Print["No es un número"]];
escribe escribe

Printed by Wolfram Mathematica Student Edition


Métodos.nb 7

Sintaxis if-else → If[condición,verdadero,falso,ninguna de las anteriores]

No es un número

for y break
Print["Sintaxis for →
escribe
for[inicio,condición,incremento,instrucciones a ejecutar en cada iteración]"];

(* Imprime los números del 1 al 10 *)


For[i = 1, i ≤ 10, i ++, Print[i]];
para cada escribe

(* Obtiene las primeras 10 integrales indefinidas de f(x) *)


m[x] = f[x];
For[i = 1, i ≤ 10, i ++,
para cada
m[x] = Integrate[m[x], x];
integra
Print["Integral número " , i , " de f(x)", m[x]];
escribe integral
];

(* Imprime los números del 1 al 5, porque break rompe el ciclo cuando i⩵6*)
For[i = 1, i ≤ 10, i ++,
para cada
If[i > 5, Break[], Print[i]]
si finaliza it⋯ escribe
];

Printed by Wolfram Mathematica Student Edition


8 Métodos.nb

Sintaxis for → for[inicio,condición,incremento,instrucciones a ejecutar en cada iteración]

10
x3
Integral número 1 de f(x)x - x2 +
3
x2 x3 x4
Integral número 2 de f(x) - +
2 3 12
x3 x4 x5
Integral número 3 de f(x) - +
6 12 60
x4 x5 x6
Integral número 4 de f(x) - +
24 60 360
5 6
x x x7
Integral número 5 de f(x) - +
120 360 2520
6 7
x x x8
Integral número 6 de f(x) - +
720 2520 20 160
x7 x8 x9
Integral número 7 de f(x) - +
5040 20 160 181 440
x8 x9 x10
Integral número 8 de f(x) - +
40 320 181 440 1 814 400
x9 x10 x11
Integral número 9 de f(x) - +
362 880 1 814 400 19 958 400
10 11
x x x12
Integral número 10 de f(x) - +
3 628 800 19 958 400 239 500 800
1

Printed by Wolfram Mathematica Student Edition


Métodos.nb 9

Método de Bisección
(* Función y su gráfica para determinar el intervalo*)
f[x_] := x ^ 2 - 2 x - 5;
a = - 2;
b = - 1;
Plot[f[x], {x, a, b}, PlotStyle → Red]
representación gráfica estilo de repr⋯ rojo

(*Número máximo de iteraciones y tolerancia*)


maxiter = 100;
ϵ = 0.1;

(*Algoritmo*)

xma = a;(*Estimado anterior de la raíz para el cálculo del error *)


Fori = 1, i ≤ maxiter, i ++,
para cada
(*Punto medio*)
xm = Rounda + b  2, 0.001;
entero más próximo
(*Cálculo del error y condición de parada del método*)
ea = Abs[xm - xma];
valor absoluto
Print[i, " ", xm, " ", ea];
escribe
If[N[ea] ≤ ϵ, Break[],
si valor numér⋯ finaliza iteración
(*Determinación del nuevo intervalo*)
If[f[a] * f[xm] < 0, b = xm, a = xm, Break[]];
si finaliza iteración
];
(*Actualización de la aproximación anterior de la raíz*)
xma = xm;
;
Print["La raíz de f(x) es:", Round[xm, 0.001], " y se obtuvo con ", i, " iteraciones"]
escribe entero más próximo

-1.8 -1.6 -1.4 -1.2 -1.0

-1

-2

Printed by Wolfram Mathematica Student Edition


10 Métodos.nb

1 -1.5 0.5

2 -1.25 0.25

3 -1.375 0.125

4 -1.438 0.063

La raíz de f(x) es:-1.438 y se obtuvo con 4 iteraciones

Método de Regla Falsa


(* Función y su gráfica para determinar el intervalo*)
Clear[f]
borra
f[x_] := Exp[x] - 2;
exponencial
a = - 1;
b = 2;
Plot[f[x], {x, a, b}, PlotStyle → Magenta]
representación gráfica estilo de repr⋯ magenta

(*Número máximo de iteraciones y tolerancia*)


maxiter = 100;
ϵ = 0.001;

(*Algoritmo*)

xma = a;(*Estimado anterior de la raíz para el cálculo del error *)


Fori = 1, i ≤ maxiter, i ++,
para cada
(*Punto medio*)
f[b] * b - a
xm = Roundb - , 0.001;
f[b] - f[a]
entero más próximo

(*Cálculo del error y condición de parada del método*)


ea = Abs[xm - xma];
valor absoluto
Print[i, " ", xm, " ", ea];
escribe
If[N[ea] ≤ ϵ, Break[],
valor numér⋯ finaliza iteración
(*Determinación del nuevo intervalo*)
If[f[a] * f[xm] < 0, b = xm, a = xm, Break[]];
finaliza iteración
];
(*Actualización de la aproximación anterior de la raíz*)
xma = xm;
;
Print["La raíz de f(x) es:", Round[xm, 0.001], " y se obtuvo con ", i, " iteraciones"]
escribe entero más próximo

Printed by Wolfram Mathematica Student Edition


Métodos.nb 11

-1.0 -0.5 0.5 1.0 1.5 2.0

-1

-2

1 -0.303 0.697

2 0.134 0.437

3 0.39 0.256

4 0.532 0.142

5 0.609 0.077

6 0.649 0.04

7 0.67 0.021

8 0.681 0.011

9 0.687 0.006

10 0.69 0.003

11 0.692 0.002

12 0.693 0.001

La raíz de f(x) es:0.693 y se obtuvo con 12 iteraciones

Printed by Wolfram Mathematica Student Edition


12 Métodos.nb

Método de Punto Fijo


(*Función y su gráfica para determinar el punto inicial*)
f[x_] := Exp[- x] - x + 1;
exponencial
Plot[f[x], {x, - 5, 5}, PlotStyle → Green]
representación gráfica estilo de repr⋯ verde
g[x_] := f[x] + x;
(*Valor Inicial*)
x0 = 1.5;
(*Número máximo de iteraciones y tolerancia*)
maxiter = 100;
ϵ = 0.001;
(*Algoritmo*)
xa = x0(*Estimado anterior de la raíz*)
For[i = 1, i ≤ maxiter, i ++,
para cada
(*Aproximación de la raíz*)
xn = Round[g[xa], 0.0001];
entero más próximo
(*Cálculo del error y condición de parada*)
ea = N[Abs[xn - xa]];
⋯ valor absoluto
Print[i, " ", xn, " ", ea];
escribe
If[ea ≤ ϵ, Break[], xa = xn];
si finaliza iteración
];
Print["La raíz de f(x) es : " , N[xn], " y se obtuvo con " , i, " iteraciones"]
escribe valor numérico

80

60

40

20

-4 -2 2 4

1.5

Printed by Wolfram Mathematica Student Edition


Métodos.nb 13

1 1.2231 0.2769

2 1.2943 0.0712

3 1.2741 0.0202

4 1.2797 0.0056

5 1.2781 0.0016

6 1.2786 0.0005

La raíz de f(x) es : 1.2786 y se obtuvo con 6 iteraciones

Printed by Wolfram Mathematica Student Edition


14 Métodos.nb

Método de Newto Raphson en


Mathematica
(*Función y su gráfica para determinar el punto inicial*)
f[x_] := Exp[- x] - x + 1;
exponencial
Plot[f[x], {x, - 5, 5}, PlotStyle → Blue]
representación gráfica estilo de repr⋯ azul

(*Valor Inicial*)
x0 = 1.5;
(*Número máximo de iteraciones y tolerancia*)
maxiter = 100;
ϵ = 0.001;
(*Algoritmo*)
xa = x0(*Estimado anterior de la raíz*)
Fori = 1, i ≤ maxiter, i ++,
para cada
(*Aproximación de la raíz*)
f[xa]
xn = Roundxa - , 0.0001;
f '[xa]
entero más próximo

(*Cálculo del error y condición de parada*)


ea = N[Abs[xn - xa]];
⋯ valor absoluto
Print[i, " ", xn, " ", ea];
escribe
If[ea ≤ ϵ, Break[], xa = xn];
finaliza iteración
;
Print["La raíz de f(x) es : " , N[xn], " y se obtuvo con " , i, " iteraciones"]
escribe valor numérico
80

60

40

20

-4 -2 2 4

1.5

Printed by Wolfram Mathematica Student Edition


Métodos.nb 15

1 1.2736 0.2264

2 1.2785 0.0049

3 1.2785 0.

La raíz de f(x) es : 1.2785 y se obtuvo con 3 iteraciones

Printed by Wolfram Mathematica Student Edition

También podría gustarte