www.isutp.fora.pl

energy69 · Gość

hey! ma ktos zadania ze scilaba z rozwiazaniem?? ja mam ale cos mi niezbyt dobrze wychodza:( a juz niedlugo zaliczenie:(

Stara

daj maila to podesle.

alastor

Mógłbym też prosić Razz

lejcz

ja tez chce:P znasz mojego maila kama:)

fordonutp

za bardzo zaniedlugo Very Happy

ja mam we wtorek.musze na 3.5 zaliczyc . tak mi powiedzial :/

Stara

zad. 4

xbasc() // wyczyszczenie okna graficznego

a=linspace(-25, 25, 100) // zdefiniowanie x'ów

// funkcja 1
y1 = 42 - 20*a - 2*a^2
plot2d(a, y1) // wykreślenie funkcji 1

// funkcja 2
y2 = 7 - 2*a
plot2d(a, y2)

//42 - 20*a - 2a^2=7 - 2*a
y3=-18*a-2*a^2+35

a1=-2
b1=-18
c1=35

delta1 = b1^2 - 4*a1*c1

printf("delta = %f\n", delta1)

// obliczenie 1 punktu przecięcia
x1 = (18 - sqrt(delta1)) / (2 * a1)
// obliczenie 2 punktu przecięcia
x2 = (18 + sqrt(delta1)) / (2 * a1)

// wstawic obliczanie y'ów

dopiszecie sobie y dla x1 i y dla x2 i bedzie całe pieknie

neo · Gość

kto sie jak poprawic blad w 3 zadaniu? odnosnie macierzy B. Bo po jej wywolaniu wychodza liczbyod 1 do 20. A powinny byc od 0 do 3.8? Da rade ktos poprawic?

fordonutp

1
x=linspace (0,10,21) ;
y = 2*x^3-32*x^2+188*x-88 ;
minimum = y(1) ;
maksimum = y(1) ;
for i =2:size ( y , "*" ) do
if y(i) > maksimum then
maksimum = y(i);
end
end
minimum
maksimum

fordonutp

2
a=linspace(5,10,25);
f=2*a.*sin(a.*0.7)-4*sin(a.*0.7)-1;
c=abs(f);
d=min(c)

fordonutp

3
n=10; //wymiar macierzy glownej

function y=l(x,a,b)
y=a*x+b;
endfunction

A=zeros(n,n);
B=zeros(n,1);

for i=1:1:n
B(i,1)=i;
end

for i=2:1:n-1
A(i,i-1)=1;
A(i,i)=-2;
A(i,i+1)=1;
end

A(n,n)=1;
A(1,1)=1;

X1=A/(B');
X2=A^(-1)*B

R1=B-A*X1
R2=B-A*X2
//X

for i=1:1:n
printf("f%d(x)=%.35f*x+%f\n",i,X(i),B(i))
end

//function y=f(x)
//y=x^2
//endfunction

fordonutp

4
a=linspace(0,10,51);
f=2*a^3-32*a^2+118*a-88;
x=min([f])
y=max([f])

fordonutp

5b
n=10; //ilość liczb do wylosowania

Va=grand(1,n,'uin',-10,10)
Vb=Va>2

fordonutp

5 custom
n=10; //ilość liczb do wylosowania

for i=1:1:n
Va(i)=round(rand()*20)-10;
//rand zwraca liczbę od 0 do 1, trzeba przemnożyć przez 20
//żeby uzyskać amplitudę 20 (od 0 do 20) i odjąć 10 żeby
//przesunąć wszystko do przedziału od -10 do 10
if (Va(i)>=2) then
Vb(i)="T";
else
Vb(i)="F";
end

end
Vb'

fordonutp

7
n=6; //wymiar macierzy - potrzebny do ładnego drukowania
A=[5 3 1 0 0 0; //macierz główna
1 1 0 0 -3 -8;
2 0 1 -1 -5 1;
0 -1 -3 0 -5 -3;
0 1 -3 0 -1 +5;
-2 0 0 0 0 7];

B=[7 -5 0 -5 0 5]; //macierz wyrazów wolnych

X=(A^(-1)*B')'; //macierz rozwiązań
E=(A*X')-B'; //macierz błędów
printf("Rozwiązanie\n");
for i=1:1:n
printf("x%d=%f Błąd: xe1=%.20f\n",i,X(i),E(i));
end

fordonutp

8
//pole pod okręgiem danym parametrycznie
//wzory na podstawie [link widoczny dla zalogowanych]
function x=x(t) //pierwsze równanie parametryczne
x=2*cos(t)+3
endfunction
function y=y(t) //drugie równanie parametryczne
y=2*sin(t)+1
endfunction
t1=0; //minimalny parametr
t2=2*%pi; //maksymalny parametr
//koniec danych z zadania
//-------------------------------------------

ymin=y(1.5*%pi); //"dołek" okręgu
ymax=y(0.5*%pi); //"góra" okręgu
R=abs(ymax-ymin)/2; //promień okręgu

d=0.000000000001; // precyzja
//szukanie pierwszego przecięcia z osią OX metodą bisekcji
a=0; //początek przeszukiwanego przedziału
b=1.5*%pi; //koniec przeszukiwanego przedziału
while (abs(a-b)>d)
t0a=(a+b)/2;
if (y(a)*y(t0a)<0) then
b=t0a; else a=t0a; end
end

//szukanie drugiego przecięcia z osią OX metodą bisekcji
a=1.5*%pi; //początek przeszukiwanego przedziału
b=t2; //koniec przeszukiwanego przedziału
while (abs(a-b)>d)
t0b=(a+b)/2;
if (y(a)*y(t0b)<0) then
b=t0b; else a=t0b; end
end

// ! ! ! TU SIĘ DZIEJE CAŁA MAGIA ! ! !
// czytać linka z drugiej linii Razz

c=abs(x(t0a)-x(t0b)); //długość cięciwy
alpha=2/sin(0.5*c/R); //kąt między promieniami przecinającymi oś OX
S=0.5*R^2*(alpha-sin(alpha)); //pole powierzchni pod osią OX
// ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! !

printf("Pole wynosi %f\n",S);