Skip to content

Instantly share code, notes, and snippets.

View Observatorio-de-Matematica's full-sized avatar

Observatorio-de-Matematica

View GitHub Profile
@Observatorio-de-Matematica
Observatorio-de-Matematica / Description.txt
Last active November 30, 2023 16:18 — forked from MSoegtropIMC/Description.txt
Example of main stream linear program solvers failing on small problems
The attached files contain a piecwise degree 3 polynomial L-inf norm approximation for the sine and cosine function
in the interval 0...1 with 4 and 64 pieces.
The functions are sampled at the extrema of a Chebyshev polynomial of degree 4 since the residuals of an L-inf polynomial
approximation are by construction typically close to Chebyshev polynomials. The small deviations in the extrema positions usually
don't lead to substantial deviations in the extrema values.
This leads to rather small problems with only 5 variables (4 coefficients and one maximum error variable) and 10 inequalities
(5 extrema points, each positive and negative bounded).
The solution is implemented in:
@Observatorio-de-Matematica
Observatorio-de-Matematica / gist:0cbb78dea7440b7a5d8fe16f9b422fda
Created November 30, 2023 16:19 — forked from vext01/gist:1437651
"image file is not of type 15"
/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/
/* [ Created with wxMaxima version 0.8.7 ] */
/* [wxMaxima: input start ] */
kill(all);
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
chull(z, a, b) := (a+b)*z - (a*b);
/* [wxMaxima: input end ] */

Kurvendiskussion für das Abi NRW mit Maxima

under development

Dies richtet sich an deutschsprachige Lehrerinnen und Lehrer und Schülerinnen und Schüler, die eine automatische Kurvendiskussion für Übungsaufgaben haben wollen. Dabei wird das freie Algebrasystem Maxima (wxMaxima) verwendet.

Funktionen/Routinen zur Kurvendiskussion

Eingabe einer Funktion

@Observatorio-de-Matematica
Observatorio-de-Matematica / Lagrange.wxm
Created November 30, 2023 16:40 — forked from FooBarrior/Lagrange.wxm
NM expirement. Author: Pousqie
/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/
/* [ Created with wxMaxima version 11.08.0 ] */
/* [wxMaxima: input start ] */
load(draw)$
Lagrange(x, fn, X, N) :=
sum(
fn(X[i]) * product(
if equal(i, j)
then 1
@Observatorio-de-Matematica
Observatorio-de-Matematica / q3.wxm
Created November 30, 2023 16:40 — forked from fernandomalmeida/q3.wxm
Resolução do Problema XOR utilizando função de base radial
/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/
/* [ Created with wxMaxima version 11.08.0 ] */
/* [wxMaxima: comment start ]
Essa primeira parte calcula os multiplicadores de lagrange
[wxMaxima: comment end ] */
/* [wxMaxima: input start ] */
rbf(x,i,y) := (%e^(-((x[i]-y).(x[i]-y))/2)) $
/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/
/* [ Created with wxMaxima version 12.01.0 ] */
/* [wxMaxima: input start ] */
load(diag);
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
a1:matrix([7,-1,6],[-10,4,-12],[-2,1,-1]);
/* [wxMaxima: input end ] */
/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/
/* [ Created with wxMaxima version 13.04.2 ] */
/* [wxMaxima: input start ] */
/*
get_M(_VARS,_F):=block([X,N,r,c,B,y,coeff,flag],
X:_VARS,
f:_F,
N:length(X),
M:[],
A.F.Beardon Complex Analysis page 13:
```maxima
(%i1) z: sqrt( 1/2*( sqrt(u^2+v^2) + u) ) + %i*(v/sqrt( 2*( sqrt(u^2+v^2) + u) ) );
2 2
sqrt(sqrt(v + u ) + u) %i v
(%o1) ----------------------- + -------------------------------
sqrt(2) 2 2
sqrt(2) sqrt(sqrt(v + u ) + u)
(%i2) z^2;
(%i73) assume(l>0);
(%o73) [l > 0]
(%i74) integrate( exp(-%i*l*t^2),t,-inf,inf);
1 %i
sqrt(%pi) (------- - -------)
sqrt(2) sqrt(2)
(%o74) -----------------------------
sqrt(l)
(%i75) polarform (%);
%i %pi
(%i13) limit(integrate( ev( (-%i*x)^(n-1)*exp(%i*x*(D+%i*e))/gamma(n), n:3), x,0,inf), e,0);
%i
(%o13) --
3
D
(%i14) limit(integrate( ev( (-%i*x)^(n-1)*exp(%i*x*(D+%i*e))/gamma(n), n:4), x,0,inf), e,0);
%i
(%o14) --
4
D