两个M文件内容分别是:
1.
function [ output_args ] = Untitled2( input_args )
%UNTITLED2 Summary of this function goes here
% Detailed explanation goes here
clc;
clear all;
w=[];
z=[];
figure;
h=ezplot('(w^(4*z)*(1821744317201706*w^(2/5 - 2*z) + 1125899906842624*w^(4/5 - z)*(2*cos((pi*z)/2) + cos((pi*(z + 4/5))/2)) + 910872158600853*w^(6/5))^2*(2*log(w) + pi*w^(z + 2/5)*sin((pi*(z + 2/5))/2) + 2*w^(z + 2/5)*log(w)*cos((pi*(z + 2/5))/2)))/(63382530011411470074835160268800*(w^(2*z + 4/5) + 2*w^(z + 2/5)*cos((pi*(z + 2/5))/2) + 1)^3) - (20000*w^(4*z)*((910872158600853*w^(2/5 - 2*z))/562949953421312000 + (w^(4/5 -z)*(2*cos((pi*z)/2) + cos((pi*(z+ 4/5))/2)))/1000 + (910872158600853*w^(6/5))/1125899906842624000)*(w^(2/5 - 2*z)*log(w^(910872158600853/281474976710656000)) + (pi*w^(4/5 - z)*(2*sin((pi*z)/2) + sin((pi*(z + 4/5))/2)))/2000 + w^(4/5 - z)*log(w^(1/1000))*(2*cos((pi*z)/2) + cos((pi*(z + 4/5))/2))))/(w^(2*z+ 4/5) + 2*w^(z + 2/5)*cos((pi*(z + 2/5))/2) + 1)^2 - (20000*w^(4*z)*((529429888639651*w^(2/5 - 2*z))/450359962737049600 - (w^(4/5 -z)*(2*sin((pi*z)/2) - sin(1/2*pi*(z + 4/5))))/1000 + (5294298