dsv prak 5

Aufgabe 1:

function z=addprimes(s,e)
z=0;
for i=s:1:e
   if isprime(i) == 1
       z = z + i;
   end
end


Aufgabe 2:

function zero_diag(A)
dim = size(A);
if dim(1) <= dim(2)
    for i=1:1:dim(1)
        A(i,i) = 0;
    end
elseif dim(2) <= dim(1)
   for i=1:1:dim(1)
        A(i,i) = 0;
    end
end
A


Aufgabe 3:

function z=fibo(n)
z = zeros(1,n);
z(1) = 1;
z(2) = 1;
for i=3:1:n
    z(i) = z(i-1) + z(i-2);
end
z(n)


Aufgabe 4:

function r=factorial(n)
z=1;
for i=1:1:n
  z = z*i;
end
r=z;


Aufgabe 5:

function z3=auf5()
n=input('Geben sie ein n > 2 ein:');
x=input('Geben sie ein x ein:');
z2=zeros(1,n);
for k=2:1:n
   z2(k) = ((x/factorial(k)) + 3^k);
end

z3 = sum(z2)


Aufgabe 6:

function auf6()
a=input('Geben sie a ein:');
b=input('Geben sie b ein:');
c=input('Geben sie c ein:');
d=input('Geben sie d ein:');
e=input('Geben sie e ein:');
f=input('Geben sie f ein:');

%ax+by=e
%cx+dy=f

%daraus folgt:

%y=(e-ax)/b
%cx+d((e-ax)/b)=f
%==> x = (f-((d*e)/b))/(c-((da)/b))

%x = (f-((d*e)/b))/(c-((d*a)/b))
%y = (e-a*x)/b

A=[a,b;c,d]
B=[e;f]
A\B


Aufgabe 7:

function z1=bodediag(n)
y=logspace(0,3,200)
z=[2 5 7];
n=[3 1 4 9];
g=tf(z,n)
w=2*pi*y
bode(g,w)