En general, la mayor parte de las ecuaciones se resuelven en Maxima con la función solve.
Empezamos por calcular una ecuación de primer grado, 4x−610+2x=21−3(x+1)12
solve( (4*x-6)/10 + 2*x = 21 - (3*(x+1))/12 );
[x=42753]
Una ecuación de segundo grado, x2−512x+6000=0
solve( x^2 - 512*x +6000=0);
[x=12,x=500]
Una ecuación de tercer grado, x3+3x2+x−6=0
s: solve( x^3 + 3*x^2 + x-6=0);
[x=2(√3i2−12)3(√6432332+52)13+(√6432332+52)13(−√3i2−12)−1,x=(√6432332+52)13(√3i2−12)+2(−√3i2−12)3(√6432332+52)13−1,x=(√6432332+52)13+23(√6432332+52)13−1]
Si la interpretación se hace dificultosa, siempre se puede pasar a formato decimal,
float(expand(s));
[x=−1.135939889088928i−2.047275740771164,x=1.135939889088928i−2.047275740771164,x=1.094551481542327]
Resolvemos un sistema lineal de dos incógnitas {u3+v−12=−33(u+v)+v=10−v
solve([u/3 + (v-1)/2 = -3, 3*(u+v) + v = 10 - v], [u, v] );
[[u=−105,v=65]]
Un sistema no lineal con parámetro, {xy=7x2+y2=a
solve([x*y = 7, x^2 +y^2 = a], [x,y]);
[[x=−√√a2−196+a√2,y=−7√2√√a−14√a+14+a],[x=√√a2−196+a√2,y=7√2√√a−14√a+14+a],[x=−√a−√a2−196√2,y=−7√2√a−√a−14√a+14],[x=√a−√a2−196√2,y=7√2√a−√a−14√a+14]]
Estudiamos algebraica y gráficamente los puntos de corte entre la parábola y=x2+3x+2 y la recta y=2x+3,
sol: solve([y=x^2 + 3*x + 2, y=2*x + 3], [x,y]);
[[x=−√5+12,y=2−√5],[x=√5−12,y=√5+2]]
Transformamos a decimales,
float(sol);
[[x=−1.618033988749895,y=−0.2360679774997898],[x=0.6180339887498949,y=4.23606797749979]]
Comprobamos gráficamente,
draw2d( grid = true, background_color = "#f0e68c", explicit(x^2 + 3*x + 2, x, -3, 2), color = red, explicit(2*x + 3, x, -3, 2) ) $
Para las inecuaciones, necesitamos cargar la librería fourier_elim. Empezamos con una inecuación de una sola incógnita, x4+5x3+5x2−5x−6≥0
load("fourier_elim") $ fourier_elim([x^4+5*x^3+5*x^2-5*x-6 > 0], [x]);
[1<x]∨[−2<x,x<−1]∨[x<−3]
Resolvemos un sistema de inecuaciones con dos incógnitas, {3x−5y<2x+y>1+x3y>2
des1 : 3*x - 5*y < 2 $ des2 : x+y> (1+x)/3 $ des3 : y < 2 $ fourier_elim([des1, des2, des3],[x,y]);
[12−3y2<x,x<5y3+23,−119<y,y<2]
Si el sistema es muy complicado siempre podemos echar mano de la resolución gráfica. Por ejemplo, {cos(2x)+cos(3y)<12x2+y2≤9
Cargamos el paquete draw y utilizamos el objeto gráfico region,
in1: cos(2*x)+cos(3*y)<1/2 $ in2: x^2+y^2 <=9 $ load("draw") $ draw2d( x_voxel = 30, y_voxel = 30, fill_color = navy, background_color = "#f0e68c", grid = true, region(in1 and in2 , x, -4, 4, y, -4, 4))$
© 2011-2016, TecnoStats.