临风暖阳的博客

package mainimport (“fmt”)func main() {//共买100只，所以cock（鸡翁）最多20只，hen(鸡母)最多33只，chicken（鸡雏）最多100只for cock := 0; cock < 20; cock++ {for hen := 0; hen < 33; hen++ {for chicken := 0; chicken <= 100; chicken += 3 {if (5cock+3hen+chicken/3 == 100

clear all;load carbigX=[Acceleration Displacement Horsepower MPG Weight];X=X(all(~isnan(X),2)????;%估计因素负荷矩阵[Lambda,Psi,T,stats,F]=factoran(X,2,‘scores’,‘regression’);Lambda %输出因素负荷矩阵inv(T’T) %F的相关矩阵LambdaLambda’+diag(Psi)%X的相关矩阵Lambdainv(T) %未经旋转的

clear all;strength=[82 86 79 83 84 85 86 87 74 82 78 75 76 77 79 79 77 78 82 79];alloy={‘st’,‘st’,‘st’,‘st’,‘st’,‘st’,‘st’,‘st’,…‘all’,‘all’,‘all’,‘all’,‘all’,‘all’,…‘al2’,‘al2’,‘al2’,‘al2’,‘al2’,‘al2’};[p,a,s]=anova1(strength,alloy);[c,m,h,nums]=mul

3.5765174983199888E913 5.786926873629978E913 9.363444371949958E913 1.5150371245579919E914 2.451381561752986E914 3.966418686310971E914 6.4178002480639496E914 1.0384218934374918E915 1.6802019182438859E915 2.718623811681376E915 4.398825729925255E915 7.1174495

odefun=inline(‘y-2*t/y’,‘t’,‘y’);[t,y]=ode45(odefun,[0,4],1);[t,y]ans =0 1.00000.0502 1.04900.1005 1.09590.1507 1.1408…… ……3.8507 2.95033.9005 2.96723.9502 2.98394.0000 3...