WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Mein MATLAB Forum - goMatlab.de

Mein MATLAB Forum

  
Gast > Registrieren       Autologin?   

Partner:




Forum
      Option
[Erweitert]
  • Diese Seite per Mail weiterempfehlen
     


Gehe zu:  
Neues Thema eröffnen Neue Antwort erstellen

Problem beim Lösen eines Gleichungssystems mit sym und solv
 

Lihaku
Forum-Newbie

Forum-Newbie

Beiträge: 5
Anmeldedatum: 27.08.16
Wohnort: ---
Version: 2015a
     Beitrag Verfasst am: 27.08.2016, 21:58     Titel: Problem beim Lösen eines Gleichungssystems mit sym und solv
  Antworten mit Zitat      
Hi zusammen,

ich habe ein Problem mit dem Lösen eines Gleichungssystems. Ich habe 4 Gleichungen für 4 Unbekannte. Nun bekomme ich das Ganze leider nicht so aufgelöst, dass ich alle Unbekannten unabhängig von den anderen Unbekannten lösen kann. Meine Idee: Symbolisch das System über den solve Befehl lösen lassen (siehe code):

Code:

function [] = herleitung()

syms c1 c2 y1 y2 q1 q2 eps

y1 = sym('c1,c2');
y2 = sym('c1,c2');
c1 = sym('y1,q1,eps');
c2 = sym('y2,q2,eps');
   
sys = [q1 == c1*eps+(1-eps)*y1,...
       q2 == c2*eps+0.47*c2+y2*(1-eps),...
       y1 == 0.27*c1+0.1*c1*c1+0.12*c1*c2,...
       y2 == 0.47*c2+0.1*c2*c2+0.27*c1*c2];

[c_p1,c_p2,y_ad1,y_ad2] = solve(sys,c1,c2,y1,y2)
end
 


Als Antwort gibt mir Matlab folgende Fehlermeldung:

Code:

Warning: Support of strings that are not valid variable names or define a number will be removed in a
future release. To create symbolic expressions, first create symbolic variables and then use operations on
them.
> In sym>convertExpression (line 1536)
  In sym>convertChar (line 1441)
  In sym>tomupad (line 1198)
  In sym (line 177)
  In herleitung (line 5)
Warning: Support of strings that are not valid variable names or define a number will be removed in a
future release. To create symbolic expressions, first create symbolic variables and then use operations on
them.
> In sym>convertExpression (line 1536)
  In sym>convertChar (line 1441)
  In sym>tomupad (line 1198)
  In sym (line 177)
  In herleitung (line 6)
Warning: Support of strings that are not valid variable names or define a number will be removed in a
future release. To create symbolic expressions, first create symbolic variables and then use operations on
them.
> In sym>convertExpression (line 1536)
  In sym>convertChar (line 1441)
  In sym>tomupad (line 1198)
  In sym (line 177)
  In herleitung (line 7)
Warning: Support of strings that are not valid variable names or define a number will be removed in a
future release. To create symbolic expressions, first create symbolic variables and then use operations on
them.
> In sym>convertExpression (line 1536)
  In sym>convertChar (line 1441)
  In sym>tomupad (line 1198)
  In sym (line 177)
  In herleitung (line 8)
Error using sym.getEqnsVars>checkDuplicates (line 100)
One or more variables are used multiple times. Cannot distinguish between variables and equations. Use a
vector of equations and a vector of variables to avoid ambiguity.

Error in sym.getEqnsVars>checkVariables (line 95)
checkDuplicates(vars);

Error in sym.getEqnsVars (line 62)
    checkVariables(vars);

Error in solve>getEqns (line 450)
[eqns, vars] = sym.getEqnsVars(argv{:});

Error in solve (line 225)
[eqns,vars,options] = getEqns(varargin{:});

Error in herleitung (line 15)
[c_p1,c_p2,y_ad1,y_ad2] = solve(sys,c1,c2,y1,y2)
 


Anwendungsfehler? Komme auf jeden Fall hier nicht weiter :/ Vllt hat einer von euch eine Idee?
Private Nachricht senden Benutzer-Profile anzeigen


Harald
Forum-Meister

Forum-Meister

Beiträge: 24.501
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 27.08.2016, 23:16     Titel:
  Antworten mit Zitat      
Hallo,

was bezweckst du denn mit
Code:
y1 = sym('c1,c2');

?

Grüße,
Harald
Private Nachricht senden Benutzer-Profile anzeigen
 
Lihaku
Themenstarter

 Forum-Newbie

Forum-Newbie

Beiträge: 5
Anmeldedatum: 27.08.16
Wohnort: ---
Version: 2015a
     Beitrag Verfasst am: 28.08.2016, 08:31     Titel:
  Antworten mit Zitat      
Hatten den Part (siehe code) ursprünglich nicht drin. Hatte dann ne Idee vom Kumpel bekommen, dass diese Definition das Problem schon lösen könnte.

Sind die vier Zeilen unnötig, wenn ich weiter unten die Gleichungen definiere?

Code:

y1 = sym('c1,c2');
y2 = sym('c1,c2');
c1 = sym('y1,q1,eps');
c2 = sym('y2,q2,eps');
 
Private Nachricht senden Benutzer-Profile anzeigen
 
Harald
Forum-Meister

Forum-Meister

Beiträge: 24.501
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 28.08.2016, 09:42     Titel:
  Antworten mit Zitat      
Hallo,

Zitat:
dass diese Definition das Problem schon lösen könnte.

Welches Problem denn?
Die "Definition" ist eine ungültige Syntax, und ich habe keine Ahnung, was sie aussagen oder bezwecken soll.

Ich hatte gerade R2012b zur Hand, und da läuft der Code ohne diese "Definitionen" durch.
In aktuelleren Releases muss man bei solve die Variablen zu einem Vektor zusammenfassen.
Code:
[c_p1,c_p2,y_ad1,y_ad2] = solve(sys,[c1,c2,y1,y2])


Bei weiteren Problemen bitte auch das verwendete Release angeben.

Grüße,
Harald
Private Nachricht senden Benutzer-Profile anzeigen
 
Lihaku
Themenstarter

 Forum-Newbie

Forum-Newbie

Beiträge: 5
Anmeldedatum: 27.08.16
Wohnort: ---
Version: 2015a
     Beitrag Verfasst am: 28.08.2016, 17:50     Titel:
  Antworten mit Zitat      
Besten Dank schon mal soweit. Noch zwei kurze Fragen zu Matlabs Ergebnis:

Code:

c_p1 =
 
 (258775035589095428420141437620035475806759488469470188254908750233600*q1*q2 - 35202634128080312381624956704158978048840718093435710158642744393728*q2 - 258775035589095428420141437620035475806759488469470188254908750233600*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 1)*q1 - 206669587422584338751391431904392260854074955337207739598683032780800*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 1)*q2 + 107986236574204146516419852605326126889492041182203678532252467200000*q2*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 1)^2 - 107986236574204146516419852605326126889492041182203678532252467200000*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 1)*q2^2 + 63577860945889305201407868197128951702267844323500082303388457369600*q2^2 + 35995412191401382172139950868442042296497347060734559510750822400000*q2^3 + 143091726476695033549983563707263309151807111013707657295294575411200*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 1)^2 - 35995412191401382172139950868442042296497347060734559510750822400000*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 1)^3 + 70405268256160641394217971472183782694216335525498697542665921801231*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 1))/(45046024713657373635950747230750913349030661329598510016666491944960*q2)
 (258775035589095428420141437620035475806759488469470188254908750233600*q1*q2 - 35202634128080312381624956704158978048840718093435710158642744393728*q2 - 258775035589095428420141437620035475806759488469470188254908750233600*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 2)*q1 - 206669587422584338751391431904392260854074955337207739598683032780800*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 2)*q2 + 107986236574204146516419852605326126889492041182203678532252467200000*q2*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 2)^2 - 107986236574204146516419852605326126889492041182203678532252467200000*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 2)*q2^2 + 63577860945889305201407868197128951702267844323500082303388457369600*q2^2 + 35995412191401382172139950868442042296497347060734559510750822400000*q2^3 + 143091726476695033549983563707263309151807111013707657295294575411200*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 2)^2 - 35995412191401382172139950868442042296497347060734559510750822400000*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 2)^3 + 70405268256160641394217971472183782694216335525498697542665921801231*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 2))/(45046024713657373635950747230750913349030661329598510016666491944960*q2)
 (258775035589095428420141437620035475806759488469470188254908750233600*q1*q2 - 35202634128080312381624956704158978048840718093435710158642744393728*q2 - 258775035589095428420141437620035475806759488469470188254908750233600*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 3)*q1 - 206669587422584338751391431904392260854074955337207739598683032780800*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 3)*q2 + 107986236574204146516419852605326126889492041182203678532252467200000*q2*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 3)^2 - 107986236574204146516419852605326126889492041182203678532252467200000*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 3)*q2^2 + 63577860945889305201407868197128951702267844323500082303388457369600*q2^2 + 35995412191401382172139950868442042296497347060734559510750822400000*q2^3 + 143091726476695033549983563707263309151807111013707657295294575411200*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 3)^2 - 35995412191401382172139950868442042296497347060734559510750822400000*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 3)^3 + 70405268256160641394217971472183782694216335525498697542665921801231*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 3))/(45046024713657373635950747230750913349030661329598510016666491944960*q2)
 (258775035589095428420141437620035475806759488469470188254908750233600*q1*q2 - 35202634128080312381624956704158978048840718093435710158642744393728*q2 - 258775035589095428420141437620035475806759488469470188254908750233600*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 4)*q1 - 206669587422584338751391431904392260854074955337207739598683032780800*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 4)*q2 + 107986236574204146516419852605326126889492041182203678532252467200000*q2*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 4)^2 - 107986236574204146516419852605326126889492041182203678532252467200000*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 4)*q2^2 + 63577860945889305201407868197128951702267844323500082303388457369600*q2^2 + 35995412191401382172139950868442042296497347060734559510750822400000*q2^3 + 143091726476695033549983563707263309151807111013707657295294575411200*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 4)^2 - 35995412191401382172139950868442042296497347060734559510750822400000*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 4)^3 + 70405268256160641394217971472183782694216335525498697542665921801231*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 4))/(45046024713657373635950747230750913349030661329598510016666491944960*q2)
 
 
