Run the Two Matlab Files Below and Voilà

Total Page:16

File Type:pdf, Size:1020Kb

Run the Two Matlab Files Below and Voilà

Run the two Matlab files below and voilà:

2

1 . 8

1 . 6

1 . 4

1 . 2

1

0 . 8

0 . 6

0 . 4

0 . 2

0 0 1 2 3 4 5 6 7

Here are two Matlab M-files that together write my name.

M-file “signature.m” % Signature hold off [T,Y] = ode45(@sig2,[0 10],[1 .8 2.77 .815 4.05407 1.3562 .92 1 . 75407 .770915 2 1.02 1.71 .126 3.5 .2 4.37 .425 6.45407 .966215 4.6 . 5901 5 .6 1.3 .13]); %plot (Y(:,1), Y(:,2)) %plot (Y(:,3), Y(:,4)) %plot (Y(:,5), Y(:,6)) plot (Y(:,7), Y(:,8)) hold on plot (Y(:,9), Y(:,10)) plot (Y(:,25), Y(:,26)) plot (Y(:,11), Y(:,12)) plot (Y(:,13), Y(:,14)) plot (Y(:,15), Y(:,16)) plot (Y(:,17), Y(:,18)) plot (Y(:,19), Y(:,20)) plot (Y(:,21), Y(:,22)) plot (Y(:,23), Y(:,24)) M-file “sig2.m” function dy = sig2(t,y) dy = ones(26,1); dy(1) = y(1)^2 * y(2) - y(1) * y(2)^3; dy(2) = y(1)^2 * y(2) - y(1) * y(2)^2; dy(3) = 0.0178666*y(3)^3 * y(4) - 0.1072*y(3)^2 * y(4)^3; dy(4) = 0.0893333*y(3) * y(4)^2 - 0.536*y(4)^3; dy(5) = .0266666*y(5)^3 * y(6) - .16*y(5)^2 * y(6)^3; dy(6) = 0.13333333*y(5) * y(6)^2 - 0.8*y(6)^3; dy(7) = 1.05*y(7)^3 * y(8) - 1.05*y(7)^2 * y(8)^3; dy(8) = 0.75*y(7) * y(8)^2 - 0.75*y(8)^3; dy(9) = .016*y(5)^2.6 * y(6) - .096*y(5)^1.6 * y(6)^3; dy(10) = 0.106666*y(5) * y(6)^2 - 0.64*y(6)^3; dy(11) = 0.15*y(11)^2 * y(12)^3 - 0.075*y(11)^3 * y(12); dy(12) = 1.5*y(11)^0 * y(12)^3 - 0.75*y(11)*y(12)^2; dy(13) = 0.0625*y(13)^3 * y(14) - 0.1111111*y(13)^2 * y(14)^3; dy(14) = 0.625*y(13) * y(14)^2 - 1.666666*y(14)^3; dy(15) = y(1)^2 * y(2) - y(1) * y(2)^3; dy(16) = 0.5*y(1)^2 * y(2) - 0.5*y(1) * y(2)^2; dy(17) = 0.0178666*y(3)^3 * y(4) - 0.1072*y(3)^2 * y(4)^3; dy(18) = 0.0893333*y(3) * y(4)^2 - 0.536*y(4)^3; dy(19) = 0.02666666*y(5)^3 * y(6) - 0.16*y(5)^2 * y(6)^3; dy(20) = 0.1333333*y(5) * y(6)^2 - 0.8*y(6)^3; dy(21) = 0.0787649*y(21)^3 * y(22) - 1.006441*y(21)^2 * y(22)^3; dy(22) = 0.3623188*y(21) * y(22)^2 - 2.777777*y(22)^3; dy(23) = 0.132; dy(24) = 0; dy(25) = .1*(y(1)^2 * y(2) - y(1) * y(2)^3); dy(26) = .1*(0.5*y(1)^2 * y(2) - 0.5*y(1) * y(2)^2); end

Recommended publications