Aquí viene un dibujo de un creíble (Delphi Pascal) fuente (-:
Viewport:
xmin := -4; xmax := +4;
ymin := -4; ymax := +4;
La asignación de Color:
$$
\begin{cases}
\color{red}{3x^3+4y^3=7}\\
\color{green}{4x^4+3y^4=16}
\end{casos}
$$
Espero que esto ayude. (Ahora pensando en el resto ..)
EDIT. Un parámetro de la representación de la curva de $\;3x^3+4y^3=7\;$ es:
$$
x(t) = -\sqrt[3]{\frac{t}{3}} \qquad ; \qquad y(t) = \sqrt[3]{\frac{t+7}{4}}
$$
Esto convierte el problema en la búsqueda de ceros de una función unidimensional:
$$
f(t) = 4 x^4(t) + 3 y^4(t) - 16
$$
Estándar de métodos numéricos puede ser usada para este propósito.
Newton-Raphsony
Regula Falsi
han tenido éxito aquí. Nuestro (doble precisión) de los resultados son:
$$
t = 3.56874527617414 \quad \Longrightarrow \quad
(x,y) = ( -1.05957433921527 , +1.38246606572819 ) \\
t = -8.23503156460535 \quad \Longrightarrow \quad
(x,y) = ( +1.40017183631184 , -0.675884813969605 )
$$
Dar, respectivamente:
$$
x+y = 0.322891726512912 \\
x+y = 0.724287022342236
$$
ALGORITMO (Delphi Pascal) :
programa de RF;
tipo de
funktie = function(t : doble) : doble;
función regula_falsi(d1,d0,eps : doble; F : funktie) : doble;
{
Regula Falsi
}
var
OK : boolean;
d2 : doble;
comenzar
d2 := d0;
mientras abs(F(d2)) > eps hacer
comenzar
d2 := d1 F(d1)*(d0 - d1)/(F(d0) - F(d1));
si F(d2) = 0 then Break;
OK := (F(d2)*F(d1) < 0);
si no, a continuación, OK
comenzar
d1 := d2;
end else begin
d0 := d1; d1 := d2;
end;
Escribir(d2,' : ',F(d2),' ; Pulse Enter'); Readln;
end;
regula_falsi := d2;
end;
función de la potencia(x,r : doble) : doble;
var
M : doble;
comenzar
M := 0;
si x > 0 entonces
M := exp(r*ln(abs(x)));
si x < 0 entonces
M := -exp(r*ln(abs(x)));
potencia := M;
end;
la función original(t : doble) : doble;
{
La función en sí
}
var
x,y,f : doble;
comenzar
x : = potencia(t/3,1/3);
y := potencia((t+7)/4,1/3);
f := 4*sqr(sqr(x))+3*sqr(sqr(y))-16;
original := f
end;
procedimiento de Calcular(x1,x2 : double);
const
eps : double = 1.E-14;
var
t,x,y,t1,t2 : doble;
comenzar
t1 := 3*sqr(x1)*x1;
t2 := 3*sqr(x2)*x2;
t := Regula_Falsi(t1,t2,eps,original);
Writeln(t);
x : = potencia(t/3,1/3);
y := potencia(t/4+7/4,1/3);
Writeln(x,' +',y,' =',x+y);
end;
comenzar
Calcular( 1.0, 1.5);
Writeln;
Calcular(-1.5,-0.5);
final.
Tenga en cuenta que el único lugar que es realmente dependiente de la máquina es
el criterio de parada con 'eps'. Como para el resto, el algoritmo es
bastante general. También tome nota de la más cruda a partir de los valores de $(1.0,1.5)$ $(-1.5,-0.5)$ $x$ al final del código, haciendo hincapié en la solidez de este algoritmo.