Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- %Zadanie 2
- imie='Kamil';
- nazwisko='Sobolewski';
- A=[0,1; -length(imie), -length(nazwisko)];
- B=[0;1];
- C=[1,0;0,1];
- h=1;
- t=sym('t');
- Aplus=expm(h*A);
- Bplus=double( int( expm( h*A )*B, t, 0, h ) );
- Cplus=C;
- figure();
- sim('model2');
- i=1:h:11;
- hold on;
- plot(i, x_d_2(:,2), 'o');
- plot(x_c_2.Time, x_c_2.Data(:,2));
- hold off;
- %Zadanie 3
- K_0=1000; %kapitał początkowy
- n=5; %czas trwania inwestycji przy założeniu że jest wielkrotnością okresów kapitalizacji
- r=0.05; %nominalna stopa procentowa
- figure();
- for m=1:12
- hold on;
- K=K_0*(1+r/m)^(m*n);
- plot(m, K, 'o');
- end;
- hold off;
- xlabel('Liczba kapitalizacji odsetek w roku(m)');
- ylabel('Kapitał końcowy (K)');
- %Zadanie 4
- liczba_iteracji=15; %maksymalna liczba iteracji
- epsilon_1=10^-6; %dokładność z jaką funkcja jest bliska zeru w wyznaczonym punkcie
- epsilon_2=10^-6; %dokładność określająca minimalną długość kroku w kolejnej iteracji
- x0=-4; %punkt startowy
- f=@(x)x^3-2*x^2-11*x+12; %funkcja f(x)
- g=@(x)3*x^2-4*x-11; %pochodna funkcji f(x)
- i=0; %licznik iteracji
- while i<liczba_iteracji
- i=i+1;
- x=x0-(f(x0)/g(x0));
- if(abs(f(x))<epsilon_1)
- break;
- end;
- if(abs(x-x0)<epsilon_2)
- break;
- end;
- x0=x;
- end;
- disp(x); %wypisuje znalezione miejsce zerowe na ekran
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement