Kjhdffddsfghjk

1) Qual a função da subrotina SEARCH no M1
A subrotina procura por uma raiz de uma função através de um método numérico, chamado método da bisseção.

2) Se Vx varia com X na prática, será necessário usar uma integral em X. Simplifique para tiras estreitas de 0,5m e explique como obter a força cisalhante em cada tira.

A função que rege o efeito que queremos escolher é a seguinte: f* d²fdn² + 2* d³fdn³ (I)

Com as seguintes condições Iniciais: f(0) = 0; df0dn=0 df∞dn=1
A solução numérica é constituída dividindo vetorialmente a velocidade, no qual cada ponto da malha representa sua cordenada (x,y) e apresenta uma velocidade V, bem como suas componentes Vx, Vy, sendo assim possível traçar um perfil de velocidades de escoamento com base na equação acima.

Para solucionar o problema, utilizaremos a rotina RKF45 (Runge-Kutta-Fehlberg Method), fornecida junto com o código do programa. A rotina utiliza o método de Runge-kutta de quarta e quinta ordem para a solução de equações diferenciais. Para tal, é necessário transformarmos a equação ( I ) em equações diferenciais de primeira ordem. Modificando as variáveis, alcançamos tais funções:

g1 = f g2 = dfdn g3 = d²fdn²

Sendo o fluido à 20ºC, temos: µ = 1,002 x 10-3 Kg/m.s (viscosidade do fluido) p = 1,465 Kg/m³ (densidade do fluido)
ᵘ = 1,300 x 10-3 (viscosidade cinemática)

Inserindo os dados no programa e pegando o resultado final (que acontece após 20 iterações do método RKF45), temos a seguinte tabela-resultado:
Para tiras de 0,5m: N | G1 | G2 | G3 | 0,00 | 0,00000 | 0,00000 | 0,33615 | 0,50 | 0,04200 | 0,16793 | 0,33498 | 1,00 | 0,16761 | 0,33382 | 0,32688 | 1,50 | 0,37466 | 0,49265 | 0,30596 | 2,00 | 0,65787 | 0,63713 | 0,26932 | 2,50 | 1,00812 | 0,75966 | 0,21897 | 3,00 | 1,41297 | 0,85497 | 0,16195 | 3,50 | 1,85836 | 0,92208 | 0,10767 | 4,00 | 2,33090 | 0,96441 | 0,06380 | 4,50 | 2,81967 | 0,98817 | 0,03352 | 5,00 | 3,31709 | 1,00000 | 0,01556 | 5,50 |