Problem with fplot function - displays grey zone instead of functions

1 回表示 (過去 30 日間)
Marek Páleník
Marek Páleník 2021 年 10 月 21 日
コメント済み: Marek Páleník 2021 年 10 月 26 日
Hello All,
I would like to ask if You could help me with this problem described below?
I want to plot graphs defined by very long equations and I receive white and grey display field after script execution instead of 4 curves.
What am I doing wrong, may this be connected to gaps between the digits / letters in the equations?
Thank You for help in advance.
clear all
close all
clc
syms m k_11 k_12 k_22 Ia Ip W_v w
A=[k_11-m*W_v^2 -i*k_12;
i*k_12 -Ia*W_v^2+Ip*w*W_v+k_22];
Det=det(A)
f=solve(Det,W_v,'MaxDegree',4)
F=subs(f,[k_11, k_22, Ia, Ip, m, k_12], [100, 30, 5, 10, 3, 20])
fplot(@(w) w/2 - (6*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3) + 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3) + 282*w^2 - (9*w^4)/4 + ((3*w^2)/2 + 118/3)^2 + 2080)^(1/2)/(6*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/6)) - (3*6^(1/2)*(- w^3 + (82*w)/3)*(3*3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2) + 72*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) + 27*(- w^3 + (82*w)/3)^2 - 2*((3*w^2)/2 + 118/3)^3)^(1/2) - 12*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)*(6*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3) + 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3) + 282*w^2 - (9*w^4)/4 + ((3*w^2)/2 + 118/3)^2 + 2080)^(1/2) - ((3*w^2)/2 + 118/3)^2*(6*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3) + 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3) + 282*w^2 - (9*w^4)/4 + ((3*w^2)/2 + 118/3)^2 + 2080)^(1/2) - 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3)*(6*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3) + 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3) + 282*w^2 - (9*w^4)/4 + ((3*w^2)/2 + 118/3)^2 + 2080)^(1/2) + 12*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3)*(6*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3) + 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3) + 282*w^2 - (9*w^4)/4 + ((3*w^2)/2 + 118/3)^2 + 2080)^(1/2))^(1/2)/(6*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/6)*(6*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3) + 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3) + 282*w^2 - (9*w^4)/4 + ((3*w^2)/2 + 118/3)^2 + 2080)^(1/4)))
hold on
fplot(@(w) w/2 - (6*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3) + 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3) + 282*w^2 - (9*w^4)/4 + ((3*w^2)/2 + 118/3)^2 + 2080)^(1/2)/(6*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/6)) + (3*6^(1/2)*(- w^3 + (82*w)/3)*(3*3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2) + 72*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) + 27*(- w^3 + (82*w)/3)^2 - 2*((3*w^2)/2 + 118/3)^3)^(1/2) - 12*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)*(6*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3) + 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3) + 282*w^2 - (9*w^4)/4 + ((3*w^2)/2 + 118/3)^2 + 2080)^(1/2) - ((3*w^2)/2 + 118/3)^2*(6*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3) + 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3) + 282*w^2 - (9*w^4)/4 + ((3*w^2)/2 + 118/3)^2 + 2080)^(1/2) - 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3)*(6*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3) + 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3) + 282*w^2 - (9*w^4)/4 + ((3*w^2)/2 + 118/3)^2 + 2080)^(1/2) + 12*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3)*(6*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3) + 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3) + 282*w^2 - (9*w^4)/4 + ((3*w^2)/2 + 118/3)^2 + 2080)^(1/2))^(1/2)/(6*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/6)*(6*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3) + 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3) + 282*w^2 - (9*w^4)/4 + ((3*w^2)/2 + 118/3)^2 + 2080)^(1/4)))
fplot(@(w) w/2 + (6*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3) + 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3) + 282*w^2 - (9*w^4)/4 + ((3*w^2)/2 + 118/3)^2 + 2080)^(1/2)/(6*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/6)) - (12*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3)*(6*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3) + 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3) + 282*w^2 - (9*w^4)/4 + ((3*w^2)/2 + 118/3)^2 + 2080)^(1/2) - 12*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)*(6*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3) + 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3) + 282*w^2 - (9*w^4)/4 + ((3*w^2)/2 + 118/3)^2 + 2080)^(1/2) - ((3*w^2)/2 + 118/3)^2*(6*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3) + 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3) + 282*w^2 - (9*w^4)/4 + ((3*w^2)/2 + 118/3)^2 + 2080)^(1/2) - 3*6^(1/2)*(- w^3 + (82*w)/3)*(3*3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2) + 72*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) + 27*(- w^3 + (82*w)/3)^2 - 2*((3*w^2)/2 + 118/3)^3)^(1/2) - 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3)*(6*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3) + 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3) + 282*w^2 - (9*w^4)/4 + ((3*w^2)/2 + 118/3)^2 + 2080)^(1/2))^(1/2)/(6*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/6)*(6*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3) + 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3) + 282*w^2 - (9*w^4)/4 + ((3*w^2)/2 + 118/3)^2 + 2080)^(1/4)))
fplot(@(w) w/2 + (6*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3) + 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3) + 282*w^2 - (9*w^4)/4 + ((3*w^2)/2 + 118/3)^2 + 2080)^(1/2)/(6*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/6)) + (12*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3)*(6*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3) + 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3) + 282*w^2 - (9*w^4)/4 + ((3*w^2)/2 + 118/3)^2 + 2080)^(1/2) - 12*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)*(6*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3) + 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3) + 282*w^2 - (9*w^4)/4 + ((3*w^2)/2 + 118/3)^2 + 2080)^(1/2) - ((3*w^2)/2 + 118/3)^2*(6*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3) + 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3) + 282*w^2 - (9*w^4)/4 + ((3*w^2)/2 + 118/3)^2 + 2080)^(1/2) - 3*6^(1/2)*(- w^3 + (82*w)/3)*(3*3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2) + 72*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) + 27*(- w^3 + (82*w)/3)^2 - 2*((3*w^2)/2 + 118/3)^3)^(1/2) - 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3)*(6*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3) + 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3) + 282*w^2 - (9*w^4)/4 + ((3*w^2)/2 + 118/3)^2 + 2080)^(1/2))^(1/2)/(6*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/6)*(6*((3*w^2)/2 + 118/3)*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(1/3) + 9*((3^(1/2)*(128*((3*w^2)/2 + 118/3)^2*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^2 - 16*((3*w^2)/2 + 118/3)^4*(- (3*w^4)/16 + (47*w^2)/2 + 520/3) - 256*(- (3*w^4)/16 + (47*w^2)/2 + 520/3)^3 - 4*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)^3 + 27*(- w^3 + (82*w)/3)^4 + 144*(- w^3 + (82*w)/3)^2*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))^(1/2))/18 + (4*((3*w^2)/2 + 118/3)*(- (3*w^4)/16 + (47*w^2)/2 + 520/3))/3 + (- w^3 + (82*w)/3)^2/2 - ((3*w^2)/2 + 118/3)^3/27)^(2/3) + 282*w^2 - (9*w^4)/4 + ((3*w^2)/2 + 118/3)^2 + 2080)^(1/4)))

