Ago77

A very basic version: wxdraw2d(yrange=[-5/2, 5/2], xaxis=true, yaxis=true, transparent=true, point_type=filled_circle, color=black, polygon([[0, 2], [2, 0], [0, -2], [-2, 0]]), color=blue, polygon([[0, 2], [4, 0], [0, -2], [-4, 0]]), color=red, points([[1, 1], [1, -1], [-1, -1], [-1, 1]]), color=green, points([[2, 1], [2, -1], [-2, -1], [-2, 1]]));

That indicates that Maxima cannot find a solution. Your main options are (1) try to see if you can simplify the equation, (2) give it a go with to_poly_solve, (3) use a numerical procedure, such as the one in the commands "newton" or "find_root". In any case, you would find local solutions, for a start, then for a global solution you would need to compare the local solution, or check if the function is concave, ...

coefmatrix works on a system of linear equations

solve wants lists as inputs, and matrices are not lists, so you have to build a new list: (A: matrix([1], [x]), B: matrix([1, 2],[3,4]), C: matrix([y], [2]), fx: A-B.C); solve([fx[1][1], fx[2][1]], [x,y]);

One solution is a parametric plot: wxdraw2d(xaxis=true, yaxis=true, explicit(sin(x), x, 0, 10), parametric(sin(x), x, x, 0, 10))

I am not sure about what problem you are tackling, but this would give you a solution, for example: solve(P(w,v)-[1,2,3], [a, v, w]); Basically, your definitions already give you a system, so if you put square brackets around them you obtain a strange object.

You should explicitly as the kernel to simplify the output trigsimp(integrate(exp(x)...

define(f(x), taylor(cosh(x), x, 1, 2)); wxdraw2d(key="approximated value", explicit(cosh(x),...