c_p2 =
 
 (112589990684262400*q2)/52917295621603353 - (112589990684262400*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 1))/52917295621603353
 (112589990684262400*q2)/52917295621603353 - (112589990684262400*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 2))/52917295621603353
 (112589990684262400*q2)/52917295621603353 - (112589990684262400*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 3))/52917295621603353
 (112589990684262400*q2)/52917295621603353 - (112589990684262400*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 4))/52917295621603353
 
 
y_ad1 =
 
 (35278197081068857*q1)/35278197081068902 - (23694563539503032705*q2)/75606484749412357631816958442750443 + (709884336127808464838153795010560000*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 1)^3)/(4000890704394899234226445659490326653887018911035379*q2) - (709884336127808464838153795010560000*q2^2)/4000890704394899234226445659490326653887018911035379 + (root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 1)*(35995412191401382172139950868442042296497347060734559510750822400000*q2^2 + 68889862474194779583797143968130753618024985112402579866227677593600*q2 + 86258345196365142806713812540011825268919829489823396084969583411200*q1 - 23468422752053547131405990490727927564738778508499565847555307267077))/(67623086704983167228326207107135434778202699189235212556928635491285076376519966720*q2) - (5*(141976867225561692967630759002112000*q2 + 188132448130071777305568346411234727)*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 1)^2)/(1333630234798299744742148553163442217962339637011793*q2) + 11165549376158282233/64345944467584954819270139455733760
 (35278197081068857*q1)/35278197081068902 - (23694563539503032705*q2)/75606484749412357631816958442750443 + (709884336127808464838153795010560000*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 2)^3)/(4000890704394899234226445659490326653887018911035379*q2) - (709884336127808464838153795010560000*q2^2)/4000890704394899234226445659490326653887018911035379 + (root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 2)*(35995412191401382172139950868442042296497347060734559510750822400000*q2^2 + 68889862474194779583797143968130753618024985112402579866227677593600*q2 + 86258345196365142806713812540011825268919829489823396084969583411200*q1 - 23468422752053547131405990490727927564738778508499565847555307267077))/(67623086704983167228326207107135434778202699189235212556928635491285076376519966720*q2) - (5*(141976867225561692967630759002112000*q2 + 188132448130071777305568346411234727)*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 2)^2)/(1333630234798299744742148553163442217962339637011793*q2) + 11165549376158282233/64345944467584954819270139455733760
 (35278197081068857*q1)/35278197081068902 - (23694563539503032705*q2)/75606484749412357631816958442750443 + (709884336127808464838153795010560000*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 3)^3)/(4000890704394899234226445659490326653887018911035379*q2) - (709884336127808464838153795010560000*q2^2)/4000890704394899234226445659490326653887018911035379 + (root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 3)*(35995412191401382172139950868442042296497347060734559510750822400000*q2^2 + 68889862474194779583797143968130753618024985112402579866227677593600*q2 + 86258345196365142806713812540011825268919829489823396084969583411200*q1 - 23468422752053547131405990490727927564738778508499565847555307267077))/(67623086704983167228326207107135434778202699189235212556928635491285076376519966720*q2) - (5*(141976867225561692967630759002112000*q2 + 188132448130071777305568346411234727)*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 3)^2)/(1333630234798299744742148553163442217962339637011793*q2) + 11165549376158282233/64345944467584954819270139455733760
 (35278197081068857*q1)/35278197081068902 - (23694563539503032705*q2)/75606484749412357631816958442750443 + (709884336127808464838153795010560000*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 4)^3)/(4000890704394899234226445659490326653887018911035379*q2) - (709884336127808464838153795010560000*q2^2)/4000890704394899234226445659490326653887018911035379 + (root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 4)*(35995412191401382172139950868442042296497347060734559510750822400000*q2^2 + 68889862474194779583797143968130753618024985112402579866227677593600*q2 + 86258345196365142806713812540011825268919829489823396084969583411200*q1 - 23468422752053547131405990490727927564738778508499565847555307267077))/(67623086704983167228326207107135434778202699189235212556928635491285076376519966720*q2) - (5*(141976867225561692967630759002112000*q2 + 188132448130071777305568346411234727)*root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 4)^2)/(1333630234798299744742148553163442217962339637011793*q2) + 11165549376158282233/64345944467584954819270139455733760
 
 
y_ad2 =
 
 root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 1)
 root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 2)
 root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 3)
 root(z^4 - z^3*(4*q2 + 564397344390215331916705039233704181/141976867225561692967630759002112000) + z^2*((2041375088234133656059057877954261961*q1)/283953734451123385935261518004224000 + (303915428452917597867262642625774157*q2)/28395373445112338593526151800422400 + 6*q2^2 - 70405268256160641394217971472183782694216335525498697542665921801231/35995412191401382172139950868442042296497347060734559510750822400000) - z*((51218998470714384076084713987190365605842836661210903903186815348411*q2)/35995412191401382172139950868442042296497347060734559510750822400000 + (2041375088234133656059057877954261961*q1*q2)/141976867225561692967630759002112000 + (192280321622647140417501615508089861*q2^2)/20282409603651670423947251286016000 + 4*q2^3) + (2041375088234133656059057877954261961*q1*q2^2)/283953734451123385935261518004224000 + q2^4 + (38378914534751141138088494387778696545738465829883293*q2^2)/31970348316613501335705605715668691054973524377600000 + (390782453484157325502903134661462423*q2^3)/141976867225561692967630759002112000, z, 4)
 


Erste Verständnisfrage: Wieso führt Matlab das Ergebnis für eine der Variablen erneut als Vektor aus - also insgesamt als Matrix (4x4 für 4 Variablen und 4 mal demselben Ergebnis für eine Variable)?

Zweite Verständnisfrage: Matlab führt die Variable "z" ein, da es unendlich viele Lösungen gibt, richtig?
Private Nachricht senden Benutzer-Profile anzeigen
 
Harald
Forum-Meister

Forum-Meister

Beiträge: 24.501
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 28.08.2016, 19:46     Titel:
  Antworten mit Zitat      
Hallo,

1. Für mich sieht das nach jeweils einem länglichen Ausdruck aus. Falls dem tatsächlich nicht so ist, dann gibt es vielleicht mehrere Lösungszweige.

2. Genau.

Grüße,
Harald
Private Nachricht senden Benutzer-Profile anzeigen
 
Lihaku
Themenstarter

 Forum-Newbie

Forum-Newbie

Beiträge: 5
Anmeldedatum: 27.08.16
Wohnort: ---
Version: 2015a
     Beitrag Verfasst am: 28.08.2016, 21:31     Titel:
  Antworten mit Zitat      
Alles klar. Besten Dank für die Hilfe! Smile
Private Nachricht senden Benutzer-Profile anzeigen
 
Neues Thema eröffnen Neue Antwort erstellen



Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:

 Du kannst Beiträge in dieses Forum schreiben.
Du kannst auf Beiträge in diesem Forum antworten.
Du kannst deine Beiträge in diesem Forum nicht bearbeiten.
Du kannst deine Beiträge in diesem Forum nicht löschen.
Du kannst an Umfragen in diesem Forum nicht mitmachen.
Du kannst Dateien in diesem Forum posten
Du kannst Dateien in diesem Forum herunterladen
.





 Impressum  | Nutzungsbedingungen  | Datenschutz | FAQ | goMatlab RSS Button RSS

Hosted by:


Copyright © 2007 - 2025 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.