matalb-ODE45验证_龙格库塔_periodafw_ode45振动_matlab_非线性振动

ezplot('(-x^2+1+0.8668339830e-5*y^2)^2+3.031487314*10^(-15)*x^2 = .7958731346/(y^2*(1+5.34806421*10^(-8)*y^2/x^4+3.354847104*10^(-12)*y^4/x^8+1.462004114*10^(-23)*y^8/(x^4*(x^2+1+0.8668339838e-5*y^2)^4)+5.503426430*10^(-14)*y^4/x^4-1.069612838*10^(-7)*y^2/x^6-3.536880025*10^(-12)*y^4/x^6-2.120247581*10^(-18)*y^6/x^6+1.596070755*10^(-22)*y^8/(x^4*(x^2+1+0.8668339838e-5*y^2)^5)+9.689143957*10^(-21)*y^8/(x^8*(x^2+1+0.8668339838e-5*y^2))-5.528011082*10^(-12)*y^4/(x^2*(x^2+1+0.8668339838e-5*y^2)^3)-9.396672956*10^(-29)*y^10/(x^8*(x^2+1+0.8668339838e-5*y^2))-1.527536481*10^(-21)*y^8/(x^6*(x^2+1+0.8668339838e-5*y^2)^3)-1.039841538*10^(-17)*y^6/(x^8*(x^2+1+0.8668339838e-5*y^2)^2)+6.845352599*10^(-23)*y^8/(x^2*(x^2+1+0.8668339838e-5*y^2)^5)+7.514677917*10^(-24)*y^8/(x^6*(x^2+1+0.8668339838e-5*y^2)^4)+1.918766385*10^(-15)*y^6/((x^2+1+0.8668339838e-5*y^2)*x^8)+1.105602216*10^(-11)*y^4/(x^4*(x^2+1+0.8668339838e-5*y^2)^3)-6.203967184*10^(-23)*y^8/(x^6*(x^2+1+0.8668339838e-5*y^2)^2)+4.891934604*10^(-14)*y^4/(x^2*(x^2+1+0.8668339838e-5*y^2))+1.091008166*10^(-28)*y^10/(x^2*(x^2+1+0.8668339838e-5*y^2)^7)+7.385688311*10^(-40)*y^14/(x^2*(x^2+1+0.8668339838e-5*y^2)^7)+0.4198230391e-5*y^2/(x^4*(x^2+1+0.8668339838e-5*y^2))-0.1399410130e-5*y^2/(x^2*(x^2+1+0.8668339838e-5*y^2))-0.4198230392e-5*y^2/(x^6*(x^2+1+0.8668339838e-5*y^2))+2.156132099*10^(-23)*y^8/x^8-1.828519217*10^(-15)*y^6/(x^6*(x^2+1+0.8668339838e-5*y^2))+9.966267304*10^(-11)*y^4/(x^4*(x^2+1+0.8668339838e-5*y^2))+1.184803427*10^(-17)*y^6/(x^4*(x^2+1+0.8668339838e-5*y^2)^5)-4.804104786*10^(-46)*y^16/(x^2*(x^2+1+0.8668339838e-5*y^2)^7)-1.546256751*10^(-34)*y^12/(x^6*(x^2+1+0.8668339838e-5*y^2)^4)-1.654112392*10^(-23)*y^8/(x^4*(x^2+1+0.8668339838e-5*y^2)^2)+1.361409376*10^(-23)*y^8/(x^6*(x^2+1+0.8668339838e-5*y^2))-2.173253771*10^(-18)*y^6/(x^4*(x^2+1+0.8668339838e-5*y^2))+1.548967341*10^(-33)*y^12/(x^4*(x^2+1+0.8668339838e-5*y^2)^5)-9.04948151*10^(-30)*y^10/(x^6*(x^2+1+0.8668339838e-5*y^2)^4)-1.185404100*10^(-28)*y^10/(x^4*(x^2+1+0.8668339838e-5*y^2)^4)-1.155916435*10^(-28)*y^10/(x^2*(x^2+1+0.8668339838e-5*y^2)^5)+2.004893778*10^(-16)*y^6/(x^4*(x^2+1+0.8668339838e-5*y^2)^3)-4.261515214*10^(-22)*y^8/(x^4*(x^2+1+0.8668339838e-5*y^2)^3)-5.528011082*10^(-12)*y^4/(x^6*(x^2+1+0.8668339838e-5*y^2)^3)-2.018622531*10^(-10)*y^4/(x^6*(x^2+1+0.8668339838e-5*y^2))-4.961084548*10^(-17)*y^6/(x^6*(x^2+1+0.8668339838e-5*y^2)^2)+1.054605549*10^(-17)*y^6/(x^4*(x^2+1+0.8668339838e-5*y^2)^2)-1.209855743*10^(-12)*y^4/(x^6*(x^2+1+0.8668339838e-5*y^2)^2)+9.370950850*10^(-17)*y^6/x^8+1.299951384*10^(-12)*y^4/(x^4*(x^2+1+0.8668339838e-5*y^2)^2)+2.090829576*10^(-18)*y^6/(x^4*(x^2+1+0.8668339838e-5*y^2)^4)-9.760926200*10^(-23)*y^8/(x^8*(x^2+1+0.8668339838e-5*y^2)^2)-1.380589551*10^(-28)*y^10/(x^8*(x^2+1+0.8668339838e-5*y^2)^2)+5.348064193*10^(-8)*y^2/x^8+6.969431921*10^(-19)*y^6/(x^6*(x^2+1+0.8668339838e-5*y^2)^4)-6.650553453*10^(-41)*y^14/(x^4*(x^2+1+0.8668339838e-5*y^2)^5)+0.1399410130e-5*y^2/(x^8*(x^2+1+0.8668339838e-5*y^2))+4.603340601*10^(-35)*y^12/(x^2*(x^2+1+0.8668339838e-5*y^2)^5)-1.080989355*10^(-27)*y^10/(x^4*(x^2+1+0.8668339838e-5*y^2)^5)+1.726252725*10^(-34)*y^12/(x^4*(x^2+1+0.8668339838e-5*y^2)^4)+6.265658040*10^(-35)*y^12/(x^6*(x^2+1+0.8668339838e-5*y^2)^3)+1.031803251*10^(-10)*y^4/(x^8*(x^2+1+0.8668339838e-5*y^2))-1.265918665*10^(-34)*y^12/(x^8*(x^2+1+0.8668339838e-5*y^2)^2)+9.189513286*10^(-25)*y^8/(x^2*(x^2+1+0.8668339838e-5*y^2)^3)+3.850128166*10^(-29)*y^10/(x^6*(x^2+1+0.8668339838e-5*y^2)^2)-3.186312965*10^(-29)*y^10/(x^4*(x^2+1+0.8668339838e-5*y^2)^3)-9.009564044*10^(-14)*y^4/(x^8*(x^2+1+0.8668339838e-5*y^2)^2)-4.257951866*10^(-34)*y^12/(x^2*(x^2+1+0.8668339838e-5*y^2)^7)-1.048302262*10^(-23)*y^8/(x^2*(x^2+1+0.8668339838e-5*y^2)^7)-2.151514446*10^(-16)*y^6/(x^6*(x^2+1+0.8668339838e-5*y^2)^3)+2.992749054*10^(-40)*y^14/(x^6*(x^2+1+0.8668339838e-5*y^2)^4)+3.378713013*10^(-27)*y^10/(x^6*(x^2+1+0.8668339838e-5*y^2)^3)+1.109023300*10^(-17)*y^6/(x^2*(x^2+1+0.8668339838e-5*y^2)^3)-1.184803426*10^(-17)*y^6/(x^2*(x^2+1+0.8668339838e-5*y^2)^5)+2.720007655*10^(-23)*y^8/(x^2*(x^2+1+0.8668339838e-5*y^2)^5)-7.077039402*10^(-29)*y^10/(x^2*(x^2+1+0.8668339838e-5*y^2)^5)+5.275258127*10^(-17)*y^6/(x^4*(x^2+1+0.8668339838e-5*y^2)^3)-1.621936746*10^(-22)*y^8/(x^4*(x^2+1+0.8668339838e-5*y^2)^3)+0.3048396590e-5*y^2/(x^4*(x^2+1+0.8668339838e-5*y^2))-0.1524198295e-5*y^2/(x^2*(x^2+1+0.8668339838e-5*y^2))-0.1524198295e-5*y^2/(x^6*(x^2+1+0.8668339838e-5*y^2))+1.069612839*10^(-7)*y^2/x^6+5.088594195*10^(-18)*y^6/x^6-2.987820419*10^(-12)*y^4/x^6-2.935160763*10^(-13)*y^4/x^4-1.069612838*10^(-7)*y^2/x^4+5.371293872*10^(-18)*y^6/(x^4*(x^2+1+0.8668339838e-5*y^2))+3.282055473*10^(-12)*y^4/(x^4*(x^2+1+0.8668339838e-5*y^2)^3)+5.791862600*10^(-13)*y^4/(x^4*(x^2+1+0.8668339838e-5*y^2)^2)-3.282055473*10^(-12)*y^4/(x^2*(x^2+1+0.8668339838e-5*y^2)^3)+1.930620867*10^(-13)*y^4/(x^6*(x^2+1+0.8668339838e-5*y^2)^2)-3.484715961*10^(-18)*y^6/(x^2*(x^2+1+0.8668339838e-5*y^2)^5)+1.042309231*10^(-17)*y^6/(x^2*(x^2+1+0.8668339838e-5*y^2)^3)-1.388637563*10^(-23)*y^8/(x^6*(x^2+1+0.8668339838e-5*y^2))+6.137787468*10^(-35)*y^12/(x^2*(x^2+1+0.8668339838e-5*y^2)^5)+7.080695477*10^(-30)*y^10/(x^4*(x^2+1+0.8668339838e-5*y^2)^3)-1.837902657*10^(-23)*y^8/(x^4*(x^2+1+0.8668339838e-5*y^2)^2)-2.445967302*10^(-13)*y^4/(x^2*(x^2+1+0.8668339838e-5*y^2))-4.901073753*10^(-24)*y^8/(x^2*(x^2+1+0.8668339838e-5*y^2)^3)-3.186312965*10^(-29)*y^10/(x^6*(x^2+1+0.8668339838e-5*y^2)^2)+2.583978611*10^(-18)*y^6/(x^6*(x^2+1+0.8668339838e-5*y^2)^2)+5.556889574*10^(-18)*y^6/(x^4*(x^2+1+0.8668339838e-5*y^2)^2)+4.216291924*10^(-24)*y^8/(x^6*(x^2+1+0.8668339838e-5*y^2)^2)-5.331502240*10^(-16)*y^6/(x^6*(x^2+1+0.8668339838e-5*y^2))+6.326195740*10^(-11)*y^4/(x^4*(x^2+1+0.8668339838e-5*y^2))-6.440499442*10^(-11)*y^4/(x^6*(x^2+1+0.8668339838e-5*y^2))+-1.447965650*10^(-12)*y^4/(x^2*(x^2+1+0.8668339830e-5*y^2)^3)+7.534765512*10^(-18)*y^6/(x^2*(x^2+1+0.8668339830e-5*y^2)^3)-9.802147473*10^(-24)*y^8/(x^2*(x^2+1+0.8668339830e-5*y^2)^3)+2.546798240*10^(-11)*y^4/(x^4*(x^2+1+0.8668339830e-5*y^2))+0.1363756369e-5*y^2/(x^4*(x^2+1+0.8668339830e-5*y^2))-0.1363756369e-5*y^2/(x^2*(x^2+1+0.8668339830e-5*y^2))-1.130798707*10^(-18)*y^6/(x^4*(x^2+1+0.8668339830e-5*y^2))+7.827095355*10^(-13)*y^4/(x^2*(x^2+1+0.8668339830e-5*y^2))+1.174064303*10^(-12)*y^4/x^4+1.604419258*10^(-7)*y^2/x^4+(-0.1203314443e-5*y^2/x+3.130838142*10^(-12)*y^4/x)/(x*(x^2+1+0.8668339830e-5*y^2)))^2)',[0.7,1.3],[0,100]); hold on ezplot('(-x^2+1+0.4334169915e-4*y^2)^2+7.578718285*10^(-14)*x^2 = .7958731346/(y^2*(1+-0.2674032091e-5*y^2/x^6+0.3498525326e-4*y^2/(x^8*(x^2+1+0.4334169919e-4*y^2))+1.203165052*10^(-25)*y^10/(x^6*(x^2+1+0.4334169919e-4*y^2)^2)+6.591284685*10^(-15)*y^6/(x^4*(x^2+1+0.4334169919e-4*y^2)^2)-9.957228014*10^(-26)*y^10/(x^4*(x^2+1+0.4334169919e-4*y^2)^3)+5.239410172*10^(-23)*y^10/(x^6*(x^2+1+0.4334169919e-4*y^2)^3)-1.327128747*10^(-18)*y^8/(x^4*(x^2+1+0.4334169919e-4*y^2)^3)-2.885193661*10^(-14)*y^6/(x^8*(x^2+1+0.4334169919e-4*y^2)^2)+4.261750648*10^(-23)*y^10/(x^2*(x^2+1+0.4334169919e-4*y^2)^7)-1.033820245*10^(-20)*y^8/(x^4*(x^2+1+0.4334169919e-4*y^2)^2)-1.806119431*10^(-24)*y^10/(x^2*(x^2+1+0.4334169919e-4*y^2)^5)+5.743445804*10^(-22)*y^8/(x^2*(x^2+1+0.4334169919e-4*y^2)^3)-7.447217702*10^(-13)*y^6/(x^6*(x^2+1+0.4334169919e-4*y^2)^3)-2.918546203*10^(-17)*y^8/(x^6*(x^2+1+0.4334169919e-4*y^2)^3)+1.238901614*10^(-28)*y^12/(x^4*(x^2+1+0.4334169919e-4*y^2)^5)+1.069586344*10^(-18)*y^8/(x^2*(x^2+1+0.4334169919e-4*y^2)^5)-4.094930710*10^(-18)*y^8/(x^2*(x^2+1+0.4334169919e-4*y^2)^7)+1.461233923*10^(-8)*y^4/(x^8*(x^2+1+0.4334169919e-4*y^2))+6.910013852*10^(-9)*y^4/(x^4*(x^2+1+0.4334169919e-4*y^2)^3)+0.133701604e-5*y^2/x^4+6.819804704*10^(-13)*y^6/(x^4*(x^2+1+0.4334169919e-4*y^2)^3)-2.352267844*10^(-19)*y^8/(x^6*(x^2+1+0.4334169919e-4*y^2)^2)+0.1049557598e-3*y^2/(x^4*(x^2+1+0.4334169919e-4*y^2))-0.9203127131e-4*y^2/(x^2*(x^2+1+0.4334169919e-4*y^2))-0.1049557598e-3*y^2/(x^6*(x^2+1+0.4334169919e-4*y^2))+0.1337016048e-5*y^2/x^8+5.750831903*10^(-14)*y^6