回答 (1 件)

Alan Stevens
Alan Stevens 2021 年 10 月 23 日
Change all the ^ to .^ (i.e. dot^)and the * to .*
  1 件のコメント
Marek Páleník
Marek Páleník 2021 年 10 月 26 日
Thanks. It would be too much of work to change this for each relevant mathematical operation, so I used a different way of computation and now it worked.
clear all
close all
clc
syms m k_11 k_12 k_22 Ia Ip W_v w
A=[k_11-m*W_v^2 -i*k_12;
i*k_12 -Ia*W_v^2+Ip*w*W_v+k_22];
Det=det(A)
f=solve(Det,W_v,'MaxDegree',4)
F=subs(f,[k_11, k_22, Ia, Ip, m, k_12], [1000, 30, 5, 10, 3, 20])
W_vector=-20:1:20;
for j=1:1:length(W_vector)
d=W_vector(j);
R(:,j)=subs(F,w,d);
end
R_abs=abs(R);
w_star=sqrt(1000/3)
asympt=W_vector*10/5;
figure(1)
grid on
plot(W_vector,R,'LineWidth',2)
hold on
plot(W_vector,W_vector,'.')
plot(W_vector,w_star*ones(1,length(W_vector)),'r*')
plot(W_vector,-w_star*ones(1,length(W_vector)),'r*')
plot(W_vector,asympt,'--')

サインインしてコメントする。

カテゴリ

Help Center および File ExchangeSymbolic Math Toolbox についてさらに検索

タグ

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by