Respondida
2D information to 1D colorplot
Simple example x = 0:10; y = x.^2; temp = rand(1,10); % temperature cm = jet(10); % colormap ind = 1+rou...

más de 6 años hace | 0

Respondida
Quiver with a system of ODEs.
For vector field u = x'; v = y'; To plot the curve [T,Y] = ode45 ... plot(Y(:,1),Y(:,2)) % (x,y) I suppose x and y depen...

más de 6 años hace | 0

Respondida
Solving nonLinear equation and Linearized Eqaution using ODE45
YOu have to express , and . Solve the system function du = nonlin(t,y) fi = y(1); psi = y(2); theta = y(3); ...

más de 6 años hace | 0

| aceptada

Respondida
How to plot magnitude and directions between 2 points?
try this (not tested) data = load ('Input.txt'); [row, col] = size (data); Y2 = data (:,1); X2 = data (:,2); Y1 = data (:,3...

más de 6 años hace | 0

| aceptada

Respondida
surface fixed point value Z
I want to know the value of Z from the surfc function Here is it

más de 6 años hace | 0

Respondida
Where are the bugs for this ODE finite difference problem that solve using Newton Raphson method?
Here is my attempt clc,clear N = 10; % 10 Set parameters L = 16*0.3048; %meter b = 19.625 * 0.0254; d = 1.625 * 0.0254; ...

más de 6 años hace | 1

| aceptada

Respondida
I am receiving the error message Index exceeds the number of array elements (1) when trying to use ode45.
Try @(V,T)

más de 6 años hace | 0

| aceptada

Respondida
I have a 14 equation system of second order differential equations which I am trying to solve using dsolve. The program just keeps running without finishing ?
alpha and E should be numerical function main m = 15000; k = 150000; n = 14; A = zeros(n); A(n,n) = -k/m; A(n,n-1) = k/m;...

más de 6 años hace | 0

Respondida
Euler 3D rotation between two vectors
Here is rotation using Rodriguez matrix Here is rotation about Y axis and Z. Here is the difference between them See ...

más de 6 años hace | 0

| aceptada

Respondida
Simultaneous differential equations - derivative of y^3 wrt t
Just get dydt(1)=(-3*D*Cs/rho/r0^2)*y(1)*(1-y(2)); dydt(1)=nthroots(dydt(1),3);

más de 6 años hace | 0

Respondida
Numerically derive a continous, non-symbolic function
Derivative is (if you have numerical data) dy = (y(i)-y(i-1)) / (t(i)-t(i-1)); Maybe diff (if you have a function) ? syms x ...

más de 6 años hace | 0

| aceptada

Respondida
3D Point Cloud - Gaussian Curvature
one way clc,clear % generate some data r = 3; t = linspace(0,10)'; x = r*cos(t); y = r*sin(t); z = sin(1*t); t1 = ...

más de 6 años hace | 0

| aceptada

Respondida
Runge Kutta for system of eqs
I wrote function in the way similar to ode45 y_dot = @(x,dx)[dx; -4*dx-5*x]; %insert function to be solved for i = 1:n k1...

más de 6 años hace | 0

Respondida
The problem of finding intersection point of triangle mesh
To check every edge: for j = 1:size(face,2) i1 = 1; for i2 = [2 3 1] % checking edges [1 2], [2 3], [3 1...

más de 6 años hace | 1

| aceptada

Respondida
How can i solve a system which contains matrices
Here is the best way i know: A = [5 3 1 1 1 1 3 1 2]; B = [c1(:) c2(:) c3(:)]'; C = inv(A) * B;

más de 6 años hace | 1

| aceptada

Respondida
Using a for-loop to score values in different intervals
Correct form if 0 <= x(i) && x(i) < 1

más de 6 años hace | 0

| aceptada

Respondida
Intersection between line and circle when line ends inside circle (using geom2d)
Someone should copy their comment to the answers section so it can be accepted. It should be me % Position of line: x1 = [0...

más de 6 años hace | 1

| aceptada

Respondida
Removing points from meshgrid (mask)
Try this simple script [X,Y] = meshgrid(1:20); n = [1 -1]; n = n/norm(n); % normal of a line dist = -7; ...

más de 6 años hace | 0

Respondida
How can I find the surface normal of a 8 by 3 matrix?
Use griddata to create a surface Use surfnorm to find normal vectors [nx,ny,nz] = surfnorm(x,y,z); surf(x,y,z) hold on quiv...

más de 6 años hace | 0

Respondida
How do I increase the circle size and animate a stick figure?
Use for loop and pause t = linspace(0,2*pi); x = cos(t); y = sin(t); axis([-10 10 -10 10]) hold on for i = 1:10 h = p...

más de 6 años hace | 0

| aceptada

Respondida
Finding specific data point x for given y
I need to find the CB and CS values at t=5 Use interp1 to find i need to find the t and CB values when CS = 0.001 Use fsolve ...

más de 6 años hace | 0

Respondida
Surface between 3D lines
You are concatenating in a wrong way XX = [D1 D1 D1]; YY = [F1 F2 F3]; ZZ = [z1 z2 z3]; surf(XX,ZZ,YY,'EdgeColor','none')

más de 6 años hace | 2

| aceptada

Respondida
two body problem using ode45
- equation of equilibrium Is this condition fullfield?

más de 6 años hace | 0

Respondida
Numerically Integrating the differential equation below?
Any guidance on it is much appreciated! Helpfull page: ODE45 I do not know what I can do if I no tf. Try tf = 5. This should ...

más de 6 años hace | 0

Respondida
Streamline plot not complete
Works ok for me N = 15; xx = linspace(min(x(:)),max(x(:)),10); yy = linspace(min(y(:)),max(y(:)),10); [startx,starty] = mes...

más de 6 años hace | 1

| aceptada

Respondida
2D Matrix revolution
Use cylinder for each curve separately: z = linspace(-1,1,20); % axis r = sqrt(1-z.^2); ...

más de 6 años hace | 1

Respondida
How do I align the vertices and the adjacency graph of a Buckyball?
Here is a way ind = full(B); % convert to ordinary matrix [rr,cc] = find(ind); % find non-zero elements ix = ...

más de 6 años hace | 0

Respondida
Find the coordinates of a point chosen on a plot
Short example x = linspace(0,2,20)'; % generate some data y = sin(x); plot(x,y,'.-b') xy = ginput(1); ...

más de 6 años hace | 3

| aceptada

Respondida
How do I create a sphere without sphere function that will inflate like a balloon utilizing the getframe function?
Use surf to create a 3D surface rr = linspace(0,0.5,20); tt = linspace(0,2*pi,40); [R,T] = meshgrid(rr,tt); [X,Y] = pol2cart...

más de 6 años hace | 0

Respondida
Differential Algebraic Equations with vectors
Try f = @(t,u) inv(Z)*(J+p*W*pi*Q)*[u(3) u(4)]' + inv(Z)*(V+p*W*pi*Q)*[u(1) u(2)]'; F = @(t,u) [u(3) u(4) ...

más de 6 años hace | 1

| aceptada