Programação Basica

1) Script método Newton-Raphison.

//script: metodo de Newton-Raphson
// f(x)= x * x – 2
N = 10; x0 = 1.0; delta = 10^(-5); xn = x0; for n=1:N xn1 = xn - (xn * xn - 2)/(2 * xn); if abs((xn1-xn)/xn1) < delta then printf('Valor da raiz = %10.7f',xn1) return end xn = xn1; end 2) Funções do método de Runge-Kutta.
Função Principal function [XY] = rk4(a, b, h, y0)
// Resolucao de equacoes diferenciais ordinarias por Runge-Kutta 4a. ordem
// Entrada : [a,b] - intervalo de integracao
// h - passo da integracao
// y0 - condicao inicial em x0 X(1) = a
Y(1) = y0
Exato(1) = f2(X(1), Y(1)) n = (b-a)/h for k=1:n xk = X(k) yk = Y(k) hf1 = h * f(xk, yk) hf2 = h * f(xk + h/2, yk + hf1/2) hf3 = h * f(xk + h/2, yk + hf2/2) hf4 = h * f(xk + h, yk + hf3)
Y(k+1) = Y(k) + (hf1 + 2*hf2 + 2*hf3 + hf4)/6
X(k+1) = X(k) + h
Exato(k+1) = f2(X(k+1), Y(k+1)) end XY = [X Y Exato] endfunction Função Auxiliar 1 function [fxy] = f(x,y)
// funcao exemplo fxy = (x - y)/2 endfunction Função Auxiliar 2 function [fexato] = f2(x,y)
// funcao solucao fexato = 3 * exp(-x/2) + x - 2 endfunction 3) Gráfico de funções exatas e aproximando da equação dy/dx = (x-y)/2.

function[z]=fx(x,y) z=(x-y)/2 endfunction

x=[-2:0.1:2]; y=x; w=feval(x,y,fx); plot3d(x,x,w) 4) Gráfico Tri-dimensional.
O comando plot3d() permite traçar gráficos de superfícies, z = f(x, y)
Na notação Scilab, as variáveis independentes x e y são vetores de dimensões n1 e n2, respectivamente, e a variável dependente z é uma matriz de dimensão n1 × n2. Então, por essa definição, o elemento z(i,j) é o valor da superfície no ponto (x(i), y(j)).
Além da função plot3d(), e de suas variações, Scilab implementa outras funções que permitem traçar gráficos tri-dimensionais. Dentre elas, destacamos:
• fplot3d - que permite traçar gráficos de superfícies definidas por funções, como no script mostrado no exemplo: